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ABSTRACT 


This  thesis  develops,  implements  and  tests  a  Tactical  Decision  Aid 
for  a  Reactive  Target  ASW  Active  Search.  The  model  uses  a  Bayesian 
Filtering  Process  to  fuse  information  from  a  real  world  search  conducted 
by  several  assets  with  information  from  a  Monte  Carlo  Simulation  that 
encompasses  five  hundred  equally- likely  different  possible  initial 
positions  and  behaviors  of  the  real  target.  A  Reactive  Target  Model 
resembles  the  behavior  of  a  target  that  is  always  aware  and  reacts  because 
of  the  presence  and  activity  of  the  searchers.  An  initial  "prior",  or  best 
estimate  of  the  location  of  the  target,  is  updated  using  the  movement  of 
the  simulated  targets,  the  negative  information  conveyed  in  an  unsuccesful 
search  over  a  period  of  time  and  the  positive  information  implied  in  a 
contact  report.  The  search  effort  is  measured  using  a  Fixed  Scan 
Stochastic  Model  that  solves  the  Sonar  Equation  limited  by  noise  and.  As 
a  result  of  updating  the  prior,  a  "posterior"  distribution  is  obtained. 
The  Law  of  Total  Probabilities  is  used  to  render  a  probability  map  of  the 
location  of  the  Target  by  mapping  color  intensities  to  probabilities.  A 
recursive  expression  for  evaluating  a  contact  report  is  also  developed. 
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I .  INTRODUCTION 


A.  BACKGROUND 

The  end  of  the  cold  war  started  a  transition  from  the 
Major  Regional  Conflicts  to  the  Lesser  Intensity  Conflicts 
(Hughes,  1993, pp.l)  as  the  more  likely  scenario  to  be  devel¬ 
oped  in  the  future  for  US  Forces.  Such  a  transition  renewed 
the  interest  and  efforts  to  study  what  is  called  "Littoral 
Warfare".  Other  allied  countries  have  had  the  Littoral  Warfare 
as  their  natural  scope  and  perroanent  frame  for  the  design  and 
development  of  the  capabilities  of  their  naval  forces. 

Above  all.  Littoral  Warfighting  implies  sensible  changes 
in  the  physical  nature  of  the  properties  of  the  threat  and  the 
environment  which  together  with  the  mission  and  the  current 
capabilities  determines  the  imbalance  that  drives  research, 
development  and  acquisition. 

1.  The  environment  in  Littoral  Warfare 

The  changes  in  the  physical  properties  of  the  environment 
*^ne  to  the  transition  from  the  Open  Sea  to  Littoral  or 
Coastal  waters  have  already  been  thoroughly  studied  and  are 
characterized  by  the  distinction  between  deep  and  shallow 
waters,  often  termed  blue  and  green  waters  respectively.  A 
particular  maritime  area  is  considered  a  shallow  water  en- 
vironment  in  terms  of  the  influence  its  boundaries  exert  in 
its  physical  properties  rather  than  in  teirms  of  a  particular 
predetermined  depth. 
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From  the  standpoint  of  the  Search  and  Detection  problem, 
which  remains  crucial  to  Littoral  Warfighting  (Hughes,  1986 ,pp- 
126) ,  the  most  remarkable  differences  between  shallow  or  deep 
water  environments  are  the  intensity  of  noise  and  clutter  and 
the  spatial  and  temporal  variability  in  the  properties  of  the 
ocean  which  become  distinctive  and  determinant  factors 
affecting  the  performance  of  the  sensors  and  the  decision  aids 
associated  with  them.  Several  factors  have  a  direct  impact  in 
increasing  the  intensities  of  clutter  and  noise.  The  most 
important  factors  are  the  wind,  the  biological  content  of  the 
water,  the  proximity  of  the  bottom,  and  human  activity. 

The  wind  increases  both  ambient  noise  and  clutter.  It  has 
been  observed  that  for  the  typical  frequencies  used  in  passive 
sonar,  the  ambient  noise  in  coastal  locations  is  in  general 
more  intense  in  shallow  waters .  The  simplest  reason  is  the 
larger  mass  of  absorption  in  the  deep  ocean.  In  shallow  waters 
the  noise  produced  by  the  wind  (white  caps,  bubbles,  waves, 
etc.)  is  partially  reflected  by  the  bottom  and  the  surface, 
hence  the  dissipation  process  is  slower.  The  increase  in 
clutter  is  due  to  the  roughness  of  the  surface  which  in  turn 
increases  as  the  wave  length  decreases  due  to  the  reduction  of 
depth  (Urick,1975,  pp.  236).  Additionally,  particular  condi¬ 
tions  of  short  waves  typical  of  shallow  waters  like  tide  rips, 
choppy  seas,  etc.,  make  detection  harder  due  to  the  movement 
transmitted  to  the  platforms  where  the  sensors  are  carried. 

Biological  activity  increases  as  depth  decreases  thus 
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leading  to  an  increase  in  both  noise  (crabs,  shrimps, etc. )  and 
clutter  (whales,  schools  of  fish,  zooplankton,  etc. ) . (Urick, 
1975, pp. 193) 

Human  activity  produces  both  noise  and  clutter.  Shipping 
activity  increases  close  to  shore.  A  fishing  vessel,  for 
instance,  can  be  considered  a  source  of  both  noise  and 
clutter. 

Finally,  the  closer  bottom  of  the  ocean  not  only  increas¬ 
es  the  backscattering  due  to  reflection  in  the  case  of  an 
active  sensor,  but  also  prevents  the  formation  of  caustics  or 
convergence  zones  typical  of  deep  waters. 

Meteorological  and  Oceanographical  changes  also  intensify 
when  transitioning  from  the  macroscale  of  the  open  ocean  to 
the  mesoscale  of  coastal  regions.  For  instance,  daily  varia¬ 
tions  in  the  temperature  profile  in  the  water  and  in  the  air 
lead  to  drastic  changes  in  the  propagation  of  electromagnetic 
and  acoustic  energy. 

The  prediction  of  the  performance  of  a  sensor  for  a 
PS’^ticular  target  and  environment  is  of  the  essence  in  the 
searching  problem.  The  more  intense  the  variability  in  the 
conditions  found  in  shallow  or  coastal  waters  the  harder  the 
prediction.  Consequently,  the  uncertainty  in  the  evaluation  of 
the  effectiveness  of  a  given  detection  effort  increases. 

2 .  Changes  in  the  threat 

With  the  end  of  the  cold  war,  the  most  likely  threat 
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became  smaller  but  not  weaker,  at  least  from  a  conventional 
warfare  point  of  view.  Small  fast  patrol  boats  can  strike  fast 
with  remarkable  lethality  and  are  not  easy  to  detect  in 
choppy  seas .  Modern  conventional  submarines  have  the  same  or 
even  better  information  systems  than  older  nuclear  submarines, 
taking  advantage  of  multiple  sensors  and  massive  signal 
processing  based  on  inexpensive  and  commercially  available 
technology.  Except  for  the  snorkel ing  cycle  where  they  become 
noisy  and  visible,  conventional  submarines  on  batteries  are 
almost  undetectable  by  current  passive  sensors .  Their  smaller 
size  makes  it  more  difficult  to  detect  them  with  active 
sensors  due  to  a  reduction  in  target  strength. 

3 .  Sensors  and  Tactical  Decision  Aids (TDA) 

The  threat  and  the  environment  are  determinant  factors  in 
the  decision  to  use  a  particular  type  of  sensor  and  whether  it 
will  be  used  in  a  passive  or  active  fashion.  Regardless  of  its 
type  or  mode,  a  sensor  normally  requires  a  Tactical  Decision 
Aid  for  a  more  intelligent  and  efficient  use. 

In  the  realm  of  the  ASW  Search  Problem,  the  fact  that  a 
diesel  submarine  on  batteries  is  almost  undetectable  with 
passive  acoustic  means,  especially  in  the  noisy  shallow  water 
environment,  drives  the  need  to  use  active  sensors.  This 
conclusion  can  be  extended  to  other  threat -environment 
combinations  in  shallow  waters,  for  instance  small  fast 
surface  combatants  because  they  normally  operate  in  silence 
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unless  they  are  about  to  strike. 

A  sensor  typically  feeds  a  Tactical  Data  System  where 
software  tools  (TDA)  help  decide  the  best  course  of  action 
given  a  particular  profile  of  mission,  environment  and  threat. 
The  essence  of  a  Tactical  Decision  Aid  is  providing  qualita¬ 
tive  and  quantitative  processing  of  the  data,  so  it  becomes 
information.  The  factors  of  the  situation  should  be  ranked, 
simplified  and  presented  in  a  way  such  that  the  determinant 
ones  are  emphasized  in  terms  of  the  mission. 

The  smaller  signal  to  noise  and  target  to  clutter  ratios 
due  to  the  more  noisy  and  cluttered  environment  will  not  only 
make  detection  much  more  difficult  but  .also,  together  with 
strong  environmental  variability,  will  generate  more  uncer¬ 
tainty  about  the  effectiveness  of  the  search.  As  a  conse¬ 
quence,  the  need  for  a  stochastic  TDA  to  keep  adequate  track 
of  the  uncertainties  is  evident.  Figure  1.1  summarizes  the 
problem  structure  leading  to  the  development  of  a  TDA.  The 
chains  of  causes  and  effects  are  displayed  as  a  flow  pattern. 

B.  THE  PROBLEM 

For  any  Navy  that  does  not  have  a  TDA  capable  of  handling 
the  uncertainty,  the  problem  is  apparent  and  the  need  to 
develop  one  is  clear. 

A  large  array  of  TDA's  in  the  U.S.  NAVY  like  VPCAS, 
ASWTDA  and  PACSEARCH  (Wagner,  1989 , pp . II . 2)  responded  to  the 
needs  of  the  cold  war  where  the  most  important  means  of 
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detection  was  the  passive  sonar.  In  a  passive  sonar  search, 
specially  from  airborne  platforms,  the  target  will  usually  not 
be  aware  of  the  presence  of  the  searchers.  As  a  consequence 
those  TDA's  do  not  consider  any  reaction  of  the  target  caused 
by  the  presence  or  activities  of  the  searcher. 

In  an  active  search,  the  target  will  almost  invariably 
know  of  the  presence,  location  and  activities  of  the  searcher 
before  detection.  A  simple  signal  intensity  budget  shows  that 
there  exists  a  distance  where  the  energy  suffices  to  produce 
a  detection  in  the  passive  receptor  of  the  target,  but  it  is 
not  enough  to  produce  the  same  effect  at  the  receptor  of  the 
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active  device  in  the  searcher .  Geometrical  spreading  and 
absorption  in  the  medium  account  for  the  difference.  The 
active  searcher  is  thus  a  beacon  with  the  choice  of  avoidance 
or  approach  up  to  the  target.  Bearings  from  the  target  to  the 
searchers  change  noticeably  due  to  the  shorter  distances  and 
higher  speed  of  the  actions,  thus  enabling  a  relatively 
accurate  Motion  Analysis  (Wagner,  1989,pp. III-2)  by  the 
target . 

As  a  consequence,  the  target  is  not  only  aware  of  the 
presence  of  the  searcher  but  also  can  have  a  very  clear 
picture  of  the  situation.  It  will  invariably  have  a  mission 
driven  reaction,  either  to  attack  or  evade.  Information 
becomes  crucial  and  the  target  has  the  advantage.  This  feature 
of  the  target  is  the  essence  of  the  problem,  hence  the  name 
RTCAS,  Reactive  Target  Computer  Assisted  Search. 

Time  and  distance  scales  change  as  well  when  transi¬ 
tioning  from  a  passive  to  an  active  search  in  the  case  of 
sonar.  A  passive  search  is  a  resource  employed  in  a  larger 
time  frame  and  over  a  more  extended  area  than  the  active 
search.  Hence  the  name  Area  ASW  for  the  actions  typically 
carried  out  in  the  open  ocean  either  to  "sanitize"  a  corridor 
or  to  survey  a  given  wide  area  for  strategic  purposes. 
Conversely,  active  search  in  ASW  is  a  focal,  more  concentrated 
detection  ranges  are  shorter  hence  the  space  scale 
of  the  search  is  smaller.  The  cost  in  terms  of  allocation  of 
resources  is  higher  in  an  active  search,  hence  its  sustainabi- 


7 


lity  decreases  and  the  time  scales  shrink.  Minutes  count. 

C.  APPROACHING  A  SOLUTION 

The  objective  of  this  thesis  is  to  develop  a  Reactive 
Target  Tactical  Decision  Aid  for  Computer  Assisted  Search. 

TDA' s  are  aimed  to  assimilate  and  conveniently  present 
the  data  (target,  own  assets  and  environment)  to  the  decision 
maker,  i.e.,  transform  the  data  into  information  and  "...ana¬ 
lyze  the  tactical  problem  beyond  what  is  possible  by  humans  in 
a  timely  fashion. . ."  (Wagner,  1989,  pp.I.l) .  That  analysis  may 
either  yield  an  evaluative  picture  from  where  the  decision 
maker  will  construct  a  course  of  action  more  easily  or  produce 
a  recommended  course  of  action.  The  first  step  is  always 
needed  and  crucial  to  the  second,  so  this  thesis  is  con¬ 
strained  to  an  evaluative  TDA  leaving  the  assessment  capabili¬ 
ties  as  a  future  enhancement. 

Given  the  uncertainties  involved  in  the  problem  and 
following  approximately  the  same  approach  as  previous  TDA' s 
did,  RTCAS  uses  a  Monte  Carlo  Simulation  to  represent  a 
Reactive  Target,  a  Fixed  Scan  Model  to  measure  the  detection 
effort  and  a  Bayesian  Filtering  Process  to  fuse  the  simulation 
with  the  real  world  search  activity.  This  section  describes 
those  models. 

1.  The  Information  Fusion  Paradigm 

Figure  1.2  is  a  pictorial  description  of  the  information 
fusion  paradigm  applied  to  a  search  TDA.  Two  main  sources  of 
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information  can  be  distinguished,  reality  and  a  Monte  Carlo 
simulation.  Reality  provides  the  information  from  the  tactical 
situation,  the  datum,  the  time,  and  searcher  activity.  The 
Monte  Carlo  Simulation  provides  a  large  number  of  possible 
locations  and  movement  of  the  target  based  on  the  datum,  each 
one  equally  likely,  at  least  at  the  beginning,  to  be  the  real 
world  target:  the  "prior." 

the  information  consists  of  updating  the  "prior" 
or  probability  map  of  the  location  of  the  target  using  all 
available  information.  There  are  three  main  causes  of  updat¬ 
ing,  namely  the  passage  of  time,  negative  information  and 
positive  information.  In  the  first  case  the  prior  will  move 
with  the  variety  of  courses  and  speeds  provided  by  the  Monte 
Carlo  Simulated  Targets.  A  large  number  of  simulated  tracks 
(500  in  RTCAS)  resemble  many  of  the  possible  positions  and 
movements  of  the  real  target.  If  no  searching  were  conducted 
the  "furthest  on  circle"  would  just  expand  continuously  until 
none  of  the  500  simulated  targets  were  left  in  the  area. 

Updating  for  negative  information  takes  into  account  the 
information  embedded  in  an  unsuccessful  search  over  a  period 
of  time.  During  a  certain  time  interval,  some  simulated 
targets  are  within  the  range  of  detection  of  some  searcher.  If 
the  real  target  is  not  detected,  then  those  simulated  tracks 
become  less  likely  and  those  further  apart  from  the  searcher 
become  more  likely. 

Updating  for  positive  information  takes  place  when  a 
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searcher  gains  a  contact .  Those  simulated  targets  within  a 
certain  distance  from  where  the  contact  is  produced  become 
more  likely  than  others  further  away.  The  probability  map,  the 
posterior  or  "new  prior, "  consists  of  a  partition  of  the 
Local  Area  into  a  number  of  cells.  Each  cell  has  a  probability 
value  which  results  from  adding  up  all  the  individual  likeli¬ 
hoods  or  weights  of  all  simulated  tracks  within  the  cell. 
RTCAS  produces  a  visual  representation  of  the  prior  by  mapping 
cell  probabilities  to  a  color  code. 

This  successive  updating  of  the  individual  weights  of  the 
simulated  targets  (Bayes  Theorem)  and  the  cell  probabilities 
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(Law  of  Total  Probability)  at  a  given  time  interval  is  called 
Bayesian  Filtering  (Stone  et  al , , 1995,pp, 6) . 

2.  Detection  Effort 

In  a  deterministic  detection  model  like  the  "cookie 
cutter"  sensor  (Washburn, 1989) ,  a  target  within  a  certain 
range  is  detected  with  probability  one,  while  those  outside 
are  not  detected  at  all.  The  uncertainties  involved  in  the 
detection  process,  intensified  in  the  case  of  a  shallow  water 
environment  and  a  small  target,  prevent  the  use  of  a  determin¬ 
istic  model  to  evaluate  the  detection  effort.  Instead, the 
Fixed  Scan  Model  (Washburn,  1989, pp. 4)  ,  is  used.  For  each  pair 
simulated  target -searcher,  the  Sonar  Equation  has  to  be  solved 
yielding  a  value  of  Signal  Excess  (Kinsler  et  al.,  1976, pp. 
411)  that  represents  how  much  above  the  Detection  Threshold 
the  signal  return  from  the  target  is.  However,  the  Signal 
Excess  is  considered  a  random  variable  due  to  the  uncertain¬ 
ties  in  the  problem,  so  a  detection  probability  is  calculated 
representing  the  detection  effort.  For  example,  a  simulated 
target  very  close  to  the  searcher  will  have  a  large  value  of 
Signal  Excess  chiefly  because  of  the  small  transmission 
losses,  thus  leading  to  a  high  detection  probability. 

The  same  detection  model  serves  to  quantitatively 
appraise- the  effort  in  an  unsuccessful  search  over  an  area  as 
well  as  the  likelihood  of  a  contact  report.  The  estimate  of 
such  an  effort  is  given  by  the  nondetection  probability  in  the 
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former  and  the  detection  probability  in  the  latter. 

3 .  The  Reactive  Target 

Active  search  turns  the  one  sided  activity  of  passive 
search  into  a  game  of  two  opposed  goals.  Parameterizing  the 
reaction  of  the  target  is  a  complicated  task  because  it 
involves  human  decisions.  Selecting  a  course,  depth,  and  speed 
or  releasing  a  decoy  or  a  weapon  are  the  result  of  a  sophisti¬ 
cated  mental  process  that  takes  place  in  the  target.  The 
information  is  integrated  in  time  to  provide  a  picture  of  the 
threat  the  target  faces  versus  its  own  resources  to  accomplish 
a  goal.  Modeling  such  a  problem  is  a  task  out  of  the  scope  of 
this  thesis.  However,  there  are  some  basic  parameters  such  as 
course  and  speed  that  can  be  considered  the  drivers  of  the 
situation.  RTCAS  uses  those  basic  parameters  to  construct  a 
simple  model  of  the  behavior  of  a  real  target. 

D.  THESIS  STRUCTURE 

Chapter  II  develops  the  core  models  for  the  Reactive 
Target,  the  Bayesian  Filtering  and  the  Detection  Effort. 

Chapter  III  explains  the  implementation  of  the  models. 

Chapter  IV  verifies  the  correctness  and  robustness  in  the 
implementation  and  the  significance  of  a  Reactive  Target  Model 
as  opposed  to  a  Non  Reactive  one. 

Chapter  V  summarizes  future  enhancements  of  the  model  and 
the  implementation. 

Chapter  VI  states  the  conclusions  and  recommendations. 
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II.  THE  ANALYTICAL  MODEL 


A.  THE  TARGET  MODEL 

Three  besic  besks  have  to  be  accomplished,  in  the  target 
model:  generation,  updating  and  reaction. 

1 .  Track  Generation 

The  initial  prior  in  RTCAS  is  provided  by  a  single 
scenario  given  by  a  lost  contact  or  a  flaming  datum,  i.e.,  the 
cue  about  the  position  of  the  target  given  by  a  ship  in  flames 
as  a  consequence  of  a  torpedo  attack. 

Other  TDA  models, (Wagner, 1989,pp. II . 3) ,  allow  the 
generation  of  multiscenario  priors  but  those  are  associated 
with  large  scale  models  that  deal  mainly  with  passive  Area 
Search  where  several  initial  positions  and  estimated  movements 
are  likely  and  should  be  taken  into  account.  In  the  case  of  a 
flaming  datum,  the  expected  movement  is  to  flee  the  area. 

This  initial  prior  is  represented  by  a  bundle  of  500 
simulated  tracks  that  encompass  the  many  possibilities  for 
target  initial  position  and  movement . 

A  Datum  can  be  defined  by  the  following  elements 

•  Datum  identity, 

•  time, 

•  position, 

•  error, 

•  course  and  speed. 
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RTCAS  will  deal  with  a  single  datum,  hence  the  identifi¬ 
cation  of  the  datum  is  irrelevant.  In  a  Local  Area  Search,  in 
general,  all  time  and  space  coordinates  are  considered  rela¬ 
tives  to  that  of  the  datum,  i.e.,  the  geographical  location  of 
the  area  and  the  real  time  are  not  part  of  problem. 

The  time  scale  starts  with  the  searchers'  arrival  in  the 
Local  Area  after  a  delay  time  t^.  (Figure  2.1) 


The  positions  of  the  tracks  representing  the  datum  are 
referred  to  the  center  of  a  grid  that  covers  the  Local  Area 
where  the  search  is  conducted.  The  coordinates  of  any  track  i 
are  [xTi,yTi]  (t)  at  time  t. 

The  error  is  taken  into  account  when  the  tracks  are 
generated  and  is  assumed  to  be  bivariate  normal, 

XTi(-T^)  ^  N  (/ix  /  (2.1) 

•Yxi  ( —  T^j)  **  N  ( /^Y  ,  C^Y^  / 

where  the  mean  =  fiy  are  assumed  to  be  equal  to  the  center 
of  the  Local  Area  and  the  standard  deviations  input 

by  the  user.  The  notation  X  ~  N()  means  that  X  is  a  random 
variable  normally  distributed. 
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The  magnitude  of  each  target's  velocity  vector  is  given 
by  a  random  variable  uniformly  distributed  between  a  lower  and 
an  upper  bound  given  by  the  user  and  does  not  change  unless 
there  is  a  reaction  because  of  the  presence  of  the  searcher. 
(See  Section  3) . 

The  course  is  also  a  random  variable  uniformly  distribut¬ 
ed  between  user  input  limits .  The  purpose  of  this  model  is  to 
resemble  a  situation  where  a  preferred  course  of  evasion 
exist s,e,g.,  windward,  towards  its  mission  essential  target, 
navigational  hazard  avoidance,  etc.  The  situation  where  no 
preferences  may  be  estimated  about  the  target's  intentions 
uses  the  same  model  including  all  possible  courses  between  0 

and  360  degrees,  thus  leading  to  an  omnidirectional  fleeing 
datum . 

The  delay  time  t  is  taken  into  account  by  updating  the 
position  of  all  tracks  using  each  one's  velocity  vector  to 
obtain  the  initial  target  distribution. 

[XTi.y^j]  (0)  =  (-T^)  +TVf  (-T^)  (2.2) 

A  typical  target  distribution  at  time  zero  can  be  seen  in 
Figure  2.2  representing  the  scenario  the  searchers  find  when 
they  arrive  at  the  Local  Area.  To  the  left  of  the  figure,  the 
target  is  assumed  to  flee  omnidirectionally.  To  the  right  is 

the  case  when  the  target's  course  is  between  10  and  100 
degrees. 


2 .  Track  Updating 

At  every  time  increment  t  the  position  of  all  the  targets 
as  a  function  of  time  is  updated  using  each  track's  velocity 
vector, 

( t+T)  =  +Tyj  ( t)  (2.3) 

The  speed  of  the  target  remains  unchanged  unless  there  is 
a  reaction. 

3.  Track  Reaction 

The  reactive  updating  of  the  tracks  consists  of  changes 
in  the  properties  of  each  track  as  a  consequence  of  the 
presence  of  the  searchers.  All  reactions  of  a  real  target  may 
be  reflected  in  a  variation  of  speed,  course  or  depth.  RTCAS 
updates  the  course  only.  The  variations  in  depth  are  more 
likely  to  occur  in  deep  water  environments.  Speed  alterations, 
in  general,  are  made  within  ranges  that  do  not  have  a  great 
impact  in  the  overall  geometry  and  kinetics  of  the  problem. 
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The  target  has  two  goals,  evade  the  datum  and  evade  the 
searchers .  In  RTCAS  the  simulated  targets  evade  the  datum 
unless  they  have  to  evade  a  searcher  because  of  its  proximity. 
The  assumption  is  made  that  the  target  knows  the  searcher's 
location.  This  is  not  too  strong  an  assumption  in  active 
search.  In  general,  the  most  common  information  the  target 
knows  about  the  searchers  is  the  bearing.  In  the  close 
distances  normally  implied  in  Local  Area  actions,  the  bearing 
varies  rapidly  thus  allowing  the  target  to  perform  an  adequate 
passive  motion  analysis  (Wagner,  1989, pp.  Ill-i)  .  it  is  assumed 
that  the  target  reacts  when  a  searcher  at  a  distance  smaller 
than  some  threshold. 

The  simplest  evasion  maneuver' is  to  adopt  a  course  that 
generates  a  direction  of  relative  motion  (DRM)  opposed  to  that 
of  collision.  Figure  2.3  ,  left  depicts  such  a  situation.  Most 
of  the  time,  this  will  be  impossible  for  the  target  because 
its  speed  will  be  smaller  than  that  of  the  searcher  (Figure 
.  A  discussion  of  all  possible  courses  of  action  a- 
vailable  to  the  target  at  this  point  is  out  of  the  scope  of 
this  thesis.  However,  there  are  some  constraints  on  what  the 
target  should  do.  It  should  not  continue  at  the  same  course  if 
the  distance  is  decreasing,  nor  should  it  turn  and  end  up  with 
a  relative  bearing  of  +90  degrees  to  the  searcher,  since  this 
would  increase  dramatically  the  target  strength  in  the  sonar 
equation.  A  180  degrees  relative  bearing  should  also  be 
avoided  since  this  may  imply  a  reduction  in  the  quality  of 
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Figure  2.3.  The  evasion  alternatives  as  a  function  of  the  velocities 
of  the  searcher  Vg  and  the  target  Vi. 

information  obtained  by  the  target  because  of  the  blind  sector 
in  the  stern.  For  simplicity,  it  is  assumed  that  target  i 
turns  a  fixed  angle  60  e (90®, 180°)  off  the  bearing  to  the 
searcher.  This  new  course  for  target  i  is, 

dB 

Qi=Bis-Adsign{-^)  (2.4) 

where  Dbj_^/dt  is  the  bearing  rate.  If  the  bearing  increases 
(decreases) ,  the  target  turns  to  a  course  60  degrees  measured 
counterclockwisely  (clockwisely)  from  It  is  the  tangential 
component  of  the  searcher's  speed  relative  to  the  target's 
velocity  vector  that  determines  the  sign  of  the  last  term  in 
(2.4)  and  it  is  calculated  in  RTCAS  with  the  following 
formula, 

dB. 

sign{-^)=sign{s±n{Q^-B^;))  (2.5) 

where  0.  is  the  searcher's  course . (Figure  2.4) 


2.4.  The  evasive  course  parameters 

When  facing  several  searchers  within  the  reaction  range, 
the  target  decides  to  evade  the  closest. 

B.  THE  BAYESIAN  FILTERING  MODEL 
1.  The  Probability  Map 

The  final  purpose  of  the  model  is  to  provide  a  "probabil¬ 
ity  map"  representing  the  location  of  the  target  as  an  aid  in 
conducting  a  Local  Area  Antisubmarine  Active  Search. 

The  search  starts  with  a  "flaming  datum"  (Cheong,  1988, 
pp.l),  hence  a  relatively  small  time  and  space  frame  contains 
the  action,  e.g.,  1  to  4  hours  and  an  area  with  a  radius  from 


10  to  20  nautical  miles. 


The  probability  map  is  a  partition  of  the  Local  Area  into 
a  grid  of  KxL  cells  each  of  which  has  a  probability  value 
associated  with  it.  We  define  the  following  events: 

I  :{  "The  target  is  within  the  Local  Area") 

Gki  :{  "The  target  is  in  the  grid  cell  k,l"} 

Then  the  probability  map  is  the  bivariate  probability 
mass  function  of  P(G]^i\I).  The  idea  of  a  flaming  datum  as  the 
starting  point  of  the  search  makes  the  condition  on  the  fact 
that  the  target  in  within  the  Local  Area  rather  irrelevant  at 
the  start.  However,  as  times  goes  by,  the  influence  of  I 
increases  as  the  target  may  have  already  left  the  box  where 
the  search  is  conducted. 

In  order  to  build  the  probability  map  it  is  necessary  to 
find  a  way  to  assign  a  value  to  P(G]^i\l) . 

This  Bayesian  analysis  begins  with  a  single  scenario 
prior  distribution  of  the  "state  of  nature"  (Washburn, 
1995, pp.l)  represented  hy  N  =  500  possible  tracks.  Intuitive¬ 
ly,  when  conducting  a  Local  Area  search,  one  make  some 
assumptions  and  evaluate  a  series  of  "what  if"  questions, 
reacting  to  the  most  evident  and  easy  to  maintain  in  one's 
mind  as  a  chain  of  causes  and  effects.  In  this  case  N  simulat¬ 
ed  tracks  try  to  keep  that  chain  alive  for  as  long  as  possi¬ 
ble. 

The  prior  is  generated  by  the  Target  Model  as  described 
in  the  previous  section.  The  output  of  this  model  is,  at  this 
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stage,  the  stochastic  location  and  movement  of  each  one  of  the 
N  simulated  tracks  according  to  the  assumed  prior.  We  define 
the  event, 

Ti:  {"The  real  target  follows  track  i",  i=l.  .n) 
then,  if  all  tracks  are  equally  likely  to  be  the  real  target, 
assumption  embedded  in  the  idea  of  a  single  scenario,  then  the 
initial  "weight"  of  each  track  is 

P[rj{0)=l  (2.6) 

and  the  initial  probability  map  can  be  depicted  as 

E  P^T,]{t)  (2.7) 

where  =  (i  :  [Xri,yTi]  (t)  z  t^o} . 

Ideally,  if  there  were  no  searchers  and  (2.7)  was 
calculated  as  a  function  of  time,  it  would  convey  almost  the 
same  information  as  the  expanding  furthest-on-circle  (Cheong, 
1988,  pp.l)  .  However,  the  purpose  of  RTCAS  is  to  integrate  the 
information  from  the  real  world  provided  by  the  assets 
conducting  the  search  together  with  the  information  provided 
by  the  prior  (Stone  et  al . , 1995,pp. 1) .  The  search  done  up  to 
a  certain  time  t  may  be  unsuccessful  or  result  in  one  or  more 
searchers  obtaining  a  contact,  which  requires  updating  the 
prior  for  negative  or  positive  information  respectively  (Wag¬ 
ner,  1989, pp.  II  .  5)  . 
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2.  Updating  for  Negative  Information. 

The  valuable  information  of  an  unsuccessful  search  during 
the  interval  t  is  that  the  real  target  is  less  likely  to  be  in 
the  cells  that  had  been  searched  as  a  result  of  the  detection 
effort  imposed  to  the  tracks  located  in  or  near  those  cells. 
The  value  of  each  track  weight  will  be  decreased  in  accordance 
with  that  effort.  The  track  weights  will  change  from  the  prior 
values  at  time  t  to  the  posterior  values  at  time  t+r,  hence 
changing  the  probability  map  of  the  location  of  the  target. 

Define  the  event, 

Dst  "The  target  is  detected  by  the  searcher  s 

during  time  interval  t"  } 

then 

cdp^i^  it)  =P[D^^\T^]  it)  (2.8) 

is  the  cumulative  probability  that  the  target  was  detected  by 
the  searcher  s  during  the  time  interval  t  starting  at  time  t 
conditioned  on  that  the  real  target  was  located  at  the 
position  of  simulated  target  i  during  that  interval.  This 
value  is  provided  by  the  Detection  Model  to  be  developed  in 
section  C. 

If  D'sr  is  the  complement  of  ,  then 

P[Di, !  rj  ( t)  =l-cdp^^,  (t)  (2.9) 

is  the  nondetection  probability  of  searcher  s  over  target  i 
during  interval  t  starting  at  time  t.  Consequently,  the 
probability  that  the  target  is  following  track  i  and  is  not 
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detected  by  searcher  s  during  the  interval  t  is. 


p  [  (t)  =p  [Ds^  I  rj  ( t)  p  [  rj  ( t) 

Using  Bayes  Theorem,  the  Negative  Information 
prior  is 

i€N 


(2.10) 

Updating  of  the 

(2.11) 


Note  the  convention  that  PETJCt+r)  is  conditioned  on  all 
information  available  up  to  time  t+r. 

The  denominator  in  (2.11)  represents  the  unconditional 
probability  (t)  obtained  using  the  Total  Probability 
Theorem.  Furtheormore  it  is  the  normalizing  factor  that  will 
allow  the  new  prior  to  be  a  "proper"  probability  mass  function 
i .  e  -  P[TJ  (t+T)  =  1. 

An  iterative  application  of  this  process  permits  the 
calculation  of  the  new  prior  for  any  time.  The  update  2.11 
must  be  applied  once  per  time  interval  for  each  searcher  that 
does  not  detect  the  target. 

3.  Updating  for  Positive  Information 

The  probability  map  should  be  updated  for  Positive 
Information  when  a  searcher  reports  a  contact.  It  may  be 
contended  that  once  a  contact  is  obtained,  the  relative  value 
of  the  information  provided  by  the  simulation  of  tracks  is 
somehow  negligible  with  respect  to  the  real  world  information 
implied  in  the  contact  report.  However,  in  a  shallow  water 
acoustic  environment  where  many  elements  of  the  "clutter" 
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space, e.g.,  a  whale,  a  rock,  an  eddies,  may  be  considered  a 
contact,  the  fusion  of  information  of  both  sources  remains  in 
place  as  a  means  to  smoothing  the  reactions  and  providing 
continuity  whenever  the  contact  is  confirmed  not  to  be  the 
desired  target . 

Define  the  following  events 

CgcT  :  {  "The  searcher  s  reports  a  contact  at  the 
position  c  during  the  interval  t",  s=1..S, 

c=l..C} 

where  S  is  the  total  number  of  searchers  and  C  is  the  total 
number  of  reported  contacts  among  all  searchers,  and, 

:{  "The  real  target  is  located  at  the  position  of 

the  contact  report  c  during  the  time  interval 
t"} 

then, 

(2.12) 

where,  as  before,  Dg^  is  the  event  that  searcher  s  detects  the 
target.  Thus, 

P[Cg„!rj  (t)=P[Dg,lrJ  it)P[EjD^,,Ti]  (t)  (2.13) 

The  first  factor  in  2.13  is  the  detection  effort  alloted 
by  searcher  s  at  the  position  of  simulated  track  i  during 
period  t, 

P  [Dg,  !  EJ  ( t)  =cdp^i^  it)  (2.14) 
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uted  the  second  factor  is 


where 


P[E^\D^,T^\  =ke 


(2.15) 


ricT  is  the  distance  between  the  position  of  the  contact 


the  sensor  and  r^i  is  the  distance  between  searcher  s  and 
tracki,  then  RTCAS  takes  a^i  to  be  Q  r^j2.  Thus  2.15  becomes. 


.Ti\  =e 


-2^jei 

Qkli 


(2.16) 


where  the  constant  k.  was  dropped  because  it  cancels  when  Bayes 
Theorem  is  applied.  Substituting  2.14  and  2.16  into  2.13, 


P  [  C^„ !  r^]  it)  =cdpg^^  it)  e  ‘ ^ 


(2.17) 


represents  the  probability  that  searcher  s  reports  a  contact 
at  position  c  given  that  the  target  follows  track  i  during  the 
interval  t.  Equation  2.17  embodies  the  mathematical  core  of 
fusing  information  between  the  real  world  contact  report  and 
the  simulation. 

The  probability  that  the  real  target  follows  track  i 
given  that  the  searcher  s  reported  a  contact  at  the  position 
of  c  during  period  t  applying  Bayes  Theorem  is  (Washburn, 
1995, pp. 6) , 


P[T^\C,„]{t)=- 


P[C^„I  rjp[rj  (t) 

N 

J2pIC,„\T,]P[T,]  (t) 


2=1 


(2.18) 


which  using  2.17  becomes 


^le 


N  -TT 

(t) 

i=l 


(2.19) 
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where  the  +  superscript  implies  that  the  contact  is  actually 
a  true  contact  on  the  target.  In  practice,  updates  for 
positive  information  must  allow  for  the  possibility  of  false 
alarms,  so  the  formula  actually  used  by  RTCAS  is, 

P{T^]  (t+T)=P[rj]*(t+T)C)+P[2t]  (t)  {1-Q)  (2.20) 

where  Q  is  the  "credibility"  (Washburn,  1995, pp. 2)  of  the 
contact  report.  When  Q=l,  2.20  is  the  same  as  2.19,  but  when 
0=0,  2.20  is  just  the  prior  P[Tj(t). 

It  might  be  contended  that  the  posterior  alternate  to 
P[Ti]*  (t+T)  should  include  the  negative  information  implied  in 
the  unsuccesf ul  search  that  a  failed  contact  report  conveys . 
However,  it  is  very  plausible  and  commonly  seen  in  the 
practice  that  a  searcher  investigating  a  contact  will  focus 
onto  it,  not  only  because  of  a  natural  tendency  in  the 
distribution  of  the  attention  of  the  sonarman,  but  also 
because  investigating  a  contact  requires  more  often  than  not 
using  modes  or  techniques  in  the  sensors  that  prevent  the 
operator  from  an  adequate  scanning  of  the  rest  of  the  picture 
even  without  intention  of  doing  so.  Examples  of  those  modes  or 
techniques  are  doppler  and  image  analysis,  noise  reduction  by 
windowing  the  signal  processing  in  the  target,  angular 
scanning  reduction  spotting  to  concentrate  acoustic  intensity 
on  the  target,  etc.  As  a  consequence  and  for  computational 
convenience,  RTCAS  considers  an  absolute  focusing  on  the 
contact,  i.e.,the  individual  track  probability  alternative  to 
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the  one  updated  for  positive  information  is  the  value  at  the 
beginning  of  the  scan  period  as  in  equation  2.20. 

4.  The  undetected  probability  mass 

Although  each  probability  map  is  the  best  estimate  of  the 
location  of  the  target,  it  says  nothing  about  how  much  of  the 
original  probability  mass  remains  undetected.  A  measure  of  the 
effect  of  the  search  is  necessary  since  the  searcher  needs  to 
know  how  much  probability  mass  is  still  undetected  in  order  to 
be  able  to  decide  when  to  terminate  the  search. 

The  denominator  in  (2.11)  is  the  probability  that  the 
real  target  is  following  one  of  the  500  simulated  tracks  and 
was  not  detected,  hence  is  the  undetected  probability  mass 
during  the  interval  t  at  time  t+r: 

5,(t)=5;p[r,nD^]  (2.21) 

i€N 

for  instance,  at  time  t=0,  before  any  search  is  conducted 
S(t)=l.  At  each  updating  step  the  searching  will  reduce  the 
original  probability  mass  such  that  at  any  time  t,  the 
Cumulative  Undetected  Probability  Mass  is 

U{t)=U{t-T)  Sjt)  (2.22) 

This  recursive  definition  allows  the  calculation  of  its  value 
at  every  time  step  in  order  to  produce  the  desired  assessment . 

C.  THE  DETECTION  MODEL 

Several  possibilities  were  taken  into  account  in  select - 


28 


ing  the  detection  model.  The  idea  is  to  model  active  search. 
The  Fixed  Scan  Model  seems  to  be  naturally  applicable  since 
(Washburn, 1994,  pp.  3): 

•  Each  pulse  can  be  thought  of  as  an  independent 
look, 

•  Detection  in  the  absence  of  a  pulse  is  impossible, 

•  There  exists  a  pulse  rate  which  can  be  associated 
with  the  pulse  repetition  frequency  of  a  sonar. 

The  first  statement  is  seemingly  too  strong  an  assumption 
but  in  RTCAS  only  the  conditional  independence  of  the  pulses 
given  the  target  track,  is  required.  The  validity  of  the 
assumption  depends  on  the  rapidity  of  fluctuations  in  the 
sonar,  equation  terms  relative  to  the  time  between  active 
pulses. 

The  pulse  repetition  frequency  (prf)  is  normally  associ¬ 
ated  with  the  maximum  range  at  which  it  is  desired  to  attempt 
a  detection.  Given  a  typical  value  of  maximum  scale  range  in 
active  sonar  of  about  18000m,  the  round  trip  of  a  pulse  takes 
24  seconds  at  a  sound  speed  of  1500m/s,  thus  yielding  a  prf  of 
approximately  2  per  minute.  In  conventional  active  sonar  the 
frequency  is  normally  greater  than  4  kilohertz.  The  strongest 
time  variability  in  the  propagation  is  given  by  the  time 
change  in  the  refraction  patterns  which  in  turn  depend  on  the 
depth-velocity  profile.  In  shallow  waters,  the  fastest  change 
in  the  upper  layer  of  the  velocity  profile  can  be  regarded  as 
semi-diurnal  as  a  consequence  of  the  heat  interchanges 


29 


provoked  by  the  daily  cycles  of  the  sun.  Nonetheless,  the 
effect  of  cloud  cover,  rain,  wind  or  other  many  factors  may 
produce  noticeable  changes  in  the  propagation  conditions. 
Furthermore,  there  are  some  changes  in  the  geometry  of  the 
detection  problem  that  may  cause  variations  in  the  time  frame 
of  interest.  For  instance,  given  a  frequency  of  4.5  kilohertz, 
the  wavelength  will  be  0.3  meters.  In  shallow  waters,  it  is 
very  likely  that  the  receiver  will  sense  a  multipath  scatter¬ 
ing  from  the  target.  If  the  difference  between  two  paths  is  a 
multiple  of  0.15  meters,  there  will  be  destructive  interfer¬ 
ence  thus  leading  to  an  eventual  lack  of  detection.  Although 
a  deterministic  geometric  model  may  account  for  each  situation 
regarding  the  relative  positions  of  target  and  searcher  with 
respect  to  the  bottom,  the  surface  and  the  layer,  the  outcome 
of  such  a  model  may  be  regarded  as  randomly  distributed  with 
respect  to  space  and  time.  Other  phenomena  may  concur  in  order 
to  increase  the  randomness  in  the  signal  at  the  receiver 
between  different  pulses,  for  instance,  the  presence  of  the 
bubbles  due  to  the  wake  of  other  ships  or  the  target,  changes 
in  the  angle  of  incidence  of  the  pulse  at  the  target,  the 
destructive  interference  with  the  signals  from  other  plat¬ 
forms,  the  macroscopic  biological  activity,  etc.  All  of  these 
sources  of  randomness  are  typical  of  the  shallow  water  active 
signal  propagation  problem. 

The  assumption  in  RTCAS  is  that  each  pulse  is  condition¬ 
ally  independent  of  the  others.  Of  course  the  computational 


30 


in  considering  nonindependence  had  to  be  taken  into 
account  and  favored  the  chosen  option. 

1.  The  Fixed  Scan  Model 

This  thesis  applies  the  Fix  Scan  Model  as  a  paradigm  of 
the  use  of  hull  mounted  and  dipping  active  sonars  and  sono- 
bouys  in  searching  for  a  conventional  submarine.  However,  it 
is  applicable  to  other  types  of  sensors  like  radar,  lidar, 
photo  electronic  devices,  etc.  as  well  as  other  types  of 
targets  like  surface  ships,  life  rafts,  land  vehicles,  etc. 

The  sensor  is  assumed  to  emit  pulses  at  a  repetition 
frequency,  prf,  given  by 

(2.23) 

where  C  is  the  speed  of  the  sound  in  the  water  and  M  is  the 
maximum  unambiguous  range,  assumed  always  to  be  equal  to  the 
maximum  range  scale  selected  in  the  device. 

The  probability  of  detection  will  be  calculated  only  for 
those  targets  within  the  maximum  range  scale  of  any  searcher 
at  the  end  of  every  updating  time  inteirval  t.  The  value  of  t 
must  satisfy  a  compromise  between  the  need  of  an  adequate 
"sampling"  of  the  search  process  and  the  computational 
difficulty  of  too  many  calculations  in  a  given  time  frame. 

As  an  example,  for  a  maximum  range  scale  of  18000m,  the 
interval  between  pulses  will  be  24  seconds,  which  may  be 
approximated  as  30  seconds.  Hence,  if  the  updating  interval  is 
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T  =  5  minutes,  10  independent  scans  are  going  to  be  made  with 
the  same  sonar  operating  characteristics,  i.e.,  same  pulse 
type,  prf,  frequency,  etc.  Selecting  5  minutes  as  the  updating 
interval  is  convenient  from  the  point  of  view  of  the  computa¬ 
tional  effort  as  well  as  tactically  reasonable;  e.g.,  the  rate 
of  turn  of  a  conventional  submarine  at  6  knots  is  such  that 
it  would  take  about  4  minutes  to  complete  a  180°  turn. 

During  one  updating  interval  t,  both  the  searcher  and  the 
target  are  assumed  to  remain  stationary  at  a  point  midway 
between  their  respective  locations  at  the  beginning  and  at  the 
end  of  T  (Figure  2.6),  for  purposes  of  detection  probability 
calculations . 

Assuming  the  above  mentioned  independence  between  each 
scan,  (Washburn, 1995, pp. 3) ,  the  probability  that  at  least  one 
of  the  pulses  k  from  searcher  s  detects  the  target  given  that 
it  follows  track  i,  at  time  t,  is, 

cdp„^  ( t)  =l-ni=i  ^^-Pksi  ( t) )  (2.24) 

where  T,  the  total  number  of  pulses  during  the  updating 
interval  t,  may  be  calculated  as, 

T=x*pxf  (2.25) 

The  detection  probability  of  each  pulse  k  is  given  by, 

SF!  ( "t) 

Pjcsi(t)  - .  )  (2.26) 

where  4>  is  the  Standard  Normal  Cumulative  Density  Function  and 
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SE^i(t),  calculated  by  the  Propagation  Model  to  be  described 
in  the  next  section,  is  the  Signal  Excess  from  target  i  at  the 
receiver  in  searcher  s  at  time  t,  and  is  the  standard 
deviation  in  the  calculation  of  due  to 

•  the  randomness  in  the  detection  process,  and 

•  the  cumulative  error  or  uncertainties  in  the  calcu¬ 
lation  of  each  one  of  the  terms  in  the  sonar  equa¬ 
tion. 

The  value  of  can  be  estimated  either  by  experience  or 
by  using  a  property  of  the  variance  of  the  sum  of  random 
variables  (Devore, 1991, pp. 102) , 


33 


(2.27) 


s 

E 

N  ic=i 

where  is  the  error  in  the  calculation  or  estimation  of  each 
k  of  the  S  terms  in  the  sonar  equation.  Wagner  (1989, pp.  II- 
28)  suggests  values  of  6-9  dB  for  a„',  RTCAS  uses  15  dB  since 
it  is  intended  for  use  with  relatively  poor  sonar  forecasts. 

2.  The  Propagation  Model 

This  model  corresponds  to  the  propagation  of  active  sonar 
signals  in  a  shallow  water  environment.  The  values  of  all 
terms  in  the  sonar  equation  will  be  either  assumed  to  be 
constants  or  approximated  by  simple  models.  It  is  the  user 
that  will  finally  input  the  values  corresponding  to  more 
accurate  calculations . 

For  each  pair  searcher  s  and  target  i  at  time  t  the 
signal  excess  is  be  calculated  using  the  active  sonar  equa¬ 
tion,  (Kinsler  et  al.  1982)  limited  either  by  noise, 

SEN^i=SLg-DNL^{t)  -DT^+DI^-2TL^^  ( t)  +TS^^  ( t)  (2.28) 

or  by  reverberation, 

SER^i  ( t)  ( t)  -DT^-SS^+TLG^j^  ( t)  -lOlog^o  ( )  (2.29) 

All  levels  in  the  formulae  above  are  in  decibels  with 
respect  to  one  microPascal.  The  Source  Level  SL,  the  Detection 
Threshold  DT  and  the  Directivity  Index  DI  are  considered 
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constant  for  each  searcher  s.  The  rest  of  the  terms  are 
variables  calculated  by  RTCAS  and  are  developed  in  the 
following  sections.  DNL^  is  the  Detected  Noise  Level  for 
searcher  s,  TL^j_(t)  is  the  Transmission  Loss  for  each  pair 
(s,i)  at  time  t,  TS^i(t)  is  the  Target  Strength  of  target  i 
with  respect  to  searcher  s,  SS  is  the  Scattering  Strength  for 
the  sonar  in  searcher  s,  TLGsi(t)  is  the  transmission  loss  due 
to  geometrical  spreading  and  the  last  term  in  (2.10)  is  the 
area  of  the  reverberating  surface  within  the  beam  of  the 
transmitted  pulse  corresponding  to  the  speed  of  the  sound  c, 
the  pulse  length  t  and  the  horizontal  angle  of  the  sonar  beam 
(Kinsler  et  al,,1892).  The  smaller  of  (2.28)  and  (2.29)  is 
to  be  selected  and  used  in  calculating  (2.26) . 

In  active  search  in  shallow  waters,  the  sonar  equation  is, 
in  general,  more  likely  to  be  limited  by  reverberation  rather 
than  by  noise.  However,  in  some  cases  the  effect  of  the  noise 
will  prevail  at  distances  within  the  maximum  range  scale  since 
the  echo  level  and  the  reverberation  level  decrease  with  range 
but  the  noise  level  remains  the  same.  A  typical  behavior  can 
be  seen  in  Figure  2.7,  where  up  to  7500  meters  the  noise 
restricted  Signal  Excess  is  higher  than  the  reverberation 
restricted  one . 

The  following  sections  discuss  the  terms  in  the  sonar 
equations  (2.9)  and  (2.10)  that  are  not  considered  a  constant. 
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Figure  2.7.  Noise  aind  Reverberation  Limited  detection. 


3 .  Target  Strength 

Target  strength  is  defined  as  the  echo  returned  by  the 
target  and  hence  is  a  factor  when  dealing  with  active  sonar. 
The  corresponding  term  in  the  sonar  equation  is: 

r5^,=10logio(i),r=l 

■‘■rO 

where  Ij,  is  the  intensity  of  return  at  1  yd  from  the  target 
and  Ijj  is  the  incident  intensity. 

The  intensity  ratio  Ir/Ii  is  a  function  of  distance 
between  the  target  and  the  source  (Urick,  1975,pp.284) ,  the 
shape  of  the  target  (Urick, 1975, pp. 274) ,  the  pulse  length 
(Urick, 1975, pp. 285) ,  the  frequency  of  the  incident  signal 
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(Urick,  1975, pp. 283)  ,  and  the  geometry  of  the  reflection  (U- 
1975 , pp .  281)  .  The  influence  of  all  those  factors  except 
the  last  can  be  considered  constant. 

In  the  case  of  a  monostatic  sonar  the  direction  of 
incidence  is  the  same  as  the  direction  in  which  the  signal  is 
received  at  the  transducer  except  for  small  variations  due  to 
the  relative  motion  target -searcher  during  the  time  the  pulse 
travels.  The  target  strength  is  a  function  of  the  angle 
between  the  heading  of  the  submarine  and  the  angle  of  inci¬ 
dence  of  the  signal,  which  is  called  aspect  angle.  If  lateral 
symmetry  is  assumed,  then  the  aspect  angle  can  be  considered 
from  0  to  ±  IT,  with  0  being  at  the  bow.  The  great  variation  of 
TS  as  a  function  of  the  aspect  angle  reguires  the  aspect  angle 
to  be  considered.  The  Target  Strength  will,  in  general,  vary 
from  its  maximum  value  at  broadside  to  a  minimum  at  the  bow 
and  stern. 

Let  r  be  the  distance  target -searcher  and  v  the  target 
speed,  then, 

cos(0)=f'i>^  (2.31) 

^^^^*-^bes  the  dot  product  between  the  corresponding  unit 
vectors  and  0  is  the  smallest  angle  between  them,  in  this  case 
the  aspect  angle  (Figure  2.8). 

The  typical  values  of  TS  vary  from  about  25  dB  at  beam 
aspect  to  about  15  dB  at  bow  and  stern  aspects  (Urick,  1975 
p.282)  .  The  following  equation  is  an  approximation  adopted  as 
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a  simplified  model  of  the  variations  between  those  extremes. 


TS(e)  =TS^+ATS*{l-cos^{d))  (2.32) 

where  TS^  is  the  value  of  Target  Strength  at  bow  and  stern, 
ATS  is  the  typical  increment  to  reach  the  value  at  beam  aspect 
and  (1-cos^  (d) )  models  the  variation  as  a  function  of  the 
aspect  angle  0.  Substituting  (2.12)  into  (2.13)  yields, 

( t)  =rS^+A rs ( 1-  (4^  ( t)  ( t) )  2)  (2.33) 

where  rgi(t)  is  the  instantaneous  distance  between  searcher  s 
and  target  i,  and  v^it)  is  the  speed  of  target  i  at  time  t. 
RTCAS  uses  a  square  of  the  cosine  as  a  model  because  of 
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computational  convenience.  Note  in  2.33  that  the  trigonometric 
function  in  2.32  is  replaced  by  a  dot  product,  which  can  be 
calculated  using  the  components  of  the  corresponding  unit  vec¬ 
tors  instead  of  the  trigonometric  function. 

4.  Transmission  Loss 

Transmission  loss  involves  the  combined  effect  of  geomet¬ 
rical  spreading  and  absorption. 

For  simplicity,  assuming  that  the  distance  (t)  is  much 
than  the  depth  H,  geometrical  spreading  is  calculated 
with  the  following  isospeed  shallow— water  channel  perfectly 
rigid  bottom  expression  (Kinsler, 1984, pp. 429) , 

rLG^^(t)  =10log[r^^(t) +-^)]  (2.34) 

71 

where  the  last  term  is  the  correction  for  the  transition 
range . 

The  total  Transmission  Loss  is 

{ t)  =TLG^^  ( t)  it)  (2.35) 


where  a,  the  absorption  coefficient,  is  calculated  with 


8xl0'5 


0.7 


+1 


^  0.04 

IPOO+l 


+4xl0-’f2 


(2.36) 


from  Kinsler  et  al.  ( 1976, pp. 3  98)  where  f  is  the  frequency  in 
Hertz  at  which  each  sonar  operates. 
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5. 


Noise  Level 


The  Noise  Level  accounts  for  the  background  ambient  noise 
(an)  and  the  platform  self  noise  (sn) , 


NL=20log^Q 


an+sn' 

l|iPa 


(2.37) 


Instead  of  calculating  the  intensities  an  and  sn  their 
corresponding  levels  ANL  and  SNL  in  decibels  are  interpolated 
from  observational  data. 

The  Ambient  Noise  Level  is  calculated  by  means  of  a 
third  order  polynomial  fit  on  windspeed  and  a  linear  fit  on 
the  logarithm  of  frequency  of  observational  data  obtained  from 
curves  representing  the  noise  spectra  for  some  coastal 
locations  in  Figure  7.7.  in  Urick  (1975, pp. 191) , 


AiVL=46.12+2.22v/-0.041;/2_i7_Oipogio(f)-3]  (2.38) 

where  w  is  the  wind  speed  in  m/s  and  f  is  the  frequency  in 
hertz.  This  approximation  is  centered  around  f  =  1000  Hz. 
Figures  2.9  and  2.10  compare  Equation  2.38  to  the  data  from 
the  above  mentioned  reference . 

The  Self  Noise  Level  is  calculated  as  a  function  of 
frequency  and  the  searcher's  speed  using  figure  11.11  in 
Urick(1975,  pp.340),  which  shows  equivalent  isotropic  self¬ 
noise  levels  at  25  khz  on  a  number  of  World  War  II  American 
and  British  destroyers.  A  polynomial  fit  of  this  data, 
including  an  extrapolation  term  for  other  frequencies  assuming 
a  slope  of  -6dB  per  octave,  yields. 
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Ambient  Noise 

Windspeed  as  a  parameter 


1  “6  knots 

Hi 

7-1 0  knots 


11-19  knots 


Fx9^^6  2  •  9  •  The  Ambient  Noise  as  a  function  of  frequency  and  the 
windspeed  as  a  parameter  between  1  and  16  knots. 


SNL=23+2.02v^{t)  +20logi-^)  (2.39) 

where  is  the  searcher's  speed  in  knots  and  f  is  the 
frequency  in  kilohertz.  Figure  2.11  compares  this  equation 
with  the  corresponding  data.  SNL  and  ANL  are  power-summed  to 
obtain  NL. 

6 .  Scattering  Strength 

According  to  Urick  ( 1975, p. 253)  ,  when  downward  refraction 
occurs,  i.e.,  negative  sound  speed  profile,  bottom  back 
scattering  is  dominant  and  the  reverberation  level  can  be 
calculated  using  the  bottom  scattering  strength.  Conversely, 
when  upward  refraction  occurs,  i.e.,  positive  sound  speed 
profile  as  in  isothennal  water,  the  surface  backscattering 
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Ambient  Noise 

Wlndspeed  as  a  parameter 


(Thousands) 


17-21  knots 

22-27  knots 

A 

28-36  knots 


Fx^uitb  2  •10*  The  Ambient  Noise  as  a  function  of  frequency 
windspeed  as  a  parameter  between  17  and  35  knots. 


and 


becomes  most  important.  This  model  follows  that  rule. 

Following  Urick  (1975, pp. 241) , the  Surface  Scattering 
strength  is  calculated  with, 

S'S'L=10logio  (fxhxsinCe)  )°'^®-45.3  (2.40) 

where  f  is  the  frequency  in  hertz,  6  is  the  grazing  angle  in 
degrees  and  h  =  0.0026  where  W  is  the  windspeed 
in  meters  per  second.  In  a  shallow  water  environment , the 
backscattering  of  the  sonar  signal  will  occur  at  all  possible 
grazing  angles  due  to  the  multiple  interactions  between  the 
surface  and  the  bottom.  Due  to  the  analytical  and  computation¬ 
al  difficulty  in  accounting  for  those  angles,  0  was  chosen  to 
be  35°,  an  average  between  the  typical  limit  angles  in  Figure 
8.21  in  Urick(1975,  pp.  240). 
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Self  Noise 

1kHz 
5kHz 
10kHz 
25kHz 
25kHz(Obs,) 


Figure  2.11.  self  Noise  data  for  25kHz,  its  polynomial  fit  and  the 
extrapolations  for  other  frequencies  using  a  slope  of  -6dB  per  octave. 

In  the  calculation  of  bottom  strength,  Urick  (1970, pp. 
396)  found  that  the  spectrum  slope  in  shallow  waters  was 
similar  to  measurements  made  in  deep  waters,  therefore  figure 
8.30  in  Urick  (1975,  pp.  248)  could  be  used  as  an  approxima¬ 
tion  for  shallow  waters .  Depending  on  three  generic  types  of 
bottom,  the  backscattering  Bs  is  calculated  using  a  third 
order  polynomial  fit  of  the  data  in  figure  8.30  as  a  function 
of  frequency  as  follows 

•  TYPE  I  :  Plain  bottom,  little  or  no  roughness,  high 
angular  and  frequency  variations. 

BS=0 . 012£3+o  .  283f  2+0 . 0062f-28 . 153  (2.41) 

•  TYPE  II:  Intermediate  between  a  very  rough,  rocky 
bottom  and  a  Type  I  bottom 
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BS=0 . 076/3_o  ,233f2+2 . 87  06f-40 .41 


(2.42) 


•  TYPE  III:  Heavily  dissected  bottom,  with  underwater 
ridges,  very  rough.  Almost  no  angular  or  frequency  dependence, 

B^=-18  (2.43) 

Figures  2 . 12  and  2 . 13  show  the  data  and  the  corresponding 
polynomial  fit  for  bottom  types  I  and  II.  The  user  must  decide 
which  type  of  bottom  to  use. 


Bottom  Backscattering  -  Type  I 


m 

Data 


Polynomial  fit 
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Scattering  coefficientidB 


Bottom  Backscattering  -  Type  II 


Pioruire  2 . 13 .  The  bottom  backscattering  for  intermediate  surfaces . 
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III.  THE  IMPLEMENTATION  OF  THE  MODEL 
A.  PROGRAM  STRUCTURE 

The  analytical  model  described  in  Chapter  II  was  imple¬ 
mented  in  a  structural  procedural  program  written  in  Turbo 
Pascal  7.0  for  DOS  (Borland,  1991)  .  The  .structure  of  the 
program  can  be  observed  in  Figure  3.1. 

The  program  consists  of  a  main  program  named  RTCAS 
(Reactive  Target  Computer  Assisted  Search)  and  seven  units 
namely,  InterManager,  ScoutManager ,  TargtManager, Detection 
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The  main  program  controls  the  sequence  in  which  each  one 
of  the  "Managers"  will  execute  their  tasks. 

A  single  data  structure  contains  all  the  information 
about  the  simulated  tracks,  the  searchers,  the  environment  and 
the  simulation  control  required  to  run  the  program. 

All  units  share  and  perform  their  tasks  in  a  common  data 
structure.  The  units  could  be  changed  provided  their  inputs 
and  outputs  remain  the  same,  thus  allowing  future  expansions 
or  modifications  in  the  software. 

The  implementation  of  the  Data  Fusion  Model  by  the  main 
program  RTCAS  is  described  in  the  algorithm  shown  in  Figure 
3.2.  The  code  of  RTCAS  is  in  APPENDIX  Al. 

RTCAS  general  algorithm 

Input  initial  simulation  data 
Input  initial  searchers  data 
Loop  until  selected  time 

Get  information  from  Searchers 
Loop  for  all  searchers 

If  searcher  is  in  contact 

Update  for  Positive  Information 

else 

Update  for  Neg-ative  Information 

end 

Normalize  Tracks'  Weights  with  Bayes  Theorem 

Update  Traks'  movements 

Update  cells  probabilities  with  Total  Probability  Law 
Output  Probability  Map 

FxgiiirG  3  •  2  •  Th.e  algroxithm  iTTiplemented.  in  the  main  program  of  RTCAS 
organizes  the  sequence  or  the  component  subprograms . 
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B.  DATA  STRUCTURE 


1 .  Tracks  Data 

The  number  of  tracks  was  fixed  at  500  for  computational 
convenience.  Processing  two  searchers  during  150  minutes  of 
activity  in  the  Local  Area  took  2  minutes  running  in  a 
Personal  Computer  with  a  486  DX2  66  MHZ  processor,  which 
implies  a  time  compression  ratio  of  1/75.  Any  future  increase 
in  computational  speed  may  permit  an  increase  in  the  number  of 
simulated  targets  in  order  to  augment  the  statistical  richness 
of  the  simulation. 

Each  track  is  represented  by  a  record  with  fields  for 
position,  course  and  speed  and  weight.  Information  common  to 
all  targets  is  stored  out  of  the  records  for  memory  space 
convenience,  e.g.,  standard  deviation  in  the  initial  posi- 
tions,  speed  and  course  ranges,  reaction  distance  and  reaction 
tuim  angle . 

All  the  infoinnation  concerning  the  tracks  is  input  by  the 
user  through  the  InterManager' s  corresponding  procedures. 

2.  Searchers'  Data 

The  program  can  process  up  to  15  searchers .  However  it  is 
very  unlikely  that  such  a  large  number  of  assets  will  be 
concentrated  for  a  Local  Area  or  Lost  Contact  search. 

Each  searcher  is  represented  by  a  record  with  fields  for 
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position,  course,  speed,  sonar  parameters  like  frequency, 
status,  source  level,  etc.  and  information  regarding  the 
contact . 

The  searchers'  data  structure  is  composed  of  static  data 
like  sonar  frequency,  detection  threshold,  etc.,  input  by  the 
user  and  controlled  by  the  InterManager  unit  and  by  dynamic 
data  like  position,  speed,  sonar  range  selected,  etc.  read 
from  a  file  by  a  procedure  in  the  ScoutManager  unit.  Reading 
from  a  file  resembles  receiving  the  dynamic  searchers' 
information  through  a  data  link  or  from  the  main  stream  of  a 
tactical  data  system. 

3.  Environmental  data 

Input  by  the  user  through  the  User  Interface  includes 
bottom  type,  sound  velocity  profile  type,  wind  speed  and 
depth.  This  part  of  the  data  structure  is  used  by  DetecManager 
in  the  determination  of  the  transmission  losses . 

4.  Simulation  control  data 

This  part  of  the  data  structure  stores  information 
required  in  the  control  of  the  program  like  the  name  and 
location  of  input /output  devices,  (files  in  this  implementa¬ 
tion)  ,  times  for  which  probability  maps  are  desired,  and  delay 
in  the  start  of  the  search  from  the  datum  time.  This  informa¬ 
tion  is  input  by  the  user  using  Intermanager  procedures. 
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C .  ALGORITHM  STRUCTURE 

This  section  outlines  the  algorithmic  structure  executed 
by  the  units  Interemanager,  Targe tManager,  ScoutManager, 
DetecManager  and  CellsManager  following  the  sequence  estab¬ 
lished  by  the  main  program  RTCAS  to  follow  the  Data  Fusion 
Model . 

1 .  Interface  Manager 

This  unit  declares  the  complete  common  data  structure, 
creates  the  user  input  or  static  part  of  that  structure  and 
provides  control  of  the  program  by  means  of  a  series  of 
interactive  menus  called  the  User  Interface. 

For  computational  convenience,  inside  the  program  the 
calculations  are  done  in  measurement  units  different  from 
those  input  by  the  user .  The  Interface  Manager  converts  the 
data  from  the  external  to  the  internal  unit  system,  i.e., 
discrete  time  unit  (dtu) ,  pixels,  and  discrete  speed  units 
(dsu)  ,  according  to  Table  3.1.  The  code  for  this  unit,  the 
file  INTERMAN. PAS  is  in  APPENDIX  B  and  the  sequence  of  menus 
is  included  in  APPENDIX  C. 


Internal  Units 

External  Units 

Time 

1  dtu 

5  minutes 

Distance 

1  pixel 

0.1  Nautical  Miles 
(185.2  meters) 

Speed 

rn.i.'U.l  ^  1 

1  dsu 

1.2  Knots 

Table  3.1.  The  measurement  units  conversion. 


51 


2 .  Target  Manager 

This  unit  implements  the  model  described  in  Chapter  II, 
Section  A  for  the  generation,  updating,  and  reaction  of  the 
simulated  targets.  Figure  3.3  depicts  the  process  to  decide  an 
evasion  course  at  the  target .  Note  that  an  evasion  course  is 
calculated  for  each  searcher  {evasion)  and  the  one  correspond¬ 
ing  to  the  closest  searcher  {d(s)<d(s-l))  is  selected  as 
course.  The  code  for  this  unit,  the  file  TARGTMAN . PAS ,  is  in 
APPENDIX  D. 
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3 .  Scout  Manager 

Implements  the  core  of  the  Data  Fusion  Paradigm  as 
described  in  Chapter  II  section  B. 

The  public  procedures  InitScouts,  CloseScoutsDataLink, 
ConnectWithScouts  and  ReadScoutDataLink  have  the  purpose  of 
controlling  the  input  of  information  from  the  searchers.  In  a 
real  world  implementation  this  procedure  would  receive  the 
information  through  a  Data  Link  device  or  from  the  main  stream 
of  the  Tactical  Data  System  of  the  platform.  In  RTCAS  the 

information  is  read  from  a  file  with  a  format  shown  in  Table 

3.2. 


The  file,  resembling  the  format  of  a  data  bus  in  a 
communications  link,  has  a  row  for  every  searcher  at  every 
time  updating  period,  fixed  in  the  program  to  5  minutes.  Each 
row  contains  information  regarding  the  identification  of  the 
searcher  as  a  ship,  helicopter  or  sonobuoy,  position,  course 
and  speed,  sonar  state,  sonar  range  scale  selected  and 
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contact  information  like  position,  course,  speed  and  credibil¬ 
ity.  Course  and  speed  are  used  to  keep  track  of  searchers  as 
well  as  contacts  and  also  used  in  the  calculations  of  the 
transmission  losses  by  DetecManager .  For  example,  in  Table 
3.2,  one  surface  ship  and  two  helicopters  are  in  the  Local 
Area.  At  time  25  the  ship  SI  is  emitting  but  does  not  have  a 
contact,  whirlwind  HI  is  dipping  with  its  sonar,  holds  a 
contact  and  assigned  a  credibility  value  of  .  75 .  Helicopter  H2 
is  in  transit  at  90  knots.  The  inclusion  of  the  sonar  condi¬ 
tions  in  the  "data  link  message"  makes  the  program  generic  for 
any  type  of  platform.  All  that  is  needed  to  fuse  the  informa¬ 
tion  with  the  simulation  is  position  and  whether  or  not  the 
sonar  is  on  or  not.  The  identification  at  the  beginning  of  the 
row  is  required  to  display  an  appropriate  symbol.  The  value  of 
credibility  may  be  considered  as  analogous  to  the  standard 
confidence  level  in  the  Allied  doctrine  (NATO,  1990) . 

The  procedure  ReadScoutsDataLink  updates  all  fields  in 
the  searcher  data  base  and  when  a  contact  is  produced,  the 
corresponding  fields  for  contact  information. 

The  basic  structure  of  the  Data  Fusion  Paradigm  performed 
by  the  ScoutManager  is  shown  in  Figure  3.4. 

The  code  for  ScoutManager  is  in  files  SCOUTMAN.PAS  and  is 
shown  in  APPENDIX  E. 
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4.  Detection  Manager 

The  sole  purpose  of  this  unit  is  the  calculation  of  the 
value  of  cdp5,  the  five  minutes  cumulative  detection  probabil¬ 
ity  as  developed  in  Chapter  2,  Section  C.  The  code  for  this 
unit  is  in  the  file  DETECMAN.PAS  and  is  shown  in  APPENDIX  F. 

5 .  Cells  Manager 

This  unit  renders  the  Probability  Map.  The  final  product 
of  the  procedure  UpdateProbinCell  is  a  bidimensional  array 
with  a  probability  value  associated  with  each  cell.  This 
product  is  the  probability  map.  It  could  be  used  in  a  differ¬ 
ent  kind  of  graphics  interface  or  transmitted  by  the  procedure 
BroadcastProbabilityMap.  Figure  3.1  shows  this  output  from 
CellsMan.  In  RTCAS  the  procedure  outputs  its  broadcasting  to 
a  file  that,  in  a  real  world  application,  would  be  replaced  by 
an  Input/Output  device.  The  purpose  of  broadcasting  is  to 
provide  the  same  information  about  the  probability  map  to  all 
searchers  and  other  recipients . 

The  code  for  this  unit  is  in  the  file  CELLSMAN. PAS  and 
can  be  seen  in  APPENDIX  G. 

6 .  Graph  Manager 

This  unit  produces  the  final  output  of  the  program  in  the 
form  of  a  color  coded  Probability  Map.  Additionally,  it 
produces  the  rest  of  the  Graphic  Information  Display  (GID) 
elements  like  grids,  color  scale,  static  and  dynamic  informa¬ 
tion  display,  etc.  Figure  3.5  show  the  layout  of  the  GID. 
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The  colors  in  the  probability  map  are  calculated  by  the 
procedure  SetProbabilityColors  using  the  color  scale  shown  in 
the  lower  left  corner  of  the  GID.  The  colors  are  relative  to 
the  largest  probability  value  on  any  cell. 

Instead  of  using  one  of  the  natural  color  scales  like  RGB 
or  RYB  (Elliot  et  al . ,  1995),  white  was  selected  as  the 
brightest  color  in  correspondence  with  the  maximum  probability 
value  because  of  the  black  background  RTCAS  uses. 

The  GID  shows  the  local  area  with  the  probability  map 
together  with  graphic  information  regarding  the  searchers 
activity.  The  symbology  found  in  most  Tactical  Data  Systems 
was  adopted  to  represent  the  different  types  and  conditions  of 
searchers  and  contacts. 

The  contact  information  is  displayed  with  a  hostile 
submarine  symbol  in  the  map  and  with  information  regarding  the 
time  of  initial  contact,  course,  speed,  searcher  holding  it, 
etc.  in  the  bottom  of  the  screen. 

The  code  for  this  unit  is  in  the  file  GRAPHMAN.PAS  and 
can  be  seen  in  APPENDIX  H. 
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Figore  3,5=  The  Graphic  Information  Display  (GID) 
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IV.  DATA  ANALYSIS  METHODS:  VERIFICATION  AND  EVALUATION  OF 
THE  MODEL 

The  Data  Analysis  Methods  described  in  this  chapter 
numerically  verify  the  correctness  in  the  implementation  of 
the  Bayesian  Filtering  Model  and  evaluate  the  significance  of 
using  a  Reactive  Target  Model  as  compared  to  using  a  Non 
Reactive  one. 

The  final  product  of  the  model  is  a  probability  map, 
hence  the  correctness  of  the  model  implementation  as  well  as 
the  relative  efficiency  of  different  methodologies  must  be 
evaluated  on  the  model's  ability  to  accurately  construct  the 
probability  distribution  of  the  location  of  the  target. 

The  verification  consists  of: 

•  Selecting  a  "dummy,"  assumed  to  be  the  real  target, 
among  the  simulated  ones. 

•  Running  several  scenarios  with  different  sets  of 
variables . 

•  Extracting  the  test  statistics  of  interest. 

•  Calculating  the  Measures  of  Effectiveness  (MOE) 
based  on  those  statistics 

•  Executing  a  test  to  verify  a  given  hypothesis. 

A.  MEASURES  OF  EFFECTIVENESS 

Accuracy,  Area  of  Uncertainty,  and  Mean  Missed  Distance 
were  selected  as  MOE's  to  respectively  verify  the  correctness, 
dispersion  and  location  of  the  distribution  represented  by  the 
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probability  map  and  the  significance  of  a  Reactive  Target 
Model  with  respect  to  a  Non  Reactive  one. 

Mean  Detection  Probability  and  Mean  Time  to  Detection 
were  selected  to  evaluate  the  performance  of  the  search 
conducted  with  a  Reactive  Target  Model  as  compared  to  a  Non 
Reactive  one  when  the  information  provided  by  the  probability 
map  was  used  in  conducting  a  search.  The  rest  of  this  section 
describes  each  one  of  the  selected  MOE's. 

1 .  Accuracy 

Accuracy  is  a  statistical  measure  of  how  well  the 
probability  map  represents  the  model's  uncertainty  in  the 
target's  location  (Widdis, 1995,pp.  14).  Furthermore,  accuracy 
is  a  measure  of  correctness  in  the  implementation  of  a 
Bayesian  Updating  algorithm.  The  containment  statistic  X  is 
the  cumulative  probability  of  all  cells  with  a  higher 
likelihood  than  that  of  the  cell  G  where  the  dummy  is. 


X  X 

(i,j)eK  {i.j)eK' 


(4.1) 


where 

K  :{  (i,j):  P(Gij\l)  >  P(G\I)  }, 

K':{  P(G,j\l)  =  P(G\I)  }, 

and  U  is  a  uniformly  distributed  random  variate  over  the 
interval  [0,1] .  The  second  term  in  (4.1)  represents  a  uniform 
proportion  of  the  all  cells  where  there  is  a  tie,i.e.,  the 
cell  probability  is  equal  to  that  of  the  cell  where  the  dummy 
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is.  Assuming  a  discrete  extension  of  the  Probability  Integral 
Transformation  (Dudewicz  et  al.,  1988,  pp674)  already  used  by 
Stone  et  al{1991,pp.  12)  and  Widdis  (1995, pp.  15) ,  the  random 
variable  X  should  be  uniformly  distributed  over  the  interval 
(0,1) . 

The  MOE  described  in  this  Section  is  a  discrete  version 
of  a  one  sample  Kolmogorof -Smirnoff  Test  Statistic  (Dudewicz 
et  al . ,  1988, pp.  671).  The  K-S  statistic  represents  the 

maximum  deviation  of  the  empirical  distribution  of  X  from  the 
Uniform, 

D=supre7nuin^g[o  |x:-F(x)  I  (4.2) 

where  the  second  term  is  the  sample  cumulative  distribution 
function  of  X.  Accuracy  is  defined  as  (Stone,  1991, p.  12) , 

A=l-D  (4.3) 

It  varies  between  0  and  1,  the  larger  the  value  the  more 
accurate . 

Performing  the  K-S  test  not  only  provides  a  measure  of 
correctness  of  the  model  by  means  of  the  accuracy  but  also 
permits  a  pictorial  evaluation  of  the  behavior  of  the  model. 
If  the  percentiles  of  the  sample  CDF  of  are  plotted  against 
the  percentiles  of  a  unifoinn  distribution,  the  discrepancy 
between  both  becomes  apparent  for  all  values  in  the  interval 
[0,1].  The  nature  of  those  discrepancies  characterize  a 
tracker  as  pessimistic  or  optimistic. 

optimistic  tracker  constructs  a  "concentrated" 
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probability  map,i.e.,  it  is  optimistic  about  its  own  ability 
to  determine  the  location  of  the  target.  As  a  consequence,  the 
position  of  the  dummy  target  is  more  often  than  not  far  from 
to  the  center  of  the  empirical  distribution,  thus  yielding 
much  more  frequently  larger  values  of  the  containment 
statistic  X  (Figure  4.1). 

A  pessimistic  tracker  yields  a  more  "dispersed" 
probability  map,  hence  the  frequency  of  lower  values  of  the 
containment  statistic  is  larger.  For  a  given  location  of  the 
target,  the  containment  statistic  X  will  be  stochastically 
larger  in  an  optimistic  tracker  than  in  a  pessimistic  one.  In 
Figure  4.1,  assume  that  both  the  pessimistic  and  the 
optimistic  trackers  represent  the  same  situation.  The  value  of 
the  cumulative  probability  or  "probability  mass"  inside  any 
given  equi-area  contour  in  the  optimistic  tracker  would  be 
greater  than  that  in  the  pessimistic  one.  As  a  consequence, 
if  a  contact  were  produced,  for  instance,  at  any  point  in  such 
a  contour, e .g. ,  the  point  indicated  with  a  chevron  in  both 
maps,  then  the  containment  statistic  X  would  tend  to  be 
greater  in  the  optimistic  tracker  than  in  the  pessimistic  one 
(0.9  and  0.2  respectively).  Such  a  tendency  translates  into 
the  shape  of  the  relative  cumulative  frequency  distribution  of 
X.  The  example  in  Figure  4.1.  shows  that  a  given  intermediate 
value  of  X  ,  for  instance  0.6,  is  much  more  frequent  in  a 
pessimistic  tracker,  i.e. ,  82%  of  the  time,  than  in  an 
optimistic  one, i.e.,  18%  of  the  time. 
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IViitaiBmMit  Sttiktie 


Figure  4.1.  Examples  of  Pessimistic  and  Optimistic  Trackers. 


2.  Area  of  Uncertainty  (AOU) 

The  Area  of  Uncertainty  (AOU)  can  be  defined  as  the  area 
of  those  cells  (i,j)  where  P(Gij\I)  >  P(G\I), 

AOU=  52 

{z,j)eK 


where  is  the  area  of  an  individual  cell  and  K  is  defined 
as  in  4.1.  The  value  of  AOU  provides  a  quantitative  idea  of 
how  spread  the  probability  map  is.  The  smaller  the  AOU  the 
better  since  the  probability  mass  is  more  concentrated,  thus 
leading  to  a  better  localization. 
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It  may  be  argued  that  the  AOU  and  the  accuracy  signify 
the  same  MOE,  but  it  is  clear  that,  for  example,  a  constant 
value  of  accuracy  may  coexist  with  different  values  of  AOU. 

3 .  Mean  Missed  Distance (MMD) 

The  MOE  selected  to  verify  the  location  features  of  the 
probability  map  with  respect  to  the  dummy  is  the  Root  Mean 
Square  Missed  Distance,  which  can  be  expressed  using  a 
straightforward  extension  from  the  one  used  in  Stone  et  al. 
(1988,  pp.  14) , 


N 


RMSMD= 


7=1 


(4.5) 


where  Pi=  P(Ti\I)  is  the  weight  of  each  simulated  target,  N  is 
the  total  number  of  simulated  targets  and  is  the  distance 
between  the  dummy  d  and  each  simulated  target  i.  This  MOE 
represents  the  standard  error  in  the  location  of  the  target. 

4.  Mean  Detection  Probability  (MDP) 

The  Mean  Detection  Probability  is 


DR=-^  (4.6) 
N 

where  is  the  number  of  times  out  of  N  identical  replica¬ 
tions  of  an  identical  experiment  in  which  a  searcher  obtains 
at  least  one  contact  within  a  period  of  time  T. 


5.  Mean  Time  to  Detection  (MTD) 

The  Time  to  Detection  TD  is  the  period  of  time  from  the 
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beginning  of  a  search,  i.e.,  time  zero  of  the  simulation  when 
the  searcher  enters  the  Local  Area,  to  the  time  at  which  the 
first  detection  is  produced  or,  if  no  contact  is  gained,  the 
end  of  the  experimental  scenario .  Then  the  Mean  Time  to 
Detection  is  the  average  among  N  replications  of  the  experi¬ 
ment  . 


B.  VERIFICATION  OP  THE  MODEL  IMPLEMENTATION  CORRECTNESS 

1 .  Test  Structure  and  Experiment  Design 

The  purpose  of  this  experiment  was  to  verify  whether  the 
implementation  correctly  models  a  Bayesian  filtering  process. 
The  MOE  or  rather  the  measure  of  correctness  was  the  Accuracy. 

For  this  experiment  a  dummy  target  was  selected  at  random 
among  the  500  simulated  targets.  The  Accuracy  was  measured 
with  respect  to  the  location  of  the  dummy  target .  The 
detections  were  sampled  based  on  the  five  minute  cumulative 
detection  probability  (cdpg) ,  corresponding  to  each  pair 
searcher-dummy . 

The  model  was  tested  using  Updating  for  Negative  informa¬ 
tion  only  and  for  Negative  and  Positive  information.  In  the 
former  case  the  simulation  ceased  at  the  first  contact .  The 
Reactive  Target  model  was  tested  in  combination  with  each 
updating  type. 
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Six  different  scenarios  or  searcher  patterns  were  played 
as  portrayed  in  Figure  4.2.  In  the  first  two  the  searcher  is 
stationary  resembling  sonobuoys  or  dipping  sonars,  in  the 
second  pair  the  searcher  moves  through  the  datum  area  in  a 
straight  line,  and  in  the  fifth  and  sixth,  one  and  two 
searchers  respectively  proceed  with  a  random  walk  type  of 
search  around  the  datum. 

The  number  of  replications  was  50  and  for  each  one  20 
time  samples  were  taken.  The  K-S  test  requires  a  random 
sample, i.e.,  the  observations  should  be  independent  and 
identically  distributed  (Dudewicz  et  al.,1988,  pp.280) .  In 
this  case,  the  1000  time  samples  are  not  independent.  Though 
only  one  sample  could  have  been  taken  from  each  replication, 
^11  time  samples  were  retained  for  they  provided  smoothing 
within  a  replication  "batch" . 

The  tests  were  identified  with  a  3  digit  label .  The  first 
digit  in  the  test  label  indicates  whether  Updating  for 
Positive  information  took  place,  1  if  yes,  2  if  not.  The 
second  digit  denotes  the  Target  Speed,  1  for  slow  (2-5  kts) 
and  2  for  fast  (6-12  kts)  the  last  digit  corresponds  to  the 
scenario. 

The  test  structure  is  shown  in  Table  4.1  and  the  set  of 
default  values  for  the  variables  used  in  the  experiment  are 
shown  in  Table  4.2. 
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NEGATIVE  INFORMATION 
UPDATING  ONLY 


POSITIVE  AND  NEGATIVE 
INFORMATION  UPDATING 


Table  4.1.  The  Experiment 
implementation  of  the  Bayesi 
the  Dummy  are  reactive. 


SLOW  TARGET 
(2-5  KTS)  . 

11 

12 

13 

14 

15 

16 


11 

12 

13 

14 

15 

16 

tructure  for  testing  the  correctness  in  the 
i  Updating.  In  all  Tests  the  Target  Model  and 
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VARIABLES 

VALUE 

TARGET  VARIABLES 

Datum  error 

2  nm 

Course  Range 

0-360° 

Minimum  Estimated  Target  Speed 

3  kts 

Maximum  Estimated  Target  Speed 

7  kts 

Reaction  Range 

5  nm 

Reaction  Tack 

150° 

SEARCHER  AND  SONAR  VARIABLES 

Source  Level 

200  Db  re  /iPa 

Ambient  Noise  Level 

80  Db  re  juPa 

Detection  Threshold 

20  Db  re  /iPa 

Directivity  Index 

30  Db  re  jiVa. 

Frequency 

5  Khz 

Pulse  Length 

0 . 05  sec 

Beam  Width 

12° 

Circular  Error  in  Contact  Reports 

1  nm 

Overall  Standard  Error  in  Sonar  Terms. 

6  Db  re  /uPa 

Credibility. 

1 

ENVIRONMENT  VARIABLES 

Type  of  Bottom 

Flat , Sandy 

Sound  Speed  Profile 

Isothermal 

Depth 

100  m 

Wind 

10  kts 

TIME  VARIABLES 

Total  Simulated  Search  Time 

100  min 

Scouts  Delay 

15  min 

4*2.  The  default  values  of  variables  used  in  the  Experiment  to 
test  the  correctness  of  the  model. 
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The  correctness  in  the  implementation  of  the  model  for 
each  scenario  was  decided  using  a  Kolmorgorv- Smirnov  One 
Sample  Test.  The  null  hypothesis  regarding  the  sample  CDF  of 
the  containment  test  statistic 

H^iF{x)=x  (4.8) 

tested  for  all  x  against, 

H^iF(x)*x 

is  rejected  at  level  a  if 

^n>^n.a  (4.10) 


where  D„  is  calculated  with  (4.2)  and  is  the  asymptotic 
critical  point  corresponding  to  the  sample  size  n,  obtained 
from  Dudewicz  et  al.  (1988, pp. 670) . 

The  accuracy  for  each  test  was  calculated  and  plotted  for 
a  visualization  of  the  pessimistic  or  optimistic  nature  of  the 
resulting  probability  map  and  then  logged  in  a  table. 

2.  Results 

Table  4.3  displays  the  results  of  this  test.  The  first 
two  columns  describe  the  tests  regarding  the  type  of  updating 
that  was  allowed,  the  target  model  in  use  and  the  type  of 
scenario  played. 

The  third  column  is  the  sample  size  n  used  in  performing 
the  tests.  The  variations  are  due  to  the  fact  that  in  the  case 
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of  updating  for  negative  information  only,  the  simulation 
stopped  at  the  first  contact.  This  happened  in  a  random 
fashion  in  time  hence  the  number  of  collected  samples  for  each 
replication  was  random  too. 

The  fourth  column  shows  the  average  nature  of  Optimistic, 
Pessimistic  or  Variable  in  the  probability  map  produced  at 
every  time  step  in  accordance  with  Figure  4.1. 


ID 

# 

UPDATING 
NEG  POS 

SAMPLE 
SIZE  n 

OPT 

PES 

KS 

TEST 

ACCURACY 

11 

Y 

N 

173 

VAR 

arr 

0.8919 

12 

Y 

N 

879 

OPT 

rrr 

0.9343 

13 

Y 

N 

671 

VAR 

aaa 

0.9625 

14 

Y 

N 

918 

VAR 

arr 

0.9538 

15 

Y 

N 

189 

VAR 

arr 

0.9052 

16 

Y 

N 

98 

PES 

aar 

0.8767 

21 

Y 

Y 

1000 

VAR 

rrr 

0.9343 

22 

Y 

Y 

1000 

OPT 

arr 

0.9269 

23 

Y 

Y 

1000 

VAR 

rrr 

0.9684 

24 

Y 

Y 

1000 

VAR 

aaa 

0.9842 

25 

Y 

Y 

1000 

VAR 

rrr 

0.9243 

26 

Y 

A  rj 

Y 

1000 

VAR 

aar 

0.9123 

Table  4  •  3  .  The  in^lementatioii  correctness  test  results . 


The  fifth  column  in  Table  4.3  shows  the  result  of  the 
Kolmorgorof-Smirnoff  test  of  hypothesis  defined  in  the 
previous  section.  The  three  letters  indicate  the  acceptance  of 
rejection  of  the  null  hypothesis  for  three  different  levels  : 
0.01,  0.05  and  0.10.  The  null  hypothesis  is  accepted  most  of 
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the  time  for  the  first  level  and  rejected  by  a  slim  margin, 
0.0152  in  the  average,  most  of  the  time  for  the  other  two. 

The  last  column  shows  the  final  values  of  Accuracy 
obtained  for  each  test.  The  results  range  from  0.8214  to 
0.9842  with  an  average  of  0.8945. 

C.  VERIFICATION  OF  THE  IMPLEMENTATION  ROBUSTNESS 

1.  Test  Structure  and  Esqperiment  Design 

A  sensitivity  analysis  was  performed  on  Accuracy  as  a 
function  of  Target  Speed,  Searcher  Speed,  Datum  Age,  Source 
Level,  Sample  Size  and  Credibility.  The  purpose  was  to  have  an 
idea  of  the  range  within  which  the  model  remained  correct .  For 
this  experiment,  with  the  exception  of  the  variable  of 
interest,  the  conditions  in  test  25  (Table  4.2)  were  used  for 
all  tests. 

2 .  Results 

Figure  4.3.  shows  the  plots  of  Accuracy  as  a  function  of 
Target  Speed,  Searcher  Speed,  Datum  Age  and  Source  Level.  No 
conspicuous  dependency  can  be  observed  within  the  explored 
ranges  of  each  one  of  those  variables,  and  accuracy  remains 
high. 

An  increase  in  Accuracy  as  a  function  of  Sample  Size  can 
be  observed  in  Figure  4.4.  This  relationship  conveys  the  idea 
that  the  quality  of  the  tests  is  noticeably  affected  by  the 
sample  size. 
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D.  THE  INFLUENCE  OF  CREDIBILITY 

1.  Test  Structure  and  Experiment  Design 

A  sensitivity  analysis  was  also  performed  on  Accuracy  as 
a  function  of  Credibility  since  the  latter  is  one  of  the 
hardest  parameters  to  estimate.  The  purpose  was  to  have  an 
idea  of  the  range  within  which  the  model  remained  correct.  In 
this  experiment  the  conditions  in  test  25  (Table  4.2)  were 
used,  except  for  the  credibility  Q.  In  the  tests  performed  in 
Section  B,  credibility  was  set  to  be  equal  to  1.  In  RTCAS,  the 
user  is  allowed  to  input  other  values.  This  experiment  tested 
values  of  Q  between  0.6  and  1.  Since  in  this  test,  as  in 
Sections  B  and  C,  the  contacts  were  generated  only  at  the 
position  of  the  dummy,  i.e.,  there  were  not  false  alarms,  the 
true  value  of  credibility  in  the  simulation  was  always  1 . 

2 .  Results 

As  can  be  observed  in  Figure  4.5,  the  Accuracy  curves 
were  plotted  for  several  values  of  Credibility  larger  than  0.5 
since  smaller  values  are  very  unlikely  to  be  used.  When  Q 
decreases  the  Accuracy  decreases  as  well.  A  tendency  towards 
a  pesimistic  probability  map  output  can  be  seen  as  decreasing 
values  of  Q  causes  the  model  to  ignore  good  information  and 
produce  a  less  "concentrated"  probability  map. 
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Aoculacy  Plots  fot  Vatbiis  CtedlbiUty  Values 


Figure  4.5.  Decreasing  the  value  of  Credibility  makes  the  probability 
map  more  pessimistic. 


E.  THE  SIGNIFICANCE  OF  A  REACTIVE  TARGET  MODEL 


1.  Test  Structure  and  Experiment  Design 

The  purpose  of  this  experiment  was  to  test  whether  or  not 
a  Reactive  Target  Model  is  operationally  better  than  a  Non 
Reactive  one.  All  of  the  MOE's  defined  in  section  A  were  used 
in  order  to  compare  the  results  from  both  Target  Models. 

The  dummy  was  selected  among  the  five  hundred  simulated 
targets  by  means  of  a  uniform  variate  drawn  and  was  set  to  be 
always  reactive.  This  implies  the  assumption  that  a  real 
target  is  reactive.  All  MOE's  were  measured  with  respect  to 
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the  location  of  the  dummy  target.  The  detections  were  decided 
as  in  Section  B.  The  model  always  updated  for  negative 
information  only,  stopping  at  the  first  contact.  The  purpose 
of  not  updating  for  positive  information  is  that  a  contact 
report  will  noticeably  improve  most  of  the  MOE's  due  to  the 
"concentration"  of  the  probability  mass  in  the  surroundings  of 
the  target  thus  obscuring  the  difference  between  a  Reactive 
and  a  Non  Reactive  Target  Model.  Furthermore,  the  differences 
between  both  Target  models  are  more  pivotal  to  the  problem 
before  the  first  detection. 

The  Reactive  Target  Model  was  tested  in  the  same 
conditions  as  the  Non  Reactive  Target  Model.  Two  different 
ranges  of  target  speed  and  three  different  values  of  datum 
delay  were  tested  as  can  be  observed  in  the  test  structure 
shown  in  Table  4.4.  The  main  purpose  of  the  test  was  to 
compare  how  good  each  model  (reactive  or  non  reactive)  deals 
with  the  same  scenario  when  constructing  the  probability  map. 
As  a  consequence  variables  with  the  greatest  impact  in  the 
kinematics  of  the  problem  were  selected. 

The  number  of  replications  was  the  same  as  in  the 
experiment  in  Section  D.  The  tests  were  labeled  with  a  3  digit 
code.  The  first  digit  in  this  case  indicates  the  speed  range, 
the  second  the  delay  and  the  last  one  the  scenario  or  search 
pattern,  which  are  the  same  as  in  the  test  in  Sections  D  and 
E  with  the  addition  of  a  seventh  pattern  representing  the 
first  270  degrees  in  an  exhaustive  spiral  search.  The  set  of 
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default  values  for  the  variables  used  in  the  experiment 
remained  the  same  as  in  Table  4.2  with  the  exception  of  the 
target  speed  and  delay  time  which  were  varied  from  test  to 
test.  For  each  running  all  MOE's  were  recorded  in  separate 
files  for  the  reactive  and  non  reactive  case  and  then  compared 
by  means  of  a  Wilcoxon  Test  for  paired  data  (Stone  et  al., 
1991,pp.28) .  The  Wilcoxon  Rank-Sum  Test  (Devore, 1991, pp. 612) 
is  a  non  parametric  test  used  to  compare  paired  observations 
without  regard  of  the  distributions  they  come  from. 


Target  Speed/ 

Delay 

2-5  kts 

5-11  kts 

0  min 

111 

211 

112 

212  , 

113 

213 

114 

214 

115 

215 

116 

216 

15  min 

121 

221 

122 

222 

123 

223 

124 

224 

125 

225 

126 

226 

127 

227 

30  min 

131 

231 

132 

232 

133 

233 

134 

234 

135 

235 

136 

236 

137 

237 

Table  4.4.  Test  structure  for  comparing  the  Reactive 
Target  Model  with  the  Non  Reactive  one.  Each  one  of 
the  tests  in  the  table  was  performed  for  each  model . 
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Let  Rj,  and  be  the  corresponding  MOE's  for  the  case  of 
a  Reactive  and  a  Non  Reactive  Target  Model  respectively  for 
each  run  i=I.  .n.  We  want  to  compare  whether  or  not  one  of  them 
has  a  tendency  to  be  larger  than  the  other,  and  are 
correlated  for  each  i,  but  the  differences  1?^=  R^-N^  are 
independent  for  all  i;  the  Wilcoxon  test  makes  their  original 
distributions  irrelevant  to  the  comparison.  A  statistically 
valid  statement  about  such  a  comparison  can  be  done  by  means 
of  a  test  of  hypothesis  of  the  mean  [1^  and  of  each  r.v. 
where 

(4.11) 

and, 

(4.12) 

are  the  null  and  alternate  hypotheses.  The  order  in  the 
differences  between  both  mean  values  can  be  varied  depending 
on  the  nature  of  the  MOE.  For  instance  if  R^  and 
represented  the  Accuracy,  we  would  want  to  know  which  is 
larger  since  the  larger  the  better  in  such  a  case.  Conversely 
if  they  were  the  Area  of  Uncertainty,  we  would  want  to  find 
out  which  is  smaller. 

In  this  test,  the  differences  =  (R^  -  N^)  are  obtained 
and  sorted  by  ascending  absolute  value.  The  rank  of  such  a 
sorting  is  r^.  Adding  up  the  rank  of  all  negative  yields  the 
test  statistic 
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(4.13) 


Di>0 

If  there  is  not  a  tendency  of  the  random  variables  R  to 
be  smaller  than  N,  then  W  is  expected  to  be  roughly  close  to 
the  value  of  its  mean  =  n(n+l)/4,  i.e.,  half  of  the  sum  of 
the  ranks  from  1  to  n.  Conversely,  if  tends  to  be  smaller, 
the  test  statistic  F\r  will  be  smaller  than  n(n+l)/4. 

For  n>30  IV  is  a  normal  random  variable  with  mean  and 
variance  a„‘  =  (n+1)  (2n+l) /24 .  As  a  consequence  the  test  of 
hypothesis  can  compare  the  value  of  W  with  the  critical  value 
of  a  normal  distribution  C(a,n)  corresponding  to  the  sample 
size  n  and  the  desired  significance  level  a.  If  IV  >  C(a,n) 
then  we  reject  the  null  hypothesis. 

2 .  Results 

Table  4.5.  sumarizes  the  results  of  the  Wilcoxon  Tests  of 
Hypothesis  conducted  to  verify  whether  or  not  using  a  Reactive 
Target  Model  improves  the  performance  of  the  tracker  from  the 
standpoint  of  Accuracy,  Area  of  Uncertainty  and  Mean  Missed 
Distance.  The  tests  were  performed  on  the  data  shown  in  Tables 
4.6  and  4.7  ,  where  for  each  MOE  and  each  run,  an  asterisk 
indicates  when  the  Non  Reactive  Target  Model  wins  for  a  visual 
appraisal  of  the  tendencies. 

On  the  average,  as  can  be  observed  in  Table  4.5,  second 
column.  Accuracy,  Area  of  Uncertainty  and  Mean  Missed  Distance 
show  a  better  result  for  a  Reactive  Target  Model.  However, the 
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Wilcoxon  Test  does  not  reject  the  null  hypothesis  that  there 
is  no  difference  between  the  RT  and  the  NRT  models. 
Furthermore,  the  null  hypothesis  should  be  accepted  for 
Accuracy  and  Area  of  Uncertainty.  This  result  is  somewhat 
disappointing  for  we  expected  the  former  to  be  the  MOE  where 
the  improvement  would  be  the  greatest.  Nonetheless,  at  a  10% 
significance/  would  be  rejected. 


The  RT  model  proved  to  be  definitely  better  than  the  NRT 
one  with  respect  to  the  Mean  Missed  Distance,  yet  by  a  slim 
margin . 


MEAN 

NRT 

MEAN 

RT 

W 

p-val- 

ue 

Ho 

TEST 
RESULT 
a=  1% 

ACC 

0.8292 

0.8567 

285 

451 

152 

0.2743 

Accep¬ 

ted 

AOU 

9.8121 

8.3223 

187 

451 

152 

0.0834 

Mr-Mn=0 

Accep¬ 

ted 

MMD 

7.0542 

5.8874 

50 

451 

152 

0.0084 

Rejec¬ 

ted 

Tcible  4.5.  significance  of  a  Reactive  Target  Model  in  the  average  ana  in 
a  Wilcoxon  Test  of  Signed  Ranks.  W  is  the  Wilcoxon  Test  Statistic. 


P.  FURTHER  EVALUATION  OF  THE  SIGNIFICANCE  OF  A  REACTIVE 
TARGET  MODEL:  A  NON  OPTIMAL  MYOPIC  SEARCH  EXPERIMENT 


1.  Test  Structure  and  Experiment  Design 

Although  the  experiment  described  in  Section  E  ascertains 
that  a  Reactive  Target  Model  performed  better  than  a  Non 
Reactive  one  in  most  of  the  cases,  the  difference  in 
effectivness  reflected  by  the  tests  is  much  less  significant 
than  the  expectations.  In  the  previous  set  of  tests  the 
searchers  proceeded  without  regard  to  the  probability  map 
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output  by  the  model .  They  followed  typical  lost  contact 
exhaustive  initial  search  patterns. 

In  order  to  further  evaluate  the  significance  of  a 
Reactive  Target  Model  as  compared  to  a  Non  Reactive  Target,  an 
experiment  was  conducted  in  which  the  searchers  were  set  to  be 
fed  back  by  the  output  of  the  model,  conditioning  their 
behavior  by  the  probability  map  rendered  at  every  time  step. 
At  all  times,  the  single  searcher  simply  headed  for  the 
brightest  cell. 

The  visualization  of  the  search  pattern  conducted  with 
the  above  mentioned  rules  was  often  irrational  from  a  tactical 
point  of  view  in  the  sense  that  the  searcher  would,  for 
instance,  move  from  one  brightest  cell  to  the  other  leaving 
unexplored  cells  that  were  close  to  his  way.  Again,  the 
purpose  was  to  have  a  standard  set  of  rules  of  behavior  of  the 
searcher  in  order  to  obtain  more  insight  in  the  comparison 
between  a  Non  Reactive  Target  Model  and  a  Reactive  one.  The 
dummy,  uniformly  drawn  from  the  total  number  of  simulated 
targets  at  each  replication,  was  set  to  be  always  reactive  and 
the  actions  were  set  to  last  until  a  contact  was  obtained  or 
150  minutes,  whichever  occurred  first.  The  detections  were 
decided  as  in  section  B.  The  set  of  tests  was  again  identified 
with  a  3  digit  code  as  shown  in  Table  4,8.  The  first  digit 
indicates  the  target  speed  range,  the  second  one  the  datum  age 
and  the  third  one  the  searcher  speed  and  initial  position. 
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ID  # 

ACC 

NRT 

RT 

111 

0.7523 

0.8500 

112 

0.9361 

0.9636 

113 

0.9295 

0.9353 

114 

0.9518 

0.9551 

115 

0.7795 

0.8541 

116 

0.7401 

0.8456 

117 

0.9089* 

0.8687 

121 

0.7134 

0.9635 

122 

0.9441* 

0.8242 

123 

0.9230 

0.9261 

124 

0.8035 

0.9012 

125 

0.8872 

0.9162 

126 

0.7731* 

0.6647 

127 

0.8810 

0.8623 

131 

0.6532 

0.7823 

132 

0.8452* 

0.8124 

133 

0.8687* 

0.8536 

134 

0.8049 

0.9304 

135 

0.7689 

0.8120 

136 

0.8897 

0.7025 

137 

0.8873 

0.9151 

avg 

0.8400 

0.8637 

AOU 


NRT 


4.4575  4.2770  3.4273 


5.9462  3.5910  2.5723 


3.9566 


6.5120  4.4159 


9.5123  4.6987 


3.6002 


2.7164 


5.4262  2.5634  3.4049 


2.8725  2.1560  2.8419  1.2563 


3.4521  3.8001  2.5800 


4.7425  4.9179 


0.8242  13.2612  10.2040  4.6058 


9.5623 


12.9650  6.0323  5.1757 


10.1212  8.5612  4.5867 


3.8525  2.3350  3.6402 


12.2564  5.3115 


5.1925*  5.5478  5.5896 


9.5412  8.4150 


9.9871  6.6959 


15.2412 


5.8650  5.6523  4.0254 


2.7475  4.6841  3.8912 


10.2547  5.9379  4.2356 


4.1629 


4.2763 


6.1347  5.4798 


5.3936 


2.8612 


4.7462 


4.5200 


5.9528  3.5600 


6.4589 


7.2086  6.0514 


7.2992 


4.95149  4.0451 


Table  4.6.  The  results  of  the  Coir^arison  between  a  Reactive  Target 
Model  (RT)  with  a  Non  Reactive  one  (NRT)  for  slow  target  speeds.  Asterisks 
indicate  when  NRT  wins. 
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ID  # 


ACC 


AOU 


MMD 


NRT 

RT 

NRT 

RT 

NRT 

- - 1 

RT 

211 

0.6947 

0.7942 

7.7325 

5.7474 

6.5314 

4.1247 

212 

0.896* 

0.7123 

13.0422 

5.4373 

6.5159 

2.8195 

213 

0.8624 

0.9237 

9.6524 

7.4512 

6.9158 

6.1266 

214 

0.893* 

0.8222 

12.2575 

8.5647 

7.9634 

4.9579 

215 

0.7242 

0.8520 

7.6672 

7.6542 

5.2720 

3.2549 

216 

0.733* 

0.7047 

3.3073 

2.0707 

4.0580 

3.1092 

217 

0.8962 

0.9014 

13.1071 

8.4475 

5.9064 

3.6117 

221 

0.7530 

0.8562 

12.5350 

9.4358 

10.0125 

8.6542 

222 

0.8444 

0.9115 

8.445* 

8.4512 

8.4972 

5.7633 

223 

0.9134 

0.9139 

13.2425 

9.5623 

9.8595 

9.1776 

224 

0.901* 

0.8425 

13.4175 

12.6275 

9.7686 

8.3196 

225 

0.8223 

0.8366 

11.8777 

7.5158 

8.8233 

7.1634 

226 

0.7219 

0.8675 

4.050* 

5.6475 

7.7192 

6.1558 

227 

0.926* 

0.9027 

12.025* 

12.9775 

9.8413 

6.2208 

231 

0.8125 

0.8457 

19.6851 

10.2125 

13.4794 

11.8224 

232 

0.7846 

0.9315 

6.267* 

8.6958 

10.9819 

10.2591 

233 

0.903* 

0.8728 

8.965* 

9.8561 

12.9581 

12.8511 

234 

0.8044 

0.8841 

9.795* 

15.6056 

12.9458 

12.7729 

235 

0.7337 

0.8412 

6.325* 

12.7345 

9.817* 

11.4833 

236 

0.7951 

0.8324 

15.1725 

7.4125 

12.1102 

12.0041 

237 

0.7711 

0.7958 

6.962* 

20.1525 

12.3220 

11.6747 

AVG 

0.81849 

0.84975 

10.26347 

9.345557 

9.157086 

7.729829 

4«7.  The  results  of  the  Cotnpeirison  between  a  Reactive  Target 
Model  (RT)  with  a  Non  Reactive  one  (NRT)  for  fast  target  speeds .  Asterisks 
indicate  when  NRT  wins . 
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TEST 

ID 

TARGET 

SPEED 

DATUM 

AGE 

SEARCHER 

SPEED 

111 

2-4  Kts 

10  min 

12  kts 

112 

2-4  Kts 

10  min 

18  kts 

121 

2-4  Kts 

60  min 

12  kts 

122 

2-4  Kts 

60  min 

18  kts 

211 

4-8  Kts 

10  min 

12  kts 

212 

4-8  Kts 

60  min 

18  kts 

221 

4-8  Kts 

10  min 

12  kts 

222 

4-8  Kts 

60  min 

18  kts 

311 

8-12  Kts 

10  min 

12  kts 

312 

8-12  Kts 

60  min 

18  kts 

321 

8-12  Kts 

10  min 

12  kts 

322 

8-12  Kts 

60  min 

18  kts 

Table  4.8.  The  non  optimal  myopic  search  experiment  structure. 


For  each  test  30  replications  were  run  and  the  average 
values  of  Detection  Rate,  Time  to  Detection  and  Search  Speed 
were  extracted. 

2 .  Results 

Table  4 . 9  shows  the  results  of  this  test .  The  Mean 
Detection  Probability  showed  the  strongest  difference  between 
a  Non  Reactive  model  and  a  Reactive  one,  favorable  to  the 
latter.  The  Wilcoxon  Test  of  Signed  Ranks  Sum  performed  on 
Mean  Detection  Probability  and  Mean  Detection  Time  yielded  the 
results  shown  in  Table  4.10.  At  a  5  %  significance  level,  the 
null  hypothesis  that  there  is  no  difference  between  a  Reactive 
Target  Model  and  a  Non  Reactive  one  can  be  Rejected  with 
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respect  to  the  Mean  Detection  Probabilities  and  should  be 
accepted  regarding  the  Mean  detection  Time. 

The  observation  of  the  evolution  of  the  search  step  by 
step  on  the  screen  output  of  the  model,  provides  more  insight 
about  the  real  gist  of  Mean  Time  to  Detection  in  the  context 
of  the  test.  In  the  Non  Reactive  Target  Model,  the  searcher 
promptly  starts  "consuming"  probability  mass  since  the 
expanding  ring  of  the  probability  map  does  not  change  because 


MDP 

MDT 

NRT 

RT 

NRT 

RT 

111 

0.700 

0.900* 

73.8 

54.3* 

112 

0.767 

0.933* 

28.1* 

28.3 

121 

0.267 

0.700* 

48.8 

39.6* 

122 

0.867 

1.000* 

55.6 

32.1* 

211 

0.067 

0.167* 

75.5 

67.5* 

212 

0.667 

0.667* 

45.1 

39.2* 

221 

0.167 

0.267* 

45.2* 

46.1 

222 

0.833 

0.900* 

61.3* 

63.3 

311 

0.000 

0.000 

150.0 

150.0 

312 

0.567 

0.633* 

1 

54.2* 

56.7 

321 

0.000 

0.000 

150.0 

150.0 

322 

0.267 

0.367* 

58.6 

39.6* 

AVG 

'P’bVvH  ^  a  a  ml- 

0.431 

0.486* 

45.5 

38.9* 

Table  4,9.  The  non  optimal  myopic  search  experiment  structure 
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MOE 

W 

Mw 

Ow 

p -value 

Ho 

Hi 

RESULT 

MDP 

0 

39 

15.7 

0 

Mr=Mn 

Rejected 

(5%) 

MDT 

18 

39 

15.7 

0.059 

Mr^Mn 

Accepted 

(5%) 

Table  4.10.  The  Reactive  Target  Model  performed  better  than  the  Non 
Reactive  one  with  respect  to  Mean  Detection  Time  but  no  difference  was 
appreciated  with  respect  to  Mean  Detection  Probabilities. 


of  the  searching  activity.  Therefore,  when  it  proceeds  towards 
the  brightest  cell,  it  tends  to  have  a  lot  of  simulated 
targets  within  its  detection  range.  This  is  reasonable  since 
the  relative  velocity  between  each  simulated  target  and  a 
searcher  will  be  high  because  of  the  average  outbound  nature 
in  the  course  of  the  latter  and  the  inbound  course  in  the 
former  (See  Figure  4.6) .  Conversely,  in  the  Reactive  Target 
Model,  the  presence  of  the  searcher  will  "deform"  the 
expanding  ring  of  the  probability  map  in  a  way  such  that  the 
relative  velocity  between  searcher  and  targets  will  be  low  in 
general,  i.e.,  a  chasing  situation  will  be  the  norm  most  of 
the  time.  As  a  consequence,  the  rate  of  consumption  of 
probability  mass  will  be  lower  than  in  the  Non  Reactive  case 
and  the  time  to  detection  will  be  larger,  thus  leading  to  a 
slower  Search  Speed.  (See  Figure  4.7) 

Coalescing  the  ideas  of  the  three  MOE's  analyzed  in  this 
section,  it  can  be  stated  that  using  a  Reactive  Target  Model 
yielded  a  slower  yet  more  effective  search.  Nonetheless  the 
margins  by  which  this  sentence  can  be  stated  are  rather  slim. 
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Figure  4o6o  The  searcher  "consumes"  probability  mass  to  the  N.E.  of 
the  expanding  datum  of  a  Non  Reactive  Target  Model.  The  darker  colors 
indicate  those  cells  where  the  likelihood  is  reduced  by  the  search. 
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Figure  4,7.  The  reaction  of  the  target  prevents  the  searcher  from  '‘consuming 
too  much  probability  mass  as  it  proceeds  with  the  search  pattern.  The  expan-. 
Sion  of  the  datum  is  somewhat  irregular. 


V.  EXTENSIONS  OF  THE  CONCEPT  AND  IMPLEMENTATION  ENHANCEMENTS 
SUGGESTED  FOR  FUTURE  RESEARCH 

A.  EXTENSIONS  OF  THE  CONCEPT 

The  Information  Fusion  Paradigm  (IFP)  can  be  the  engine 
of  several  problems  other  than  the  one  solved  in  this  thesis. 

1 .  Other  Sensors 

Extension  of  the  model  to  other  sensors  like  radar, 
lidar,  infrared,  etc.  is  straightforward.  Even  the  idea  of  a 
multisensor  search  is  plausible  provided  that  the  correspond¬ 
ing  models  to  evaluate  the  effectiveness  of  the  search  are 
implemented.  The  structural  design  of  the  program  allows  the 
interchange  of  the  sonar  detection  model  almost  without 
penalties . 

2 •  RTCAS  and  Java 

The  appearance  of  the  programming  language  Java  (Sun 
Microsystems  Corporation,  1994)  offers  a  new  dimension  to  the 
Information  Fusion  Paradigm.  Java  allows  relatively  easy 
networking  between  different  types  of  computers .  This  new 
level  of  interchange  of  information  not  only  favors  the 
handling  of  the  real  part  in  the  IFP,  but  also  brings  the 
possibility  of  a  Distributed  Simulation  (Buss,  et  al.,1996), 
as  an  option  to  deal  with  the  "virtual"  part  of  the  IFP.  Such 
a  connectivity  is  a  perfect  solution  when  for  instance  the 
search  is  conducted  with  heterogeneous  assets,  as  in  the  realm 
Search  and  Rescue  (SAR)  problem.  According  to  interna- 
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tional  regulations,  any  vehicle,  e.g.,  fishing  or  merchant 
vessels,  airplanes,  lorries,  etc.  can  be  convoked  to  partici¬ 
pate  in  a  search  by  a  SAR  agency.  Nowadays,  the  only  existing 
links  between  those  heterogeneous  assets  are  the  International 
Maritime  Radio  Circuits.  Given  that  most  of  those  platforms 
can  cariry  a  computer,  Java  brings  the  possibility  of  linking 
them  together.  Then,  it  is  possible  to  implement  a  Real  Time 
Bayesian  Filtering  in  which  each  one  of  the  "convoked" 
searchers  will  not  only  input  its  searching  activity  via,  for 
instance,  an  Internet  browser,  but  also  will  run  its  own 
Stochastically  Synchronized  Simulation,  i.e.,  the  seed  of  the 
random  generation  is  broadcasted  in  the  network.  Hence  each 
one  of  the  searchers  will  have  its  own  picture  of  the  situa- 
thus  leading  to  a  more  intelligent  allocation  of  effort. 

B .  IMPLEMENTATION  ENHANCEMENTS 

1.  Model  enhancements 

We  envision  a  lot  of  potential  growth  without  the  penalty 
of  very  expensive  research  and  development .  Some  ideas 
susceptible  to  further  research  are  outlined  in  this  section. 

A  Multiprior  target  generation  scheme  would  enhance  the 
^^flity  of  RTCAS  to  capture  reality  by  means  of  several 
weighted  options. 

An  Intelligent  Reactive  Target  Model  would  allow  each 
to  evaluate  more  factors  before  choosing  a  Reaction 
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Pattern.  One  of  those  factors  may  be  the  counterdetection 
estimated  range.  Solving  the  sonar  equation  for  the  one  way 
path  from  the  searcher  to  the  target  is  more  realistic  than 
using  a  fixed  distance  as  a  reaction  parameter.  Artificial 
Intelligence  tools  and  concepts  would  be  in  place  when 
developing  the  new  reactive  features. 

In  water  depths  in  the  transition  between  deep  and 
shallow  from  the  acoustic  point  of  view,  shadow  zones  can  be 
expected  within  several  cases  of  sound  velocity  profiles.  As 
a  consequence,  the  search  effort  is  not  uniform  over  the  water 
A.  possible  solution  to  this  problem  would  be  enhance 
RTCAS  with  the  addition  of  a  third  dimension  in  the  data 
structure.  The  probability  map  could  be  rendered  for  different 
layers  that  would  help  decide  how  to  allocate  the  search 

effort  not  only  in  time  and  xy  positions  but  in  depth  as 
well . 

2 .  Implementation  Enhancements 

Pascal  was  selected  as  a  programming  language  because  of 
its  robustness  and  maintenance  simplicity.  However,  RTCAS  may 
benefit  from  implementation  in  other  languages  like  C  or  C++ 
mainly  because  of  the  speed  provided  by  the  pointer  arithmetic 
so  useful  in  the  simulation,  or  Java,  because  of  the 
connectivity  that  comes  with  the  language.  The  decision 
should  be  done  only  after  evaluating  advantages  and  disadvan¬ 
tages  of  such  a  change  or  expansion. 
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Regardless  of  the  language,  RTCAS  must  migrate  to  a 
graphic  interactive  user  interface  in  order  to  allow  easier 
manipulation  of  its  possibilities.  The  more  the  interface  of 
the  program  resembles  other  popular  programs,  the  faster  the 
learning  curve  and,  as  a  consequence,  the  efficiency  in  the 
use  of  its  resources.  In  the  realm  of  the  user  interface,  the 
need  can  be  also  mentioned  for  zooming  and  scaling  on  the  fly. 
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VI.  CONCLUSIONS  2^  RECOMMENDATIONS 


A.  CONCLUSIONS 

RTCAS  correctly  implements  a  Bayesian  Filtering  Process 
for  Fusing  Information  from  a  real  world  active  search  with 
a  simulation  encompassing  many  possible  positions  and 
movements  of  a  reactive  target  to  produce  a  probability  map  of 
target  location. 

In  an  active  search,  the  target  will  invariably  react. 
Using  a  very  rudimentary  Reactive  Target  Model  proved  to  be 
better  than  a  non  reactive  one  in  most  of  the  cases. 

B .  RECOMMENDATIONS 

Intensive  testing  of  RTCAS  should  continue  to  further 
verify  its  behavior  and  to  validate  it  for  use  in  real 
environments . 

Extensions  of  the  Information  Fusion  Paradigm  should  be 
explored  and  enhancements  to  RTCAS  should  be  implemented  to 
increase  its  power. 

The  possibility  of  implementing  a  reactive  target  model 
in  currently  used  TDA's  should  be  explored. 
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APPENDIX  A.  RTCAS  MAIN  PROGRAM 


program  rtcas; 
uses 

G  RAPH,  CRT ,GraphManager,T argtManager,  DetecManager, 
CellsManager,Scoutmanager,lnterManager; 

begin 

cirscr; 

Graphmanager. Start; 

Scoutmanager.lnitscouts; 
while  WannaContinue  do  begin 
WannaSeeData:  =  FALSE; 
Scoutmanager.ConnectWithScoutsThruDataLink; 
Intermanager.AskOptions; 
Detecmanager.UpdateSonarParameters; 
if  WannaContinue  then  begin 
Targtmanager.GenTargets; 
Graphmanager.Pleasewait; 
ScoutManager.ResetUndetectedMass; 
if  WannaSeeData  then  begin 
Graphmanager.DisplayData; 
end  else  begin 

for  time:  =  1  to  totalsimtime  do  begin 
Scoutmanager.ReadScoutsDataLink; 
Scoutmanager.UpdatePriorWithScoutsInfo; 
Targtmanager.UpdateWeightsWithBayes; 
Scoutmanager.UpdateUndetectedMass; 
Cellsmanager.UpdateInCellProbabilities; 
Graphmanager. OutputGraph; 
end; 
end; 
end; 

Scoutmanager.CloseScoutsDataLink; 

end; 

end. 
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APPENDIX  B*  UNIT  INTERMANAGER 


unit  InterManager; 


interface 

< - 


uses 

CRT; 


^TARGETS  CONSTANTS. 
< - 


TYPES  AND  VARIABLES} 


const 

NUMBEROFTGTS  =500; 

BIGNUMBER  =1E30; 

DEFMOTIONMOD  ='OMNI' 

DEFSIGMARINIT  =20; 

DEFPORTMOSTCOURSE=0; 
DEFCOURSERANGE  =360; 

DEFMINSUBSPEED  =3; 

DEFMAXSUBSPEED  =7; 

DEFREACTRANAGE  =40; 

DEFREACTIONTACK  =150; 


OMNI 

type 

targetype 


='OMNr 


=  record 

X 

y 

course 
speed 
cel  lx 

celly  : integer; 

prob  :real; 

threatdis  :real; 

evas i oncour se : rea  L ; 


: integer; 
: integer ; 
:real; 
:real; 

: integer; 


datumx 

datumy 


end; 


: integer; 
: integer; 


CNUMBER  OF  TARGETS> 
OiTILITY  VARIABLE> 
^OMNIDIRECTIONAL  MOTION  M0DEL> 
CDATUM  ERROR> 
CLEFT  LIMIT  COURSE> 
CCOURSE  RANGE} 
CMIN  SUB  SPEED} 
CMAX  SUB  SPEED} 
CREACTION  RANGE} 
CREACTION  CHANGE  IN  COURSE} 
CMOTION  MODEL} 


CTARGET  STRUCTURE} 
CX  POSITION} 
a  POSITION} 
CTARGET  COURSE} 
CTARGET  SPEED} 
CTARGET  CELL  IN  X} 
CTARGET  CELL  IN  Y} 
CTARGET  INDIVIDUAL  WEIGHT} 
CTHREAT  DISTANCE  ACCUMULATOR} 
CTACK  ANGLE} 
CTARGET  MIND  DATUM  POSITION  IN  X} 
CTARGET  MIND  DATUM  POSITION  IN  Y} 


tgtbundletype 

=  arrayCI 

var 

mytarget 

ztgtbundl 

sigmarinit 

:byte; 

minsubspeed 

rbyte; 

maxsubspeed 

:byte; 

courserange 

: integer; 

portmostcourse 

linteger; 

motionmodel 

: St ring; 

reactrange 

: integer; 

react iontack 

: integer; 

totprob 

:real; 

CSET  OF  ALL  TARGETS} 


CSTANDARD  DEVIATION  OF  THE  DATUM} 
CESTIMATED  MIN  TARGET  SPEED} 
CESTIMATED  MAX  TARGET  SPEED} 
CESTIMATED  RANGE  OF  COURSES  FOR  THE  TARGET} 
CPORT  LIMIT  IN  THE  COURSE  RANGE} 
CTYPE  OF  MOTION  MODEL} 
CDISTANCE  OF  REACTION  FROM  SEARCHER} 
CANGLE  OF  EVASION  FROM  THE  BEARING  TO  THE  SEARCHER} 
CACCUMULATOR  OF  TARGETS  WEIGHTS  AFTER  UPDATING} 


97 


CSEARCHERS  CONSTANTS.  TYPES  AND  VARIABLES} 

C - 

const 

MAXSCOUT  =15; 

ANGLE  =360; 

NUMOFPULSES  =10; 

CDEFAULT  VARIABLE  VALUES} 

DEFSL  =200; 

DEFNL  =80; 

DEFDT  =20; 

DEFDI  =20; 

DEFABSCOEFF  =0-004; 

DEFREVRBAREA=8; 

DEFSCATTSTRG=-40; 

DEFCONT_.AOU  =40; 

DEFFREQUENCY=6; 

DEFPULSELENG=0.05; 

DEFBEAMWIDTH=14; 

DEFNOISESIG  =30; 

type 

scouttype  =  record  CSEARCHERS  DATA  STRUCTURE} 

X  : integer;  {SEARCHER  X  POSITION} 

y  : integer;  {SEARCHER  Y  POSITION} 


course 

: integer; 

{SEARCHER  COURSE} 

speed 

:real; 

{SEARCHER  SPEED} 

SONAR_ON 

: boo lean; 

{SONAR  CONDITION  ON  OFF} 

maxscale 

: integer; 

{SONAR  SCALE  SELECTED  IN  THE  SEARCHER'S  SONAR} 

IS  JOT 

:boolean; 

{WHETHER  THE  SEARCHER  HOLDS  A  CONTACT  OR  NOT} 

abscoeff 

:real; 

{ABSORPTION  COEFFICIENT} 

frequency 

:real; 

{SONAR  FREQUENCY} 

pulselength:real; 

{PULSE  LENGTH} 

reverbarea 

: integer; 

{AREA 

OF  REVERBERATION  FOR  GIVEN  LENGTH  AND  BEAMWIDTH} 

beamwidth 

:real; 

{SONAR  BEAMWIDTH} 

SL 

: integer; 

{SOURCE  LEVEL} 

NL 

:byte; 

CNOISE  LEVEL  CORRESPONDING  TO  THE  SEARCHER  SPEED> 

DT 

;byte; 

{DETECTION  THRESHOLD} 

DI 

:byte; 

{DIRECTIVITY  INDEX} 

SS 

: integer; 

{SCATTERING  STRENGTH  CORRESPONDING  TO  SONAR  C0NDITI0NS> 

subx 

:  integer; 

{X  POSITION  OF  A  CONTACT> 

suby 

:integer; 

{Y  POSITION  OF  A  CONTACT> 

subcourse 

: integer; 

{CONTACT  COURSE> 

subspeed 

: integer; 

{CONTACT  SPEED> 

cred 

:real; 

{CREDIBILITY  IN  THE  CONTACT  REPORT> 

id 

:char; 

{PLATFORM  IDENTITY,  SURFACE  AIR  BUOY> 

initialhot 

: integer; 

{TIME  OF  INITIAL  CONTACT> 

end; 

scoutteamtype  =  array [1  - -MAXSCOUT]  of  scouttype;  {STRUCTURE  TYPE  OF  ALL  SEARCHERS} 

var 

myscout  : scoutteamtype;  {INSTANCIATED  SEARCHERS} 

numberofscouts  :byte;  {TOTAL  NUMBER  OF  INSTANCIATED  SCOUTS} 

noisesigma  :byte;  {STANDARD  DEVIATION  IN  THE  SOLUTION  OF  SONAR  EQQUATION} 

undetectedmass  :real;  {UNDETECTED  PROBABILITY  MASS} 


{SOURCE  LEVEL} 
{NOISE  LEVEL} 
{DETECTION  THRESHOLD} 
{DIRECTIVITY  INDEX} 
{ABSORPTION  COEFFICIENT} 
{REVERBERATION  AREA} 
{AREA  SCATTERING  STRENGTH} 
{AREA  OF  UNCERTAINTY  IN  CONTACT  REPORTS} 

{FREQUENCY} 
{PULSE  LENGTH} 
{BEAMWIDTK} 

{UNCERTAINTY  IN  SONAR  EQN.  TERMS} 


- } 

{MAX  NUMBER  OF  SCOUTS} 
{UTILITY  CONSTANT} 
{NUMBER  OF  SONAR  PULSES  PER  SCAN  PERIOD} 
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<ENVIRONMENTAL  VARIABLES  AND  CONSTANTS} 


< - - - - 

const 

DEFBOTTOMTYP=1;. 

DEFSVTPR0FIL=1; 

DEFWIND  =10; 

DEFDEPTH  =200; 

- - - - - y 

{BOTTOM  TYPE> 
{SOUND  VELOCITY  PROFILE  TYPE> 
{WIND  INTENSITY  IN  KMOTS> 
{DEPTH  IN  M> 

var 

bottomtype  :byte; 

svtprof 1 letype  :byte; 
wind  : integer; 

depth  : integer; 

{TYPE  OF  BOTTOM  I  II  OR  III> 
{TYPE  OF  VERTICAL  SOUND  SPEED  PROFILE} 
{WIND  INTENSITY  IN  KNOTS} 
{DEPTH  IN  M} 

CSIMULATION  AND  CONTROL  FLOW  PARAMETERS} 
C - ^ - - - 

const 

DEFTOTALSIMTIME  =  15; 

DEFSCOUTSDELAY  =  6; 

- — - > 

{TOTAL  SIMULATION  TIME} 
{SCOUTS  DATUM  DELAY} 

var 

WannaContinue  : boo lean; 

WannaSeeOata  : boo lean; 
data  link  :text; 

dlinktotaltime  : integer; 

totalsimtime  :word; 

time  :word; 

scoutsdelay  :real; 

{EXIT  BOOLEAN  CONDITION} 
{GRAPHIC  DISPLAY  BOOLEAN  CONDI TON} 
{LOCATION  OF  THE  DATA  LINK  DEVICE  DRIVER} 
{TOTAL  SIMULATION  TIME  IN  DATA  LINK  DEVICE  DRIVER} 
{TIME  IN  MINUTES  FOR  WHICH  A  RENDER  IS  DESIRED} 
{DISCRETE  TIME  UNITS  COUNTER} 
{TIME  OFF  OF  THE  SEARCHERS} 

CPUBLIC  METHODS} 

C - — - 

- - - 

procedure  GetTargetData; 
procedure  GetScoutsData; 
procedure  GetEnvironmentalData; 
procedure  GetSinOata; 
procedure  AskOptions; 
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implementation 


uses  detecmanager; 


procedure  PromptOptions; 
var  i  :integer; 
begin 
clrscr; 

textcolor(Cyan); 

clrscr; 


writelnC' 
writelnC' 
writelnC' 
writelnC'  Please,  selec 
writeln; 

for  i:=1  to  16  do  writeln; 

text CO I or (Red); 

writelnC' 

writelnC' 

writelnC' 

writelnC' 

writelnC' 

writelnC' 

writelnC' 

writeln; 

write  C' 

end; 


Please,  select  one  of  the  following  options:'); 


0.  QUIT  PROGRAM'); 


1.  UPDATE  TARGET  DATA'); 

2.  UPDATE  SEARCHER  DATA'); 

3.  UPDATE  ENVIRONMENTAL  DATA'  ); 

4.  SET  TIME  VARIABLES'); 

5.  START  SIMULATION  WITH  CURRENT  DATA'); 

6.  DISPLAY  INITIAL  TARGET  POSITIONS'  ); 

Enter  your  option:  '); 


procedure  AskOptions; 

var 

answer  :char; 

origmode  : integer; 

begin 

answer:='9'; 

or i gmode : =LastMode; 

T  ext Mode ( C80+Font8x8) ; 

while  CCanswero'5')  and  Canswer<>'6'))  andCanswero'O')  do 
begin 

PromptOptions; 
answer :=readkey; 
case  answer  of 

'O':  WannaContinue:= FALSE; 

'1':  GetTargetData; 

'2':  GetScoutsData; 

'3':  GetEnvironmentalData; 

'4':  GetSimData; 

/<:/  •  . 

j  •  , 

'6':  WannaSeeData:=TRUE; 

end; 

end; 

T  extModeC  ori gmode) ; 

end; 
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procedure  GetTargetReaction; 

var 

count  :byte; 

begin 

textcolor(Red); 
clrscr; 
writeln(' 
writeln(' 
wn'telnC' 

writelnC'  Please  input  the  following  data  for  500  ^ 

for  count:=1  to  12  do  writeln; 

writeln;writeln; 

writeln('  TARGET  REACTION  PARAMETERS  /)- 

writeln;  ' 

writeln; 
writeln; 

writelnC'  Please  input  the  estimated  REACTION  RANGE  in  nautical  miles,  ')- 

writeln;  ' 

writelnC'  The  current  value  is  :  ',reactrange/1 0:3:1) 

write  ('  Enter  your  value  ;  /)- 

read I n( react range ) ; 

react  range := react  range* 1 0; 

writeln; 

writelnl'  Please  input  the  estimated  REACTION  TACK  in  nautical  degrees,  '); 

writelnC'  (Evasion  course  relative  to  the  bearing  to  the  closest  scout  ')' 

writeln;  ' 

writeln('  The  current  value  is  :  ',reactiontack:3); 

write  ('  Enter  your  value  .  /), 

read I n( react iontack);  ' 

end; 


procedure  GetDatumError; 


count  tbyte; 
begin 

text CO I or (Red); 
clrscr; 
writelnC' 
writelnC' 
writelnC' 

writelnC'  Please  input  the  following  data  for  500  targets:'); 
for  count:=1  to  12  do  writeln; 
writeln; writeln; 

writelnC'  DATUM  OMNIDIRECTIONAL  ERROR  (integer) 
writelnC'  (Standard  Deviation  in  Tgts  Position,  in  nm 


writelnC' 
writeln; 
writelnC' 
write  (' 


Current  value  is 
Enter  your  value 


readlnCsigmarinit); 

s i gmari n i t : =s i gmari ni t*1 0; 

writeln; 

writeln; 

writeln; 


',sigmarini 


t/10:3:0); 


writeln; 

end; 


101 


procedure  GetMotionModel ; 

var 

ans  :byte; 
count  :byte; 
begin 

textcolor(Red); 


clrscr; 

writelnC' 

writelnC' 

Please  input  the  following  data  for  500  targets:'); 

for  count :=1 

to  12  do  writeln; 

writeln;writeln; 

writelnC' 

writeln; 

TARGET  MOTION  MODEL 

writelnC' 

Cl)  OMNI  :  Omnidirectional  fleeing 

writelnC' 

C2)  FAN  :  Directional  movement'); 

writeln; 

writelnC' 

Current  motion  model  is 

:  ' ^motionmodel); 

write  C' 
readlnCans); 

Enter  your  model  Cl  OMNI,  2  FAN) 

:  '); 

case  ans  of 
1:  begin 

motionmodel  :=  'OMNI'; 
por tmost course :=0; 
courserange:=360; 
end; 

2:  motionmodel  :=  'FAN'; 
end; 
end; 


TARGET  DATA 


procedure  GetTargetCourses; 

var 

count  :byte; 
begin 

textcolor(Red); 
clrscr; 

wr i  t e I n( '  '  y" 

writelnC' 

writeLn('  „  „  _ 

writelh('  Please  input  the  following  data  for  500  targets:'); 

for  count :=1  to  12  do  writeln; 

writeln;writeln; 

writelnC'  TARGET  COURSE  LIMITS 
writeln; 
writeln; 
writeln; 
writelnC' 
writeln; 
writelnC' 
write  C' 


IP); 

IP); 


Please  input  the  courses  port  limit 


'); 


The  current  value  is 
Enter  your  value 
readlnCportmostcourse); 
writeln; 

writelnC'  Please  input  the  course  range 

writeln; 

writelnC'  The  current  value  is 

write  C'  Enter  your  value 

read I nC  courser ange ) ; 

writeln; 

end; 


' , por tmostcourse : 3 ) ; 
'); 


',courserange:3); 
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TARGET  DATA 


procedure  GetTargetSpeeds; 
var 

count  :byte; 
begin 

textcolor(Red); 
clrscr; 

writelnC' 

writelnC'  _ 

writeLnC'  Please  input  the  following  data  for  500  targets:'); 
for  count :=:1  to  12  do  writeln; 
writeln;writeln; 

writeln('  TARGET  SPEED  LIMITS 
writeln; 
writeln; 
writeln; 

writelnC'  Please  input  the  Minimum  estimated  Target  speeed 
writeln; 

writelnC'  The  current  value  is  :  ' ,roundCminsubspeed*1.2):3) 

write  ('  Enter  your  value  •  /)• 

read I n( mi nsubspeed ) ; 

minsubspeed:^  round(0.8333*mi nsubspeed );Cconver ting  to  discrete  speed  units> 
writeln; 

writelnC'  Please  input  the  Maximum  estimated  Target  speed 
writeln; 

HritelnC'  The  current  value  is  :  ',round<inaxsubspeecl*1 .2):3) 

write  ('  Enter  your  value  :  /); 

readln(maxsubspeed); 

maxsubspeed:=  round(0.8333*inaxsubspeed); {converting  to  discrete  speed  units> 

writeln; 

end; 


procedure  GetTargetData; 
var 

answering  :  char; 
count  :  byte; 

begin 

answer ing:='P'; 

while  answering  <>  '0'  do  begin 
textcolor(Red); 
clrscr; 

Default  values  of  the  parameters  usi^  to  generate  and  u^ate  'T; 
500  simulated  tracks:  '); 
for  count :=1  to  8  do  writeln; 
writeln; 

(0)  Back  to  Main 


writelnC' 

writelnC' 

writeln(' 

writeln(' 

writelnC' 


writelnC' 
writeln; 

write  C'  CD  DATUM  ERROR  IN  NM  : 

writeln;writeln; 

C2)  TARGET  MOTION  MODEL  ; 

Cspeed  &  course  uniformly  distributed)'); 


writelnC 

writelnC' 

writeln; 

writelnC' 

writeln; 

writelnC' 

writelnC' 

writeln; 

writelnC' 


C3)  COURSE  LIMITS  CFOR  FAN  MOTION  MODEL  ONLY) 

PORT  LIMIT  CO-360)  : 

RANGE  CNot  starbord  limit)  : 

C4)  SPEED'); 


D; 

',sigmarinit/10:3:0); 

',motionmodel); 


D; 

'  ,portmostcourse:3); 
',courserange:3); 
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wn'teln; 

writeln('  MINIMUM  ESTIMATED  TARGET  SPEED  IN  KTS 

writelnC'  MAXIMUM  ESTIMATED  TARGET  SPEED  IN  KTS 

wn’teln; 

wn'telnC^  (5)  TARGET  REACTION  MODEL 


wn’teln; 

wn’telnC'  REACTION  RANGE 

wn’teln('  REACTION  TACK 

wn’teln; 

wn’teln('  Enter  the  number  of  the  parameters  you  want 

write  to  go  back  to  main.  Return  is  not  necessary 

answer i ng:=readkey; 
case  answering  of 


'V  :  GetDatumError; 

*2*  :  GetMotionModel; 

'3'  ;  GetTargetCourses; 

'4'  :  GetTarget Speeds; 

'5'  :  GetTargetReaction; 

end; 

end;<whi le> 

end ; {GetTargetData} 


:  • , round(minsubspeed*1 .2) :3) 
:  • , round(maxsubspeed*l .2) :3) 

'); 

;  ',reactrange/10:3:0); 

:  react iontack:3); 

to  change  or  0  '); 

•  * 

-  it 


procedure  DisplayDataLinkInfo; 

var 

k  :byte; 

begin 

closer;  ^  iiiilllifp;  jjiilpi 

write  V  Connected  to  Data  Link,  #  of  scouts  present  :  '); 
textcolor(140); 

writelnC'  •>— >  number  of  scouts); 

textcolor(red); 

write  ('  Time  searching  so  far  :  •); 

textcolor(140); 

writelnC'  — >  ^dlinktotal timers,'  min'); 

textcolor(red); 

writeln; 

writelnC'  Receiving  the  following  initial  conditions:'); 

writeln; 

writeln; 

writelnC'  SCOUT  #  XPOS  YPOS  COURSE  SPEED  SONAR_^ON  SCALE  '); 

writeln; 

for  k:=1  to  nuntoerof scouts  do  begin 
with  myscoutCk]  do  begin 

writelnC'  ',k:8  x:8  ,  y:8,  course:8,  speed:8:0,SONAR_ON:8,maxscale*185:8); 

end; 
end; 

writeln; 

writeln; 

writelnC'  Hit  any  key  to  continue...'); 

readkey; 

end; 
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SCOUTS  DATA 


procedure  ChangeSonarDefaultsCscout.param.-byte) : 

var 

f  :real; 

k  :byte; 

begin 
cLrscr; 
writelnC' 
writelnC' 
writelnC' 

writelnC'  Scouts  Sonar  Parameters:'); 

writeln; 
writeln; 

writeln<'  |  /). 

writelnC'  |  (1)  (2)  (3)  (4)  (5)  (6)  |  '); 

writeln('  |SCOUT  SL  DT  DI  Frq.  PulseL.  Beam  W.  |  '); 

writeln('  |  #  (dB)  (dB)  (dB)  (Khz)  (msec)  (deg)  I  '); 

writelnC'  [  - -  | 

for  k:=1  to  nunberof scouts  do  begin 
with  myscoutCk]  do  begin 

writelnC'  |'fk:5  ',SL:5/  '  ,DT:5,'  ',01:5/  ',frequency:8:1, ' 

,pulselength*1000:8:0,'  ',beamwidth:8:1, '  |'); 

end; 

end; 

writelnC'  .  ■  n.y 

writeln; 

writelnC'  You  are  going  to  produce  modifications  to  scout  number:  '  ,scout); 
with  myscout [scout]  do  begin 
if  param  =  1  then  begin 
writeln; 

writeC'  Please  input  new  SOURCE  LEVEL  :  /)- 

readlnCSL); 
end; 

if  param  =  2  then  begin 
writeln; 

writeC'  Please  input  new  DETECTION  THRESHOLD  •  /)• 

readlnCDT);  *  ' 

end; 

if  param  =  3  then  begin 
writeln; 

writeC'  Please  input  new  DIRECTIVITY  INDEX  :  o- 

readlnCDI); 
end; 

if  param  =  4  then  begin 
writeln; 

writeC'  Please  input  new  FREQUENCYC updates  aborption  coeff)  :  '); 

read I nC frequency);  ’ 

abscoef f :=Absorpt i onCoef f i c i ent  C  f  requency ) ; 
end; 

if  param  =  5  then  begin 
writeln; 

writeC'  Please  input  new  PULSE  LENGTHCin  ms)  :  /); 

read I nCpulse length); 
pu I se I eng t h : =pu I se I eng th/ 1 000; 
end; 

if  param  =  6  then  begin 
writeln; 

writeC'  Please  input  new  BEAM  WIDTHCin  degrees)  :  #). 

readlnCbeamwidth); 
end; 
end; 
end: 
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procedure  Di spl ayandChangeSonarDef aul ts ; 
var 

k,iteni  :byte; 
scoutmod  :byte; 
option  :char; 
begin 

while  scoutmod  <>  0  do  begin 
clrscr; 


SCOUTS  DATA 


CD 

C2) 

(3) 

C4) 

C5) 

C6) 

1 ') 

SL 

DT 

DI 

Frq. 

PulseL. 

Beam  W. 

1 ') 

CdB) 

CdB) 

(dB) 

CKhz) 

Cmsec) 

Cdeg) 

1 ') 

writelnC' 
writeln(' 
writeln(' 

writelnC'  Scouts  Sonar  Parameters:'); 

writeln; 

writeln; 

writelnC'  . 

writelnC'  CD  C2)  C3)  C4)  C5)  C6) 

writelnC'  SCOUT  SL  DT  DI  Frq.  PulseL.  Beam  W. 

writelnC'  #  CdB)  CdB)  CdB)  CKhz)  Ctnsec)  Cdeg) 

writelnC'  [i— 

for  k:=1  to  numberof scouts  do  begin 
with  myscoutCk]  do  begin 

writelnC'  |',k;5  ',SL:5,'  '  ,DT:5,'  ',DI:5,'  frequency: 

,pulselength*1000:8:0,'  ',beamwidth:8:1, '  |'); 

end; 

end; 

writelnC'  ^ 

writeln; 

writelnC'  Please  input  the  number  of  the  searcher  you  want  to  modify 
write  C'  or  0  to  go  back  to  Scouts  Menu  :  '); 
read I nC scoutmod) ; 
writeln; 

if  scoutmod  <>  0  then  begin 

if  scoutmod  >  numberof scouts  then  begin 

writelnC'Please  input  scout  #  smaller  than  ',numberofscouts); 
writelnC'Hit  any  key  to  continue- .'); 
readkey; 
end  else  begin 

write  C'  Please  input  the  column  of  the  item  you  want  to  modify 

read I nCi tern); 

writeln; 

ChangeSonarDef aul tsCscoutmod, i tern); 
end; 

end;C:if> 
end;Cwhi le> 
end; 


'  ,DT:5,'  ',DI:5,'  ',frequency:8:1 

,beamwidth:8:1,'  |'); 


Please  input  the  number  of  the  searcher  you  want  to  modify  '); 
or  0  to  go  back  to  Scouts  Menu  :  '); 


procedure  Di spl ayandChangeErrorDefaul ts ; 
var 

thisanswer  : integer; 
begin 
clrscr; 
writelnC' 
writelnC' 
writelnC' 

writelnC'  ERROR 

writeln; 

writeln; 

writeln; 

writelnC'  Total  estimated  error  in  t 
writelnC'  calculation  of  the  sonar  e 


SCOUTS  DATA 


ERROR  VARIABLE  DEFAULT  VALUE 


Total  estimated  error  in  the  calculation  of  the  terms  in  the 
calculation  of  the  sonar  equation  and  in  the  uncertainties  in 
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wn’telnC'  the  detection  process  in  dB  re  A/Pascal, 
writeln; 

writelnC'  “NOISESIGMA"  current  value  is  :  noises igma:6  /  dB  re  /iPascal'); 

writeln; 

write  i*  Please  input  the  new  valaue  : 

read I n( noises igma); 

end;  {procedure} 


Please  select  one  of  the  following  options: 


procedure  GetScoutsData; 

var 

response  :char; 
begin 

response: ='9'; 

while  response  <>  'O'do  begin 
clrscr; 

writeln(' 
writeln(' 
writelnC' 
writelnC' 
writeln; 
writeln; 
writeln; 
writeln; 
writeln; 
writelnC^ 
writelnC' 
writeln(' 
writelnC' 
writeln; 
write 

response: =readkey; 
case  response  of 
M':  DisplayDataLinkInfo; 

*2*  :  DisplayandChangeSonarDefaults; 
'3' :  DisplayandChangeErrorDefaults; 
^0':  ; 
end; 
end; 
end; 


0. 

1. 

2. 

3. 


BACK  TO  KAIN  '); 

DISPLAY  INFO  RECEIVED  FROM  SCOUTS 
DISPLAY  AND  OR  CHANGE  DEFAULT  SONAR  PARAMETERS'); 
DISPLAY  AND  OR  CHANGE  DEFAULT  ERROR  VALUES  '); 

Enter  your  option:  '); 


procedure  GetSimData; 

var 

answer :char; 

begin 

clrscr; 

writelnC'  . .  '); 

writelnC'  The  simulation  will  ^  carrii^^  out  from  time  6  to  Tot  a  fs^ 

writelnC'  in  updating  steps  of  5  minutes.  Totalsimtime  should  be  smaller  '); 

writelnC'  than  the  last  updating  time  to  the  data-link  device  file.  ')- 

writelnC'  ,y 

writeln; 
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writeln; 

writeln; 

writelnC'  Total  Simulation  Time  (TotalSimtime,  a  multiple  of  5  min)'); 
write  ('  No  info  from  scouts  available  after  :  '); 

textcolor(140); 

wr  i  tel  n(d  I  i  nktota  1 1  i  me*5  /  min'); 

text CO I or (red); 

writeln; 

writeln('  Current  value  is  :  ',totalsimtinie*5,'  min'); 

write  ('  Enter  your  value  :  '); 

readln(totalsimtime); 

totalsimtime:=totalsimtime  div  5; 

if  totalsimtime  >  dli nktota It ime  then  begin 

writelnC'  Simulation  time  should  be  smaller  than  ',dlinktotaltime*5); 
writeln; wr i te I n; wr i te I n; wr i te I n; wr i te I n; 
writelnC'  Hit  any  key  to  continue...'); 
readkey; 
end; 

writeln; 

writeln; 

writelnC'  Delay  Time  of  scouts  after  DATUM  time,(DATUM  ”age»  in  min)  '); 

writelnC'  Current  value  is  :  ',scoutsdelay*5:4:0); 

write  ('  Enter  your  value  :  '); 

read ln( scoutsde I ay) ; 

scoutsde I ay: =scoutsde I ay/5 ; 

end; 


function  Profileinwords: string; 

var 

aistring; 

begin 

case  svtprof iletype  of 

1:  a:='POSITIVE  OR  ISOSPEED'; 
2:  a:='NEGATIVE'; 
end; 


Prof i leinwords:=a; 
end; 


function  Bottominwords: string; 
begin 

case  bottomtype  of 

1:  Bottominwords:='SANDY,  FLAT'; 

2:  Bottominwords:^' INTERMEDIATE  BETWEEN  VERY  ROUGH  AND  FLAT'; 
3:  Bottominwords:='ROCICY,  VERY  ROUGH'; 
end; 
end; 


procedure  Getwindspeed; 

begin 

clrscr; 


writelnC' 

writelnC' 

writelnC' 

writelnC' 


This  is  the  current  environmental  data: 


writeln; 


writeln; 


writeln; 


'); 

'); 


writelnC'  (1)  WINDSPEED  :  ',wind:4.'  Knots'  ); 
writelnC'  (2)  SOUND  SPEED  PROOFILE  :  '.Profileinwords  ); 
writelnC'  C3)  BOTTOM  TYPE  :  '.Bottominwords  ); 
writelnC'  C4)  DEPTH  :  ',depth:4.'  m'  ); 
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) 


write In; 

write  ('  Please  input  wind  speed  (knots):  • 

read I n( wind); 

end; 


_  ENyiRONME^TAf^^^^D^^^ 

This  is  the  current  environmental  data: 


procedure  GetSoundProfile; 

var 

OUTTAHERE  iboolean; 
begin 

OUTTAHERE:=FALSE; 
while  not  OUTTAHERE  do  begin 
clrscr; 
writelnC' 
writelnC' 
writelnC' 
writelnC' 
write In; 
writeln; 
writelnC' 
writeln(' 
writeln(' 
writeln(' 
writeln; 
writeln; 
writeln(^ 
writelnC' 
writeln; 
writeln(' 
writelnC' 
writeln; 

write  V  Please  enter  your  option  here  :  ' 

readlnCsvtprof i letype); 

if  (svtprofiletype=  1)  or  (svtprofi letype  =  2)  then  begin 
OUTTAHERE :=TRUE; 
end; 
end; 


II'); 
III'); 
); 
); 


(1)  UINDSPEED 

(2)  SOUND  SPEED  PROOFILE 

(3)  BOTTOM  TYPE 

(4)  DEPTH 


',wind:4,'  Knots' 
',Prof i leinwords 
',Bottominwords 
depth :4,'  m' 


); 

>; 

); 

); 


You  can  choose  the  sound  speed  profile  according  to  the  following'); 
code:'); 


(1)  POSITIVE  OR  ISOSPEED 

(2)  NEGATIVE 


>; 


procedure  GetBottomtype; 

var 


OUTTAHERE  :boolean; 
begin 

OUTTAHERE:=FALSE; 

while  not  OUTTAHERE  do  begin 

clrscr; 

writelnC'  “ "  ^ENmONMEN^^^^^  DATA 

writelnC'  This  is  the  current  bottom  &  depth  environmentaf  ciata: 


') 

') 

') 


writeln; 

writeln; 

writeln(' 

writelnC' 

writelnC' 

writelnC' 

writeln; 

writeln; 

writelnC' 

writeln; 

writelnC' 

writelnC' 

writelnC' 


(1)  WINDSPEED  :  ',wind:4,'  Knots'  ); 
C2)  SOUND  SPEED  PROOFILE  :  Prof i lei nwords  ); 
C3)  BOTTOM  TYPE  :  '  ,Bottominwords  ); 
(4)  DEPTH  :  ',depth:4,'  m'  ); 


You  can  input  bottom  type  according  to  the  following  code'); 

(1)  SANDY, FLAT  BOTTOM, NEGLECTABLE  ROUGHNEESS  '); 

(2)  INTERMEDIATE  BETWEEN  FLAT  AND  VERY  ROUGH  '); 

C3)  ROCKY,  VERY  ROUGH  /)' 
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wn’teln; 

write  Please  enter  your  option  here  :  ' 

read  I  n(  bot  t omtype) ; 

if  (bottomtype  =  1)  or  (bottomtype  =  2)  or  (bottomtype  =  3)  then  begin 
OUTTAHERE:=TRUE; 
end; 
end; 
end; 


); 


read  ln(  depth); 
end; 


procedure  Getdepth; 

begin 

clrscr; 

writelnC' 

writelnC' 

writelnC' 

writeln; 

writeln; 

writeln; 

OlliiiiiiiOyiiS  environmental  data 

This  is  the  current  environmental  data: 

writelnC' 

(1)  WINDSPEED  :  ',wind:4/  Knots' 

); 

writelnC' 

C2)  SOUND  SPEED  PROOFILE  :  ',Prof ileinwords 

); 

writelnC' 

C3)  BOTTOM  TYPE  :  ',Bottominwords 

); 

writelnC' 

writeln; 

C4)  DEPTH  :  ',depth:4/  m' 

); 

write  C' 

Please  input  depth  Cmeters):  ' 

); 

•:r*  : 

_ 

!il  ENV 

^IRONMENTAL 

data 

procedure  GetEnvironmentalData; 

var 

hisanswer  :  char; 

begin 

while  hisanswer  <>  ^0'  do  begin 
clrscr; 
writelnC' 
writeln(' 
writelnC' 
writelnC' 
writeln; 
writeln; 
writeln; 
writeln(' 
writeln(' 
writelnC' 
writelnC' 
writelnC' 
writeln; 
writelnC' 
write  C^ 


This  is  the  current  environmental  data: 


CO)  BACK  TO  MAIN 

Cl)  UINDSPEED  :  ',wind:4,'  Knots' 

C2)  SOUND  SPEED  PROOFILE  :  Prof i lei nwords 

C3)  BOTTOM  TYPE 
C4)  DEPTH 


'fBottominwords 
', depth :4,'  m' 


Please  input  the  numer  of  the 
parameters  you  want  to  change 
h i sanswer : =readkey; 
case  hisanswer  of 


'1' 

:  Getwi ndspeed; 

'2' 

:  GetSoundProfi le; 

'3' 

:  Get Bottomtype; 

'4' 

:  GetDepth; 

end; 

end; 

end; 

'); 

l'>; 


); 

); 

); 

); 

'); 

'); 
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begin 


<TARGET  INIT> 
sigmarimt 
portmost course 
courserange 
mi nsubspeed 
maxsubspeed 
react range 
react iontack 
mot ionmode I 


:=DEFSIGMARINIT; 

:=DEFPORTMOSTCOURSE; 

:=:DEFCOURSERANGE; 

:=DEFNINSUBSPEED; 

:=DEFMAXSUBSPEED; 

:=DEFREACTRANAGE; 

:=DEFREACTIOMTACIC; 

:=0MNI; 


CSEARCHER  INIT> 

no i ses i gma : =DE FNOI SES I G ; 

<SIMULATION  INIT> 

totalsimtime  :-DEFTOTALSIMTIME; 
WannaCont i nue:=TRUE; 
scoutsdelay  :=DEFSCOUTSDELAY; 


CENVIRONMENT  INIT> 
bottomtype  :=DEFBOTTOMTYP; 

svtprofiletype  :=DEFSVTPROFIL; 
Mind  :=DEFWIND; 

depth  :=DEFDEPTH; 

end. 


Ill 
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APPENDIX  C.  USER  INTERFACE  SEQUENCES 

1.0  GENERAL  MENU 


Please,  select  one  of  the  following  options: 


0.  QUIT  PROGRAM 

1.  UPDATE  TARGET  DATA 

2.  UPDATE  SEARCHER  DATA 

3.  UPDATE  ENVIRONMENTAL  DATA' 

4.  SET  TIME  VARIABLES 

5.  START  SIMULATION  WITH  CURRENT  DATA 

6.  DISPLAY  INITIAL  TARGET  POSITIONS' 

Enter  your  option; 


2.1  TARGET  PARAMETERS  MAIN  MENU 


Default  values  of  the  parameters  used  to  generate  and  update 
500  simulated  tracks: 


(0)  Back  to  Main 

(1)  DATUM  ERROR  IN  NM  : 

(2)  TARGET  MOTION  MODEL  : 

(speed  &  course  uniformly  distributed) 

(3)  COURSE  LIMITS  (FOR  FAN  MOTION  MODEL  OHLY) 

PORT  LIMIT  (0-360)  : 

RANGE  (Not  starbord  limit)  : 

(4)  SPEED 

MINIMUM  ESTIMATED  TARGET  SPEED  IN  KTS  : 

MAXIMUM  ESTIMATED  TARGET  SPEED  IN  KTS  : 

(5)  TARGET  REACTION  MODEL 

REACTION  RANGE  : 

REACTION  TACK  : 


Enter  the  number  of  the  parameters  you  want  to  change  or  0 
to  go  back  to  main.  Return  is  not  necessary  :  ' 


2.2  ERROR  INPUT 


Please  input  the  following  data  for  500  targets: 

DATUM  OMNIDIRECTIONAL  ERROR  (integer) 

(Standard  Deviation  in  Tgts  Position,  in  nm 

Current  value  is  : 

Enter  your  value  s 

2.3.  MOTION  MODEL  INPUT 


Please  input  the  following  data  for  500  targets: 


TARGET  MOTION  MODEL 

(1)  OMNI  :  Omnidirectional  fleeing 
<2)  FAN  :  Directional  movement 

Current  motion  model  is 
Enter  your  model  (1  OMNI,  2  FAN) 


2.4.  TARGET  COURSE  LIMITS  INPUT 


TARGirDATA 

sHiHHliin;! 

IS 

_ 

liii 

.si. 

Please  input  the  following  data  for  500  targets: 

TARGET  COURSE  LIMITS 

Please  input  the  courses  port  limit 

The  current  value  is  : 

Enter  your  value  : 

Please  input  the  course  range 

The  current  value  is  : 

Enter  your  value  : 
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2.5.  TARGET  SPEED  ESTIMATE  INPUTS 


2.6.  TARGET  REACTION  PARAMETERS  INPUTS 


Please  input  the  following  data  forToO  targets 
TARGET  REACTION  PARAMETERS 

Please  input  the  estimated  REACTION  RANGE  in  nautical  miles. 

The  current  value  is  : 

Enter  your  value  : 


Please  input  the  estimated  REACTION  TACK  in  nautical  degrees, 
(Evasion  course  relative  to  the  bearing  to  the  closest  scout 

The  current  value  is  : 

Enter  your  value  : 
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3. 


SCOUTS  MAIN  MENU 


SCOUTS  DATA 


Please  select  one  of  the  following  options: 


0.  BACK  TO  MAIN 

1-  DISPLAY  INFO  RECEIVED  FROM  SCOUTS 

2.  DISPLAY  AND  OR  CHANGE  DEFAULT  SONAR  PARAMETERS 

3.  DISPLAY  AND  OR  CHANGE  DEFAULT  ERROR  VALUES 

Enter  your  option: 


3.2.  SCOUTS  DATA  LINK  DISPLAY 
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3. 


SONAR  PARAMETERS  DISPLAY  AND  INPUT 


SCOUT 

# 

(1) 

SL 

(dB) 

(2) 

DT 

(dB) 

(3) 

DI 

(dB) 

(4) 

Frq. 

(Khz) 

(5) 

PulseL. 

(msec) 

(6) 

Beam  W. 
(deg) 

1 

1 

1 

You  are  going  to  produce  modifications  to  scout  number: 
(eventual) 

Please  input  new  SOURCE  LEVEL  : 
Please  input  new  DETECTION  THRESHOLD  : 
Please  input  new  DIRECTIVITY  INDEX  ; 
Please  input  new  FREQUENCY(updates  aborption  coeff)  : 
Please  input  new  PULSE  LENGTH(in  ms)  : 
Please  input  new  BEAN  UIDTH(in  degrees)  : 


4.  SONAR  EQUATION  ERROR  INPUT 


“  error"’ vMable"  default""\wil’ue 


Total  estimated  error  in  the  calculation  of  the  terms  in  the 
calculation  of  the  sonar  equation  and  in  the  uncertainties  in 
the  detection  process  in  dB  re  /ypascal. 


"NOISESIGHA"  current  value  is 
Please  input  the  new  valaue 


117 


4.1 


ENVIRONMENTAL  MAIN  MENU 


ENVIRONMENTAL 


DATA 


This  is  the  current  environmental  data: 


(0)  BACK  TO  MAIN 

(1)  UINDSPEED  : 

(2)  SOUND  SPEED  PROOFILE  : 

(3)  BOTTOM  TYPE  : 


4.2.  WIND  SPEED  INPUT 


This  is  the  current  environmental  data: 


(1)  UINDSPEED  : 

(2)  SOUND  SPEED  PROOFILE  : 

(3)  BOTTOM  TYPE  : 

(4)  DEPTH  ; 

Please  input  wind  speed  (knots): 
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.3. 


SOUND  SPEED  PROFILE  TYPE  INPUT 


This  is  the  current  environmental  ^~ta: 


(1)  WINDSPEED 

(2)  SOUND  SPEED  PROOFILE 

(3)  BOTTOM  TYPE 

(4)  DEPTH 


You  can  choose  the  sound  speed  profile  according  to  the  following 
code: 


(1)  POSITIVE  OR  ISOSPEED 

(2)  NEGATIVE 

Please  enter  your  option  here  : 


.4.  BOTTOM  TYPE  INPUT 


ENVIRONMENTAL  DATA  "  "  gisi5 

This  is  the  current  bottom  &  depth  environmental  data: 


(1)  WINDSPEED 

(2)  SOUND  SPEED  PROOFILE 
{3)  BOTTOM  TYPE 

(4)  DEPTH 


You  can  input  bottom  type  according  to  the  following  code 

(1)  SAMDy,FLAT  BOTTOM, NEGLECTABLE  ROUGHNEESS 

(2)  INTERMEDIATE  BETWEEN  FLAT  AND  VERY  ROUGH 

(3)  ROCKY,  VERY  ROUGH 


Please  enter  your  option  here  : 


5.0  TIME  VARIABLES  INPUT 


The  simulation  will  be  carried  out  from  time  0  to  TotalSimtime 
in  updating  steps  of  5  minutes.  Totalsimtime  should  be  smaller 
than  the  last  updating  time  to  the  data- I ink  device  file. 


Total  Simulation  Time  (Totalsimtime,  a  multiple  of  5  min) 
No  info  from  scouts  available  after  : 

Current  value  is  : 

Enter  your  value  : 


Hit  any  key  to  continue... 


Delay  Time  of  scouts  after  DATUM  time,<DATUM  "age"  in  min) 
Current  value  is  : 

Enter  your  value  : 
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APPENDIX  D.  TARGETMANAGER  UNIT  CODE 


unit  TargtManager; 

interface 

uses 

GRAPH,  CRT; 
van 

function  normaUvar  u:double):real; 
procedure  GenTargets; 

procedure  TargetReact i on( t i , ss: i nteger ;di stst : rea I ) ; 
procedure  UpdateTargetPositionCti : integer); 
procedure  Upda teWeight sUith Bayes ; 


implementation 

uses  GraphManager.InterManager; 


function  normaUvar  u: double): real; 

^converts  uniform  random  numbers  to  standard  normal> 

<by  Alan  Uashburn> 
const 

a=2.5155l7;b=0.802853;c=0.0l0328;d=1.432788;e=0.189269;f=0-001308; 

var 

sign: integer; 
x,y:real; 
begin 

sign:=-1; 

if  u>0.5  then  begin 
sign:=1; 
u:=1.0-u; 
end; 

y:=sqrt(-ln(u*u)); 

X  :=(  a+y*(bH-c*y)  )/<  1+y*(d+y*<e+f  *y)  )  )  ; 
normal :=sign*(y-x); 
end; {normal} 


function  bearing(xl.yl,x2,y2:integer) : integer; 

CFROM  1  TO  2> 
var 

x,y  : integer; 


begin 

x:=x2-x1; 

y:=y2-y1; 


if  x>0  then  begin 

if  y  >  0  then  bearing 
if  y  =  0  then  bearing 
if  y  <  0  then  bearing 
end; 

if  x<0  then  begin 

if  y  <  0  then  bearing 
if  y  =  0  then  bearing 


round(180  -  arctan(x/y)*180/pi); 
90; 

-round(arctan(x/y)*180/pi); 


round(360  -  arc tan(x/y)* 180/pi ); 
270; 
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if  y  >  0  then  bearing  :=  roundClSO  -  arctan(x/y)*180/pi); 
end; 

if  x=0  then  begin 

if  y  >=0  then  bearing  :=  180; 
if  y  <  0  then  bearing  :=  0; 
end; 
end; 


procedure  GenTargets; 

var 

r  : double; 

i  iword; 

datumerror  : integer; 
range  :real; 

begin 

Randomize; 

for  i:=1  to  NUMBEROFTGTS  do  begin 
r:=random; 

ntytarget[i].x  :=  CENTER  +  round(normal(r)  *  sigmarinit); 
r:=random; 

mytarget[i].y  :=  CENTER  +  round<normal(r)  *  sigmarinit); 

mytargetCi]  .datumx:=mytarget  Ci)  .x; 

mytargetCi]  .datuniy:=inytarget  Ci]  -y; 

mytarget  ti]  -  course  :=portmos  tcourse+randomC  courserange  ); 

mytargetEi]  -speed:=minsubspeecH-random*<maxsubspeed-minsubspeed); 

mytarget[i].x:=mytargetti].x  +  roundCmytargetCi)  .speed  ♦  scoutsdelay  * 

si nCmy target  Ei ] . course*pi/l 80) ); 

mytargetCi]  .y:=inytargetEi]  .y  -  roundCmytarget  Ei)  .spe^  *  scoutsdelay  * 

cosCmytargetCi] -course*pi/180)); 
mytargetEi]  .cel lx  :=Cmytarget [i] .x  div  GRIDSTEP)+1; 
mytargetEi] -cel ly  :=:(mytarget Ei3 -y  div  GRIDSTEP)+1; 
mytargetEi]  .prob:s0.002; 
mytarget  Ei]  .threatdis:=BIGNUMBER; 
end; 
end; 


procedure  UpdateTargetPosition(ti : integer) ; 

var 

r  : double; 

begin 

mytarget Et i]. x:=mytarget Et i]  .X  +  roundCmytargetEti] .speed  * 

sin(mytargetEti] -course*pi/180))  ; 
mytarget Eti].y:=mytargetEti].y  -  roundCmytargetEti] .speed  * 

cos ( mytarget  E t i 3 . course*pi /1 80 ) ) ; 
mytarget  Eti  3.  cel  lx  :=(inytarget  Eti]  .x  div  GRIDSTEP)+1; 
mytarget Eti]. cel ly  :=(mytarget Eti] .y  div  GRIDSTEP)+1; 
end; 


procedure  TargetReact i on ( t i , ss : i nteger ; di st st : real ) ; 

var 

filehandle  :text; 

r  :double; 

hearts  : integer; 

begin 

if  distst  <  reactrange  then  begin 

bearts:=bearing(mytarget Eti] .x,mytarget Eti] .y,myscout Css] .x,myscout [ss] .y); 
if  hearts  >  360  then  hearts :=bearts-360; 
if  hearts  <  0  then  bea rts: =360- hearts; 
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if  sin((iiiyscout[ss3.course-bearts)*pi/180)  >  0  then  begin 
m/target[ti3.course:=bearts  -  reactiontack; 
end  else  begin 

mytargetiti] .course:=bearts  +  reactiontack; 

end; 

end  else  begin 

if  motiormodel  =  'OMNI'  then  begin 

mytargetiti] . course :=bearing(mytarget [ti]  .datumx,mytarget  [ti]  .datumy, 

mytarget [ti] .x^mytarget [ti]  .y); 

end  else  begin 

mytarget  [ti] . course :=portmostcourse+random( courserange) ; 
end; 
end; 

if  mytarget  [ti] .course>360  then  begin 

mytarget [ti] . course :=my target [ti] .course-360; 
end; 

if  mytargetiti] .course<0  then  begin 

mytarget [ti] . course :s360+mytarget[ti] .course; 
end; 

if  mytarget [ti].threatdis  >  distst  then  begin 
mytargetiti] .threatdis:=distst; 
mytargetiti] .evasioncourse:=mytarget [ti] .course; 
end; 

mytarget [ti] -course:=mytarget  [ti] .evasioncourse; 
if  ss=numberofscouts  then  mytarget It i] .threatdis:=B I GNUMBER; 
end ; {TargetReacti on} 


procedure  UpdateWeightsWithBayes; 
var 

i  : integer; 
begin 

for  i;=l  to  NUMBEROFTGTS  do  mytarget  [i 3  .prob:=inytargetIi]  .prob/totprob; 
end; 


begin 

end. 
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APPENDIX  E.  SCOUTMANAGER  UNIT 


unit  ScoutManager; 

interface 

uses 

GRAPH,  CRT; 

procedure  InitScouts; 
procedure  CLoseScoutsDataLink; 
procedure  ConnectWi thScoutsThrtOataLink; 
procedure  UpdatePr i orWi thScoutsInf o; 
procedure  ReadScoutsDataLink; 
procedure  ResetUndetectedMass; 
procedure  UpdateUndetectecMass; 

implementation 

uses  GraphManager .TargtManager . Intermanager »DetecManager; 


procedure  Ini ti a1 i zeScoutsVarsWi thOef aul ts ; 

var 


k  : integer ; 

begin 

for  k:=1  to  numberof scouts  do  begin 


end; 


myscout [k] 
myscout [k] 
myscout [k] 
myscout  Ck] 
myscout [k] 
myscout [k] 
myscout  Ek] 
myscout [k] 
myscout [k] 
myscout [k] 
myscout [k] 
end; 


.abscoeff  :=DEFabscoeff; 
.SL  :=DEFSL; 

-NL  :=DEFNL; 

-DT  :=DEFDT; 

-DI  :=DEFDI; 

-SS  :=DEFSCATTSTRG; 

■IS^HOT  :=FALSE; 

.  frequency  :=DEF FREQUENCY; 
.beamwidth  :=DEFBEAHWIDTH; 
.pulselength:=DEFPULSELENG; 
.  reverbarea:=DEFREVRBAREA; 


procedure  Determi  neNumberOf Scout sAndTi  me ; 

var 

1  : integer; 

begin 

read  Ln(datal i nk, numberof scouts, d I i nktota 1 1 ime) ; 
for  i:=1  to  numberofscouts  do  begin 
readln<datalink,myscout£i]  .id); 
end; 
end; 


procedure  InitScouts; 
begin 

assignCdatal ink, 'c:\lnk1 .dat'); 
reset(datalink); 

Determi  neNumberOf  Scout  sAndT  i  me; 
close(datalfnk); 

Initial! zeScoutsVa  rsWi thDef aul ts; 
end; 
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procedure  CloseScoutsDataLink; 
begin 

close(datalink); 

end: 


procedure  ConnectWi  t hScout sThruDat aLi  nk ; 
begin 

reset (data I ink); 

Determi  neNumberOf  Scout  sAndT  i  me; 

ReadScoutsDatal ink; 

end; 


procedure  ReadScoutsDataLi nk ; 
var 

sonarinfo  :byte; 
contact info: byte; 
k  :byte; 

begin 

for  k:=l  to  numberof scouts  do  begin 
read(datalink,myscout [k] .x); 
read(datalink); 
read(datalink,myscout  £k] .y); 
read(datal i nk^myscout [k] .course) ; 
read(datalink,myscout [k] -speed); 
readCdata link, sonar i nf o) ; 
readCdatal ink,myscout [k] .maxscale); 
readCdata link, contact i nf o) ; 
if  sonarinfo  =  1  then  begin 
myscout [k] - SONAR_ON : =TRUE ; 
end; 

if  contact info  =  1  then  begin 

if  not  myscout [k] . I S_H0T  then  begin 
myscout [k] . i ni t i a  I hot : =t i  me*5 ; 
end; 

myscout [k] . I S_HOT  : aTRUE ; 

soundC 1 500 ) ; de I ay ( 40 ) ; 
soundC 1 600 ) ;de I ay ( 10); 
sound( 800 ) ; de  L  ay ( 200 ) ; 
soundC 1 600 ) ; de I ay ( 60 ) ; 
nosound; 

readC  data link, myscout [k] . subx ) ; 
readCdatalink,myscout Ck3 -suby); 
readCdata I ink, myscout  Ck] -subcourse); 
readCdata I ink,myscout [k] .subspeed) ; 
readCdata I ink, myscout [k] .cred); 
end  else  begin 

myscout  Ck3 . 1 S_HOT  : =FALSE ; 

end; 

readlnCdatalink); 

end; 

end; 


procedure  TempUpdateNegativeInfo(s:byte) ; 
var 

i  : integer; 

cdp5  :real; 

d  : rea I ; 

begin 

for  i:=l  to  NUMBEROFTGTS  do  begin 


d:=sqrt(distancer(mytarget  [1]  .x^mytargetCi]  .VfinyscoutEs]  .x,m/scout[s]  -y)); 
if  d  <  myscout [s] .maxscate  then  begin 
cct)5 :  =Detec  t  i  onE  f  for  t  ( i ,  s ,  d ,  depth ) ; 
mytarget[i]-prob:=inytargetCi3.prob  *  (1-cdp5); 
end; 

if  s=numberof scouts  then  begin 
UpdateTargetPpsitionC i ); 
totprob:=:totprob  +  mytarget  Ci]  .prob; 
end; 
end; 
end: 


procedure  TempUpdatePositivelnfoCsrbyte) ; 

var 

i  : integer; 

cdp5  :reaL; 

C  :reaL; 

d  rreal; 

radweight:real; 
distoff  rreal; 

error  rreal; 

begin 

for  ir=1  to  NIMBEROFTGTS  do  begin 

dr=sqrt(distancer(inytarget  Ci]  .x,mytargetCi]  .y, myscout  ts] , x, myscout  Es]  -y)); 
cdp5  r  =Detect  i  onEf  fort  ( i ,  s ,  d ,  dept  h )  ; 

Cr=myscout  Cs] .cred; 

distoff r=sqrt(distancer(mytargetEi] .x,mytarget Ci] .y,myscout Cs] -subx, myscout Cs] .suby)) 
error r=  d  *  myscout Cs] .beamwidth  *  pi  /180; 
radweightr=exp(-sqr(distoff/error>/2); 
mytarget Ci]-probr=  mytarget Ci] .prob  *  <1-C) 

+mytargetCi] .prob  *  cdp5  *  radHeight*C; 

TargetReaction(i,s,d); 
if  s=numberof scouts  then  begin 
UpdateTargetPos i t i onC i ) ; 
totprobr=totprob  +  mytarget Ci] .prob; 
end; 
end; 
end; 


procedure  UpdatePr 1 orWi t hScouts Inf  o ; 
var 

scout  rinteger; 

begin 

totprobr =0-0000000000001 ; 

for  scout r=l  to  numberof scouts  do  begin 
if  odd(time)  then  begin 
UpdateNL(scout); 
end; 

if  not  myscout Cscout] -I S_HOT  then  begin 
T  empUpda teNegat i vel nf o(scout ) ; 
end  else  begin 

T  empUpdatePos i t i vel nf o( scout  > ; 
end; 
end; 
end; 


procedure  ResetUndetectedMass; 
begin 
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ivKletectedtnass:=1 ; 

end; 


procedure  UpdateUndetectedHass; 
begin 

undetectecknass  :=undetectedmass*totprob; 

end; 


begin 

end. 
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APPENDIX  F*  DETECMANAGER  UNIT 


unit  DetecManager; 

interface 

uses 

CRT,  Intermanager, scoutmanager; 

f unct i on  Det ec t i onE f f or t ( i d , sd  : i nteger ; di st t s , dp: rea I ) : rea I ; 

f unc t i on  Absorpt i onCoef f i c i ent ( f : rea I ) : rea I ; 

function  UpdateNL<sc:byte): integer; 

funct i on  Scat ter ingAreaCtau, theta: real) : integer; 

function  Scatter i ngSt rengthCsvtprof i le: integer ;freq: real; winspeed, bottom: integer): integer; 
function  distancer(x1,y1,x2,y2:real):real; 
procedure  UpdateSonarParameters; 

implementation 

{  TOOLS} 

{ - - - ) 

function  distancer(xl.yl,x2.y2:real):real ; 

var  d1,d2:real; 

begin 

d1:=(x1-x2>*(x1-x2); 

d2:=(yl-y2>*(y1-y2); 

distancer:=d1+d2; 

end; 

function  power2(  base  :  real; 

exponent  :  real  )  :  real; 

Tbase  can  not  be  negative> 
begin 

power2;=:  exp(  exponent  *  InC  base  )  ); 
end;  {  function  power  } 


function  normalcdfWashb(mu,sig,x;real ) rreal ; 

Cbased  in  a  program  by  Alan  Washburn> 
const 

A=0.319381530;  B=-0.356563782;  C=1 .781477937;  D=-1 .821 255978; 

E=1. 330274429;  G=  0.231641900; 

var 

y,f,p  :real; 
positive  :boolean; 
begin 

if  sig<=0  then  sig:=1; 
x:=(x-mu)/sig; 
positive:=(x>=0); 
x:=abs(x); 

if  x>=100  then  x:=100; 

y:=1/Cl+G*x); 

f:=C+y*(D+y*E>; 

p:=1-exp(-0.5*x*x)/sqrt(8*arctan(1))*y*(A+y*(B+y*f)); 
if  positive  then  begin 
norma I cdf washb : =p; 
end  else  begin 

norma  I cdf washb : =1  - p; 
end; 
end; 


129 


function  sumcBCa , b: real ) : real ; 
begin 

sumclB:=(10/ln(10))*ln(  power2(10,a/10)  +  power2(10,b/10)  ); 

end; 


{SCATTERING  METHODS} 

{ - } 

function  TargetStrengthCtgx. tgy ,  sex.  scy . tgc : real ) :  i nteger ; 

const 

LOWTS=10; 

INCTS=8; 

var 

cospsisqr  :real; 

num,den  :real; 

tstemp  : integer; 

begin 

if  (sex  <>  tgx)  or  (scyo  tgy)  then  begin 

num:=sqr((scx-tgx>*sin<tgc*pi/l80)+(scy-tgy)*cos(tgc*pi/180)); 

den:=((scx-tgx)*(scx-tgx)+(scy-tgy)*(scy“tgy)); 

cosps i sqr : =num/den ; 

tstemp:=round<L0WTS+INCTS*(1 -cospsisqr)); 

T  argetSt rength : =tstemp; 
end  else  begin 
Ta rget St rength : =LOWTS; 
end; 
end; 


f uncti on  Scatteri ngArea (tau . theta : real ) ; i nteger ; 

const 

C=1500; 

begin 

theta:=theta*pi/180; 

ScatteringArea:=round(  (10/ln(10)  )  *  ln(tau*theta*C/2)  ); 
end; 


f  uncti  on  Surf  aceScatteri  ng  (f  r :  real ;  wi  nd :  i  nteger) :  i  nteger ; 

const 

THETA=37*pi/180; 

var 

h  :real; 

begin 

h:=power2(wind,2.5)*0.026; 

SurfaceScattering:=  round((l0/ln(10))  *  power2(ln(fr*h*sin(THETA)),  0.99)  -  45.3); 

end; 


f uncti on  BottomScatteri ng(f r : real ; bot : i nteger) : i nteger ; 
begin 

case  bot  of 

1  :  BottomScatteri ng;=round(  0.0120*fr*fr*fr  +  0-0283*fr*fr  +0-0062*fr  -  28.153); 

2  :  BottomScatteri ng:=round(  0.0076*fr*fr*fr  -  0.2337*fr*fr  +2.8700*fr  -  40.410); 

2  :  BottomScattering:=-18; 

end; 

end; 
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function  Scatteri ngStrength( svtprof i 1 e : i nteger ; f req : real ; wi nspeed . bottom : i nteger) : i nteger ; 

<svtprofile=1  (positive  velocity  gradient» 

“Csvtprof Ue=2  (negative  velocity  gradient)> 
begin 

case  svtprof ile  of 

1  :  ScatteringStrength:=:SurfaceScattering(freq,winspeecl); 

2  :  ScatteringStrength:=BottomScatteringCfreq,bottom>; 
end; 

end; 

procedure  UpdateSonarParameters; 
var 

irinteger; 

begin 

for  i:=1  to  numberofscouts  do  begin 
myscout [i] ■reverbarea:=ScatteringArea(myscout [i] .pulselength, 
myscout [i] .beamwidth); 

myscout  Ci] . SS:=ScatteringStrength (svtprof i letype, 

myscout  Ci 3 . f  requency , wi nd, bot tomtype) ; 

end; 

end; 


{NOISE  METHODS} 


-} 


function  AmbientNoiseCf ;rea1 ;w: integer): real ; 
begin 

AmbientNoise:=46. 121+2. 221*w-0.041*w*w-17-012*(  ln(f)/ln(10)  -  3  >; 
end; 


functi on  Sel f Noi se( v : real ; f : real ) : real ; 
begin 

SelfNoise:=  23+2.02* v+( 20/1 n( 10))* I n(25/f); 
end; 

function  NoiseLevel (scoutvel :real ;freq;real ;wi ndspeed:i nteger) ; integer; 
var 

An  :real; 

Sn  :real; 

begin 

An : =Ambi ent No i se( f  req , wi ndspeed) ; 

Sn: =Se I f Noi se( scoutve I , f  req) ; 

Noi seLevel :=round(sumdB(An, Sn) ) ; 
end; 


function  UpdateNKsc : byte) : integer; 
begin 

with  myscout Csc]  do  begin 

NL :=  Noi seLeve I (speed, frequency, wi nd) ; 
end; 
end; 


131 


{PROPAGATION  METHODS} 


■} 


function  AbsorptionCoeff i cient(f : real ) : real ; 
begin 

AbsorptionCoefficient:=  Se-S  /(  (0.7  /f*f)+1  ) 

+  0.04/«6000/f*f)+1)  +  4e-7  *  f*f  ; 

end; 


function  TransLossesGeorn{d,H: real ) : real ; 

var 


TLG  :real; 
begin 

if  d  <=  0  then  begin 
TLG:=1; 

end  else  begin 

TLG:=(10/LnC10))*(ln(d)+ln(H/pi)); 

end; 

T  ransLossesGeom:=TLG; 
end; 


{DETECTION  MODEL  METHODS} 

{ - — - } 

function  Detect i onEf fort (i d . sd: i nteger ; di stts . dp: real ) : real ; 
var 

p  ireal; 

TL  :real; 

TLGEOM  :real; 

FOM  :real; 

SE  :real; 

TS  ireal; 

SErevbr  :real; 

SEnoise  ireal; 

tx  : integer; 

ty  : integer ; 

sx  : integer; 

sy  : integer ; 

tc  :real; 

begin 

tx:=mytargetCid3 -x; 
ty:=mytarget[id]  .y; 
sx:=myscout [sd] .x; 
sy :  =inyscout  [sd] .  y; 
tc:=mytarget[id] .course; 

distts:=185*distts;C185  accounts  for  the  scale  factor> 

TL6eom:=  TransLossesGeom(distts,dp); 

TL : =TLGeom+di st ts*myscout [sd] *  abscoef f ; 

TS : =  TargetSt  rength ( tx, ty, sx, sy, tc) ; 

FOM:=  myscout  [sd]  .SL-myscout  [sd]  .NL-myscout  [sd]  .DT+TS+myscout  [sd]  .DI; 

SEnoise:=F0M-2*TL; 

SErevbr:=TS-TLGeom-myscout[sd]  -DT-myscout  [sd]  .reverbarea- myscout  [sd]  .SS; 
if  SErevbr  <  SEnoise  then  begin 
SE:=  SErevbr; 
end  else  begin 
SE:=  SEnoise; 
end; 
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p:=pouer2(1-normalcdfuashb(0,  NOISESIGMA.SE),  NUHOFPULSES); 
Detect i onEf f ort : =1 -p; 


begin 

end. 
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APPENDIX  G.  CELLSMANAGER  UNIT 


unit  CellsManager; 

interface 

uses 

GRAPH,  CRT  ; 


type 

probcelltype  =  array II.. 80,1.. 80]  of  single;. 


var 

probincel I  iprobcelltype; 

maxprob  :real; 

numberout  : integer; 

procedure  Update! nCe I I Pr obabi I i t i es; 
procedure  UpdateOutAreaProbabi lities; 
procedure  Wr i teProbi nCel I ; 


implementation 

uses  GraphManager,  TargtmanagerJnterManager; 


procedure  Clearance! Is; 
var 

i,j  : integer; 

begin 

for  i:=1  to  MAXGRID  do 
begin 

for  i:=1  to  MAXGRID  do 
begin 

probincelLIi, j] :=0; 

end; 

end; 

end; 


procedure  Update! nCel 1 Probabi 1 i ti es ; 

var 

i  : integer; 
begin 

maxprob; =0; 
numberout :=0; 

ClearallCells; 
for  i:=1  to  NUMBEROFTGTS  do 
begin 

if((mytargetCi] .X  >  GRIDSTEP)  and  Cmytarget li] .y  >  GRIDSTEP  ))  and 
((mytargetCi] .X  <MAXGRID*GRIDSTEP)  and  (mytarget li] .y  <MAXGRID*GRIDSTEP  )) 
then 
begin 

probincell [mytargetCi] -cellx, mytarget [i] .cel ly] := 
probincell Cmytarget [i] -cel lx,  mytargetCi] .celly]+mytargetCi] .prob; 
if  maxprob  <  probincell [mytargetCi] -cel lx,mytargetCf] -cel ly]  then 
begin 
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inaxprob:=  probincell  Dnytarget Ci] .cellx^mytarget Ci3 -celly] 
end; 

end  else 
begin 

inc(numberout); 

end; 

end; 

end: 


procedure  BroadcastProbinCell : 

var 

ij  :  integer; 

begin 

ass  i  gn<  outdata  I  i  nk/*c :  \  I  nk2  -  dat  ** ) ; 
for  i:=:1  to  MAXGRID  do 
for  j:=1  to  MAXGRID  do 

writeLn(outdatalink,probincell [i J] ); 

end; 

end; 

end; 


begin 

end. 
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APPENDIX  H.GRAPHMANAGER  UNIT 


unit  GraphManager; 

interface 

uses 

GRAPH,  CRT; 

const 

BOXSI2E=100; 

SCALECOLORS=10; 

BGIdir='-.\BGI'; 

GRIDC0L0R=32; 

BASEPATTERN=1; 

TEXTC0L0R=8; 

TARGETSCOLOR=30; 

SCOUTSCOLOR=24; 

PALETTE1=11; 

PALETTE2=12; 

MAXGRID=:80; 

GRIDSTEP=5; 

CENTER=MAXGRID*6RIDSTEP  div  2+6RIDSTEP; 

STEP  =  5; 

FRAME  =80; 

type 

colortype  =  arrayC1..SCALEC0L0RS]of  byte; 

cel  I type  =  array Cl.. MAXGRID,1..MAXGRID]  of  integer; 

const 

COLORCODE  :  colortype  =  (8,1,33,12,4,36,52,62,55,63); 


procedure  InitGraphics; 
procedure  SetGrid; 
procedure  paint(x,y,pat,col:word); 
procedure  SetScale; 
procedure  SetTimeCt: integer); 
procedure  DisplayTargets; 
procedure  DisplayScouts; 
procedure  PleaseWait; 
procedure  SetProbDistrCoLors; 
procedure  OutputGraph; 
procedure  Start; 
procedure  DisplayData; 

procedure  D  i  spl ayContac tSurmiary( srch r :  byte) ; 


implementation 

uses  TargtManager , ScoutManager ,Cel 1 sManager, Intermanager; 


procedure  InitGraphics; 
var 

GraphDriver, 

GraphMode, 

i  :  integer; 
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begin 

GraphD river  :=  GRAPH. Detect; 

GRAPH. Ini tGraph(GraphDriver,GraphMode,BGIdi r); 
for  i:=1  to  SCALECOLORS  do  begin 

SetPalette(i,COLORCODECi] ); 

end; 

end;{initgraphics} 


procedure  SetGrid; 
var 

i  : integer; 

oldstyle  :TextSettingsType; 
begin 

GetTextSett i ngs<  oldstyle) ; 

Set Pa I et  teC 1 5 , GR 1 DCOLOR ) ; 
for  i:=1  to  MAXGRID  do  begin 

line(i*GRIDSTEP,GRIDSTEP,i*GRIDSTEP,GRIDSTEP*MAXGRID); 

line(GRIDSTEP,i*GRIDSTEP,GRIDSTEP*MAXGRID,i*GRIDSTEP); 

end; 

SetTextStyle(2,0,<STEP+1)  div  2  ); 

line(GRIDSTEP+MAXGRID*GRIDSTEP  div  2,0,GRIDSTEP+MAXGRID*GRIDSTEP  div  2,GRiDSTEP); 
OutTextXY(MAXGRID*GRIDSTEP  div  2,0,chr(48)); 

OutTextXY(0,MAXGRID*GRIDSTEP  div  2,chr(48)); 

OutTextXY(MAXGRID*GRIDSTEP  div  4,0,concat<chr<49),chr(48))); 
OutTextXY(0,MAXGRID*GRIDSTEP  div  4,concat(chr<49),chr(48))); 
0utTextXY(MAXGRID*GRIDSTEP*3  div  4,0,concat(chr(49),chr(48))); 
OutTextXY(0,MAXGRID*GRIDSTEP*3  div  4,concat(chr(49),chr(48»); 
0utTextXY(MAXGRID*GRIDSTEP-GRIDSTEP,0,concat(chr(50),chr(48))); 
OutTextXY(0,0,concat(chr<50),chr(48))); 

OutTextXY(0,MAXGRID*GRIDSTEP-GRIDSTEP,concat(chr(50),chr(48))); 

SetColordS); 

end; 


procedure  paint(x,y.pat,col  rworcf); 
begin 

SetFillStyleCpat,  col); 

FloodFilUx,y,  GetMaxColor); 
end; 

procedure  SetScale; 

var 

i  : integer; 

outstring: string; 
begin 

I ineCSTEP, FRAME*STEP+STEP*3, 

STEP*2*SCALEC0L0RS+STEP,FRAME*STEP+STEP*3); 

I ineCSTEP, FRAME*STEP+STEP*5, 

STEP*2*SCALEC0L0RS+STEP, FRAME*STEP+STEP*5 ) ; 
for  i:=1  to  SCALECOLORS+1  do 
begin 

I ine( i*2*STEP-STEP, FRAME*STEP+STEP*3, 
i *2*STEP - STEP , FRAME*STEP+STEP*5 ) ; 

end; 

for  i:=1  to  SCALECOLORS  do 
begin 

pai nt( i*2*STEP, FRAME*STEP+STEP*4,BASEPATTERN, i ); 
SetTextStyle(2,0,(STEP+2)  div  2  ); 

OutTextXY< i*2*STEP- 1 0, FRAME*STEP+STEP, 
concatCchr(46),chr(47+i))); 
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end; 

0utTextXY((SCALEC0L0RS+1)*2*STEP-STEP,FRAME*STEP+STEP,chr<49)); 

end; 


procedure  DisplayTargets; 

var 

1  : integer; 

begin 

Set Pa I et  t e( PALETTEl , TARGETSCOLOR ) ; 
for  i:=1  to  NUMBEROFTGTS  do 
begin 

if  (<mytargetCi3.x  >GRIDSTEP)  and  (mytarget [i] -x  <MAXGRID*GRIDSTEP))  and 
<(mytargetEi]  .y  >6RIDSTEP)  and  (mytarget  [i]  .y  <MAXGRID*GRIDSTEP))  then 
PutPixeUinytarget[il.x,iiiytargetCi].y,PALETTE1); 

end; 

end; 


procedure  DisplayScouts; 

var 

i  rinteger; 

begin 

SetPa I et  t e( PALETTE2, SCOUTSCOLOR ) ; 
for  i:=1  to  numberofscouts  do 
begin 

if  ((rayscouttij.x  >GRIDSTEP)  and  (myscoutti]  .x  <MAXGRIO*GRIDSTEP»  and 
((myscoutCil.y  >GRIDSTEP>  and  (myscoutti]  .y  <MAXGRID*GRIDSTEP))  then 
begin 

PutPixel (myscoutti]  .x^myscoutCi]  .y,PALETTE2); 

SetColor(PALETTE2); 

OutTextXY(myscoutCi].x  +  5,  myscoutti] .y  -  5,  concat<inyscout ti] -id 
,chr(48),chr(48),chr(48+i))); 

Circle(inyscoutti].x,myscoutti]  -y,myscout ti]  .maxscale); 
if  myscoutti]. I S_H0T  then  begin 

PutPixeKmyscoutti]  .subx+5,  myscoutti]  .suby+5,8); 

PutPixeUmyscoutti] .subx+6,  myscoutti] .subyM-4,8); 

PutPixeKmyscoutti] .subx+7,  myscoutti] .suby+3,8); 

PutPixeKmyscoutti]  .subx+8,  myscoutti].suby+2,8); 

PutPixeKmyscoutti]  .subx+9,  myscoutti]  .suby+1,8); 

PutPixeKmyscout  ti]  -subx+10,  myscout  ti]  .suby+0,8); 

PutPixeKmyscoutti]  .subx+11,  myscoutti]  -suby-1,8); 

PutPixeKmyscoutti]  .subx+4,  myscoutti]  .suby+4,8); 

PutPixeKmyscoutti] .subx+3,  myscoutti].suby+3,8); 

PutPixeKmyscoutti]  -subx+2,  myscoutti]  -suby<-2,8); 

PutPixeKmyscoutti]  -subx+1,  myscoutti]  .suby+1,8); 

PutPixeKmyscout  ti]  .subx+0,  myscoutti3-suby+0,8); 

PutPixeKmyscoutti]  .subx-1,  myscoutti]  .suby-1,8); 

SetColorC14); 

OutTextXY(myscoutti].subx  +5,  myscoutti] -suby- 10,  concat<chr(48+i))>; 
SetColordS); 

end; 

end; 

end; 

end; 


procedure  SetProbDistrColors; 

var 

i,j,c  : integer; 
begin 
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if  maxprob  >  0  then  begin 

for  i:=1  to  MAXGRID-1  do  begin 
for  j:=1  to  MAXGRID-1  do  begin 

for  c:=0  to  SCALECOLORS-1  do  begin 
if  <100*probincel I Ci , j]/maxprob>c*SCALECOLORS) 

and  (100*probinceU[i,j]/maxprob<=<c*SCALECOL0RS+100/SCALEC0LORS))  then  begin 
paint(i*GRIDSTEP+1,j*GRIDSTEP+1,BASEPATTERN,c+1); 
end; 

end; 

end; 

end; 

end  else  begin 

OutTextXY(40,40/No  target  in  the  area, it  is  gone'); 
end; 
end: 


procedure  PleaseWait; 
begin 

InitGraphics; 

outtextxy(MAX6RID*6RID$TEP+10,  GRIDSTEP,  'Please  wait..'); 

end; 


procedure  DisplayLogo; 
var 

i  : integer; 

MAXx  : integer; 

MAXy  : integer; 

begin 

SetPalette(15,GRIDCOLOR); 

MAXx:=GetMaxX; 

MAXy:=GetMaxY; 

I i ne( r ound< ( MAXx/2 ) - BOXS I ZE ) , roundC ( MAXy/2 ) - BOXS I ZE ) , 
roundC (MAXx/2)+BOXSI ZE) , roundC (MAXy/2) -BOXSIZE )); 

I i ne( round( <MAXx/2) -BOXS I ZE) , roundC (MAXy/2)+B0XS I ZE ) , 
roundC ( MAXx/2 )+BOXS I ZE ) , roundC  CMAXy/2 )+BOXS IZE)); 

I 1 neC roundC C MAXx/2 )- BOXS I ZE ), roundC CMAXy/2 )- BOXS I ZE ) , 
roundCCMAXx/2)-B0XSIZE),roundCCHAXy/2)+B0XSI2E)); 

I i ne  C  roundC  C  MAXx/2 )+BOXS IZE), roundC  C MAXy/2 ) - BOXS IZE), 
roundC  CMAXx/2)+BOXS IZE), roundC  CMAXy/2 )+BOXS I ZE ) ) ; 

1 1 neC  roundC  CMAXx/2 ) -BOXSI ZE ) , roundC  CMAXy/2 )+BOXS I ZE+5 ) , 
roundC  C MAXx/2 )+BOXS IZE ) , roundC  CMAXy/2)+BOXSIZE+5) ); 

I i neC  round C  C MAXx/2 ) - BOXS IZE), roundC  CMAXy/2 )+BOXS 1 2E+25 ) , 
roundCCMAXx/2)+B0XSIZE),roundCCMAXy/2)+B0XSI2E+25)); 
for  i:=0  to  SCALECOLORS  do 
begin 

lineCroundCCMAXx/2)-BOXSI2E)+i*roundCBOXSIZE/5),roundCCMAXy/2)+BOXSIZE+5), 
roundC  C  MAXx/2  )  -  BOXS  I  ZE  )+ i  *roundC  BOXS  1 2E/5  ) ,  roundC  C  MAXy/2  )+BOXS  I ZE+25  )  )  ; 

end; 

for  i:=1  to  SCALECOLORS-1  do 
begin 

paintCroundCCMAXx/2)-BOXSIZE)+1+i*roundCBOXSIZE/5),roundCCMAXy/2)+BOXSIZE+6), 

BASEPATTERM+8,i); 

delayC160); 

end; 

SetPaLetteCl,8); 

pai nt C  roundCMAXx/2) , roundCMAXy/2) , BASEPATTERN , 1 ) ; 

SettextSty leC 1,0,1); 

OuttextXYCroundCMAXx/2)-50, roundCMAXy/2)-50, '  NRTCAS  '); 

SettextStyleC2,0,4); 

0uttextXYCroundCMAXx/2)-100,roundCMAXy/2)“25,  '  Reactive  Target'); 
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Out textXYC  round(MAXx/2) - 1 00, round(MAXy/2 ) - 1 5 
SettextStyle(2,0,4); 

OuttextXY( round(MAXx/2)- 1 00 , round(MAXy/2) , 
OuttextXY(round(MAXx/2)- 100, round(MAXy/2)+35, 
SettextStyle(2,0,3); 

OuttextXYC  round(MAXx/2) -100, round(MAXy/2 )+70 , 
end; 


Computer  Assisted  Passive  Search'} 

'  Release  1.1  '); 

'  NPS  '); 

'  C.RECALDE  '); 


procedure  Di spl aylntro ; 
var 

MAXx  : integer; 

MAXy  : integer; 

begin 

cleardevice; 

SetPa I et te( 1 5 , GRIDCOLOR ) ; 

MAXx:=:GetMaxX; 

MAXy:=GetMaxY; 

line(GRIDSTEP,GRIDSTEP,MAXx-GRIDSTEP,GRIDSTEP); 

line<GRIDSTEP,GRIDSTEP,GRIDSTEP,MAXy-GRIDSTEP); 

line(MAXx-GRIDSTEP,MAXy-GRIDSTEP,GRIDSTEP,MAXy-GRIDSTEP); 

line(MAXx-GRIDSTEP,MAXy-GRIDSTEP,MAXx-GRIDSTEP,GRIDSTEP); 

SetColor(TEXTCOLOR-2); 

SetTex tS ty I e( 1 , 0 , 1 ) ; 

SetT  ext Just i f y( 1,1); 
outtextxy(round(MAXx/2),GRlDSTEP+30, 

'  Non  Reactive  Target  Computer  Assisted  Search  '); 

SetTextStyle(2,0,4); 

outtextxy(round(MAXx/2),GRIDSTEP+90, 

'  This  program  will  perform  a  Montecarlo  Simulation  to  assist  in  '); 
outtextxyC  round(MAXx/2) , GRIDSTEP+1 00, 

'  conduction  of  a  Lost  Contact  Search  of  a  Non  Reactive  Target  by  '); 
outtextxyC  round(MAXx/2) , GRIDSTEP+1 10, 

'  Passive  Searchers'); 

outtextxyC round(MAXx/2),GRIDSTEP+130, 

'  The  user  must  provide  parameters  regarding:  '); 

outtextxyC  r oundC  MAXx/2 ) , GR I DSTE  P+1 50 , 

'  THE  TARGET  '); 

outtextxyC roundCMAXx/2),6RIDSTEP+160, 

'  THE  SEARCHER  '); 

outtextxyC  roundCMAXx/2) , GRIDSTEP+1 70, 

'  THE  SIMULATION  '); 

outtextxyC  roundCMAXx/2 ) , GR IDSTEP+21 0, 

'  The  program  will  output  a  Probability  Distribution  Map  in  Color' ); 
outtextxyC  roundCMAXx/2) , GRIDSTEP+220, 

'  Codes  relative  to  the  maximum  value  of  probability  encountered'); 
outtextxyC  roundCMAXx/2 ) , GR  IDSTEP+230 , 

'  in  any  of  the  square  cells  in  wich  the  Local  Area  is  divided.  '); 
outtextxyC  roundC  MAXx/2 ) , GR I D  STEP+300 , 

'  Please,  press  any  key  to  continue'); 

readkey; 

end; 


procedure  Start; 

begin 

clrscr; 

InitGraphics; 

DisplayLogo; 

<:delayC4000);> 

Displaylntro; 
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Closegraph; 

end; 


procedure  SetTime(t: integer); 

var 

fstdigit  :byte; 

snddigit  :byte; 

trddigit  :byte; 

begin 
t:=t*5; 

SetTextStyle(2,0,5); 

0utTextXY(STEP,FRAME*STEP+STEP*5,  'time  :  '); 

SetColor(TEXTC0L0R-2); 
fstdigit:=round<t  div  100); 
snddigit:=round((t  mod  100)  div  10)  ; 
trddigit:=round((t  mod  100)  mod  10)  ; 

OutTextXY(STEP+STEP*1 1 , FRAME*STEP+STEP*5, 

concat<chr(48+fstdigit),chr(48+snddigit),chr<48+trddigit))); 

end; 


procedure  OutputResults; 

var 

: integer; 

: integer; 

:byte; 


MAXx 
MAXy 
fstdigit 
snddigit 
trddigit 
fthdigit 
fihdigit 
P 

count 

holder 

stringholder 


:byte; 

:byte; 

:byte; 

:byte; 

:real; 

:byte; 

: integer; 
istring; 


begin 

MAXx:=GetMaxX; 

MAXy:=GetMaxY; 

SetPa I et  t e( 1 5 , GR I DCOLOR ) ; 

Setcolor(15); 

I i neC  FRAME*STEP+5 , STEP , MAXx- 1 0, STEP) ; 

I i ne<  FRAME*STEP+5 , STEP , FRAHE*STEP+5 , MAXy- 10); 

line(MAXx-10,MAXy-10,STEP*FRAME+5,MAXy-10); 

line(MAXx-10,MAXy-10,MAXx-10,STEP); 

SetColor<TEXTC0L0R-2); 

OutTextXY(FRAME*STEP+15,STEP*3/  RTCAS  STATUS  '); 
SetColordS)  ; 

SetTextStyle(2,0,4); 

OutTextXY(FRAME*STEP+15,STEP*28/UNDET  PROS  MASS  '); 
p:=undetectedmass* 10000; 
fstdigit:=round(p)  div  10000; 
snddigit:=(round(p)  mod  10000)div  1000; 
trddigit:=((round(p)  mod  10000)mod  1000)div  100; 
fthdigit:=(((round(p)  mod  10000)mod  1000)mod  100)div  10; 
fihdigit:=(((round(p)  mod  10000)mod  1000)mod  100)mod  10; 

Set  Co I or ( TEXTCOLOR - 2 ) ; 

OutTextXYC  FRAME*STEP+STEP*30, STEP*28, 

concat(chr(48+fstdigit),chr(46),chr<48+snddigi t), 
chr(48+trddigit),chr(48+fthdigit),chr<48+fihdigit))); 

SetColordS)  ; 

SetTextStyle(2,0,4); 
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OutTextXY(FRAME*STEP+1 5, STEP*32/ Number  of  TGTS  OFF  '); 
fstdigit:=round(numberout  div  100>; 
snddigit:=round((numberout  mod  100)  div  10)  ; 
trddigit:=round<(numberout  mod  100)  mod  10)  ; 

SetColor(TEXTCOLOR-2); 

OutTextXY(FRAME*STEP+STEP*30,STEP*32, 

concat<chr(48+fstdigit),chr(48+snddigit),chr(48+trddigit))); 

SetPa I et teC 1 5 , GRIDCOLOR ) ; 
setcolordS); 

I i ne( FRAME*STEP+5 , STEP*40 , MAXx- 1 0 , STEP*40 ) ; 

SetTextStyle(2,0,4); 

Set Co I or ( TEXTCOLOR - 2 ) ; 

0utTextXY(FRAME*STEP+15,STEP*42/  MC  SIMULATION  INFOO; 

SetColorCIS); 

0utTextXYCFRAME*STEP+15,STEP*45/TARGET  DATA' ); 
0utTextXY<FRAME*STEP+15,STEP*48, 'DATUM  ERROR  (nm)  '); 
fstdigit:=round(sigmarimt  div  100); 
snddigit:=round((sigmarinit  mod  100)  div  10)  ; 
trddigit:=round((sigmarinit  rood  100)  mod  10)  ; 

Set Co I or  C  TEXTCOLOR- 2 ) ; 

Out TextXY ( FRAME*STEP+STEP*30 , STEP*48 , 

concat(chr(48+fstdigit),chr(48+snddigit), 

chr(48+trddigit))); 

SetColordS); 

OutTextXY(FRAME*STEP+15,STEP*50, 'REACTION  RANGE  (nm)  '); 
fstdigit:=round(reactrange  div  100); 
snddigit:=:round((reactrange  mod  100)  div  10)  ; 
trddigit:=round((reactrange  mod  100)  mod  10)  ; 

Set Co I or  <  TEXTCOLOR- 2 ) ; 

OutTextXY(FRAME*STEP+STEP*30,STEP*50, 

concat(chr(48+fstdigit),chr(48+snddigit), 

chr(48+trddigit))); 


SetColorCIS); 

OutTextXY(FRAME*STEP+15,STEP*52, 'REACTION  TACK  '); 
fstdigit:=round<reactiontack  div  100); 
snddigit:=round(( react iontack  rood  100)  div  10)  ; 
trddigit:=roundC( react iontack  mod  100)  mod  10)  ; 

Set Co I or ( TEXTCOLOR- 2 ) ; 

Out TextXY ( FRAME*STEP+STEP*30 , STEP*52 , 

concat(chr(48+fstdigit),chr(48+snddigit), 

chr(48+trddigit))); 


SetcolordS); 

OutTextXYC  FRAME*STEP+1 5 , STEP*54 , ' M I N  SPEED  ' ) ; 

fstdigit:=round(minsubspeed  div  100); 
snddigit:=roundC(minsubspeed  mod  100)  div  10)  ; 
trddigit:=round((minsubspeed  mod  100)  mod  10)  ; 
SetCoLor<TEXTCOLOR-2); 

Out  TextXY ( FRAME*STEP+STEP*30 , STEP*54 , 

concat(chr(48+fstdigit),chr(48+snddigit), 

chr(48+trddigit))); 


SetcolordS); 

OutTextXY(FRAME*STEP+1S,STEP*S6, 'MAX  SPEED 
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fstdigit:=round(fnaxsubspeed  div  100); 
snddigit:=round((inaxsul:^peed  mod  100)  div  10)  ; 
trdd1git:=round((inaxsubspeed  mod  100)  mod  10)  ; 
SetColor(TEXTC0L0R-2); 
OutTextXY(FRAME*STEP+STEP*30,STEP*56, 

concat(chr(48+fstdigit),chr(48+snddigit), 

chrC48+trddigit))); 


SetColordS); 

I i ne( FRAME*STEP+5 , STEP*60 , MAXx- 1 0 , STEP*60 ) ; 

SetColor(TEXTC0L0R-2); 

OutTextXY(FRAME*STEP+15,STEP*62,'  REAL  WORLD  SEARCH  INFO  >); 
SetColorCIS); 

OutTextXY< FRAME*STEP+1 5, STEP*65, 'SEARCHERS  DATA' ); 
0utTextXY(FRAME*STEP+15,STEP*69, 'SCOUT  COURSE  SPEED  SONAR  HOT  '); 
for  count :=1  to  8  do  begin 

if  count  <=  numberof scouts  then  begin 
SetColor(TEXTCOLOR-2); 

OutTextXY(FRAME*STEP+15,STEP*70+count*2*STEP,  concatC 

myscouttcount] -id,chr(48),chr(48),chr<48+count),'  ')); 
fstdigit:=round(myscout [count] .course  div  100); 
snddigit:=roundC(myscout [count] .course  mod  100)  div  10)  ; 
trddigit::=round( (m/scout  [count]  .course  mod  100)  mod  10)  ; 
OutTextXY ( FRAME*STEP+55 , STEP*70+count*2*STEP , conca t ( 

chr(48+fstdigit),chr(48+snddigit),chr(48+trddigit))); 
holder :=round(myscout [count] .speed); 
fstdigit:=round(holder  div  100); 
snddigit:=roundC (holder  mod  100)  div  10)  ; 
trddigit:=round( (holder  mod  100)  mod  10)  ; 

CXitTextXY( FRAME*STEP+95 , STEP*70+count*2*STEP, concat ( 

ch r ( 48+f s tdi g i t ) , ch r (48+snddi g i t ) , ch r ( 48+ t rdd i g i t ) ) ) ; 
if  myscout [count] .SONAR_ON  then  begin 
stringholder:-'  ON  '; 
end  else  begin 

stringholder:='  OFF'; 
end; 

OutTextXY(FRAME*STEP+135,STEP*70+count*2*STEP,stringholder); 
<***>  myscoutC2].IS_HOT:=TRUE; 

if  myscout [count] . I S_H0T  then  begin 
stringholder:='  YES'; 

SetColor(14); 

OutTextXY(FRAME*STEP+175,STEP*70+count*2*STEP,stringholder); 

end; 

if  myscout [count] .ini tialhot  <>  0  then  begin 
D  i  splayContactSiJtinary(  count )  ; 
end; 
end; 
end; 

SetTextStyle(2,0,4); 

SetColordS); 

OutTextXY(FRAME*STEP+15,MAXy-30, 'Hit  any  key  to  go  back  to  main  menu'); 

end; {procedure} 


procedure  Displ  ayContactSunimary(srchr :  byte) ; 

var 

holder: integer; 
f stdi gi t : i nteger ; 
snddi gi t : i nteger; 
t rddi gi t : i nteger; 
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begin 

SetColor(TEXTC0L0R-2); 

OutTextXY(30*STEP,FRAME*STEP+STEP/ INITIAL  CONTACT  REPORTS  '); 
SetCo  lords  ); 

OutTextXY(30*STEP,FRAME*STEP+STEP*3/SEAR  TIME  COURSE  SPEED  ') 
SetColor(14); 

0utTextXY(30*STEP,FRAME*STEP+STEP*3+STEP*srchr*2,concat( 

myscout Esrchr].id,chr(48>,chr(48),chr(48+srchr))); 
fstdigit:=round(myscoutCsrchr]  .imtialhot  div  100); 
snddigit:=round((inyscout[srchr].initialhot  mod  100)  div  10)  ; 
trddigft:=round<<ii!yscout[srchr] . imtialhot  mod  100)  mod  10)  ; 
0utTextXY(30*STEP+35,FRAME*STEP+STEP*3+STEP*srchr*2,concat( 

chr(48+fstdigit),chr(48+snddigit),chr(48+trddigit))); 
fstdigit:=round(myscoutCsrchr] -subcourse  div  100); 
snddigit:=round((myscout[srchr] .subcourse  mod  100)  div  10)  ; 
trddigit:=round((inyscoutCsrchr)  .subcourse  mod  100)  mod  10)  ; 
OutTextXY(30*STEP+80, FRAME*STEP+STEP*3+STEP’^srchr*2, concat( 

chr(48+fstdigit),chr<48+snddigit),chr(48+trddigit))); 
ho I der : = round( myscout [s  rch  r] . subspeed) ; 
fstdigit:=round(holder  div  100); 
snddigit:=round( (holder  mod  100)  div  10)  ; 
trddigit:=round((holder  mod  100)  mod  10)  ; 
OutTextXY(30*STEP+125,FRAME*STEP+STEP*3+STEP*srchr*2,concat( 
chr(48+fstdigit),chr(48+snddigit),chr<48+trddigit))); 
end; 


procedure  DisplayData; 
begin 

ClearDevice; 

SetGrid; 

SetScale; 

DisplayTargets; 

DisplayScouts; 

undetectedmass:=1 ; 

OutputResults; 

readkey; 

cleardevice; 

CloseGraph; 

end; 


procedure  OutputGraph; 
begin 

ClearDevice; 

SetGrid; 

SetScale; 

SetTime(totalsimtime); 

CDispLayTargets;> 

SetPr obD i st  rCo I ors ; 

DisplayScouts; 

OutputResults; 

readkey; 

ClearDevice; 

CloseGraph; 

end; 

begin 

end. 
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